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ELECTRICAL  BREAKDOWN  BETWEEN  CYLINDRICAL  ELECTRODES  IN  A  PLASMA 


ABSTRACT 


A  breakdown  voltage  of  as  low  as  40  volts  has  been  observed  between 
cylindrical  electrodes  immersed  in  a  flowing  argon  plasma  which 
indicates  that  some  mechanism  other  than  the  usual  spark  discharge 
mechanism  is  present.  There  is  strong  evidence  that  at  plasma  densities 
in  the  range  10^  to  10^  cm  ^  the  breakdown  is  triggered  by  electrons 
ejected  from  the  cathode  due  to  the  intense  electric  field  in  the  cathode 
sheath.  Another  is  the  high  concentration  of  neutral  particles  in  the 
boundary  layer  around  the  electrode  surface.  The  exact  theory  of  this 
phenomenon  is  presented  and  shows  good  agreement  with  experimental 
measurements . 

It  is  hoped  that  the  results  presented  in  this  thesis  will  be 
applicable  to  the  reignition  mechanism  after  current  zero  of  a.c. 
circuit  breakers  and  to  the  triggering  mechanisms  of  unipolar  arcs, 
vacuum  gaps  and  short  gaps . 

During  the  course  of  this  research  a  new  technique  of  measuring 
the  electrical  conductivity  as  well  as  new  theories  for  the  variation 
of  the  electromagnetically  driven  plasma  velocity  and  for  the  ion 
current  collection  of  a  probe  aligned  parallel  to  the  collision 
dominated  plasma  flow  were  developed. 
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CHAPTER  I  INTRODUCTION 


In  the  course  of  plasma  density  measurements  with  floating 
double  Langmuir  probes  in  a  high  density  plasma,  it  was  found  that 
at  an  applied  voltage  between  the  probes  greater  than  a  certain 
critical  value  depending  on  the  plasma  density,  a  sudden  increase 
in  the  probe  current, which  cannot  be  explained  by  the  existing  probe 
theories, was  observed.  From  its  current -voltage  characteristics  it 
was  deduced  that  this  current  increase  was  due  to  electrical 
breakdown  between  probes.  The  breakdown  voltage,  however,  was  very 
different  from  that  predicted  by  the  well  known  Paschen's  Law  (as  low 
as  40  volts  was  observed)  ,  which  indicated  that  some  mechanism  other 
than  the  Paschen  breakdown  mechanism  was  responsible.  The  same  low 
voltage  breakdown  phenomena  in  high  density  plasmas  have  also  frequently 
been  observed  at  many  other  laboratories. 

We  find  a  similar  and  a  very  important  example  of  this  breakdown 
phenomena  in  a  plasma  in  a  reignition  after  current  zero  of  alternating 
current  circuit  interrupting  devices.  During  the  post  arc  period,  i.e. 
during  a  few  microseconds  after  current  zero  the  gas  conductivity 
between  the  electrodes  is  quite  high  due  to  the  heat  generated  by 
the  previous  arc.  If  the  ionization  of  the  interelectrode  gas  is 
sufficient  the  recovery  voltage  may  break  down  the  gap  and  re-establish 
the  arc . 

Another  important  example  occurs  in  a  high  density,  high 


2 


t empera ture  plasma,  in  contact  with  an  insulated  metal  wall.  As  the 
temperature  increases^ the  floating  potential  of  the  wall  becomes 
^  gti  ©  it  (negatively)  than  the  ionization  potential  of  the  ambient  gas 
and  the  unipolar  arc  (arc  which  joins  a  metallic  electrode  (cathode) 
to  a  plasma  (anode))  can  be  established  (1.1). 

Electrical  breakdown  of  vacuum  insulation  is  a  problem  of 
considerable  practical  importance  owing  to  the  increasingly  widespread 
application  of  so-called  vacuum  insulation  in  high  voltage  devices 
such  as  particle  accelerators  ,  x-ray  tubes  ,  high  frequency  low-loss 
condensers,  etc.  Field  emission  and  charged  particle  accumulation  in 
front  of  the  cathode  are  often  responsible  for  this  kind  of  breakdown 
and  so  the  breakdown  mechanism  for  this  case  is  similar  to  the  one 
discussed  here.  (In  the  case  of  vacuum  breakdown  the  neutral  gas  is 
generated  by  the  vaporization  of  anode  material.) 

The  importance  of  a  study  of  the  electrical  strength  of  gases 
with  high  ionization  and  high  temperature  is,  therefore,  clear. 
Nevertheless  the  works  of  Sharbaugh  et.al.(1.2)  and  Lee  et.al.  (1.3) 
seem  to  be  the  only  existing  studies  on  the  breakdown  phenomena  in 
the  plasma.  In  their  studies,  however,  the  range  of  the  charged 
particle  density  was  not  wide  (0.1-4x10^  cm  )}  the  neutral  density 
not  known,  and  so  an  exact  theory  was  not  given.  Thus  it  is  an  urgent 
requirement  to  measure  the  breakdown  voltages  as  a  function  of  the 
plasma  density  and  also  to  develop  a  theoretical  explanation. 

As  is  usual  in  any  electrical  breakdown  experiment  ,  there  are 
many  factors  to  be  taken  into  account  (e.g.  surface  conditions  of  the 
electrodes,  gas -electrode  interactions,  charged  particle  effects, 
etc.)  for  the  development  of  a  breakdown  theory.  Consequently  this 
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thesis  deals  with  the  aspects  of  surface  chemistry  to  solid  state 
physics,  as  well  as  plasma  physics  and  gas  dynamics,  including  the 
boundary  layer  problems. 

An  electromagnetically  driven  shock  tube  was  used  for  the 
production  of  the  sample  plasma.  The  plasma  thus  produced  had  many 
desirable  properties: 

(1)  A  very  wide  range  of  plasma  density  can  be  provided  (ne=10^  - 
10^6  cm“3).  Below  ne=10^  cm"^,  the  breakdown  mechanism  approaches 
that  of  Paschen  breakdown;  i.e.,  Crowe's  theory  (1.4)  on  the  effect 
of  space  charge  on  the  breakdown  voltage  becomes  applicable.  The 
mechanism  in  this  density  region  is  therefore  already  known. 

(2)  A  wide  range  of  combinations  of  plasma  and  neutral  density  can 
be  obtained.  Since  the  plasma  density  and  the  neutral  density  are 
the  two  ruling  factors,  this  option  is  important. 

(3)  Electrons,  ions  and  neutrals  are  in  thermal  equilibrium.  They 
are  thermalized  at  the  first  stage  of  their  acceleration  in  the  shock 
driver  and  when  being  tested,  they  are  in  an  "after  glow"  period  when 
a  Maxwellian  velocity  distribution  can  be  expected. 

(4)  Rankine-Hugoniot  relations  can  be  used  under  some  restrictions 
to  evaluate  the  densities  and  temperatures  of  both  ionized  and 
neutral  particles.  This  simple  and  reliable  method  of  determining 
especially  the  neutral  density  cuts  the  experimental  effort  considerably. 

(5)  The  plasma  has  a  well  defined  flow  velocity  with  respect  to  the 
electrode  system.  In  many  plasmas  where  electrical  breakdown  occurs, 
such  as  the  plasma  in  the  circuit  breaker,  flash  over  in  an  electrical 
motor,  theta-pinched  plasma,  plasma  created  by  plasma  gun,  there  is 

a  motion  of  ionized  and  neutral  gas  relative  to  the  electrode  system 
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and  this  flow  effect  would  have  a  considerable  influence  on  the 
breakdown  voltage.  Thus  this  property  is  important. 

The  contents  of  this  book  are  as  follows: 

Chapter  II  The  theories  of  plasma  properties  and  electrostatic 
probes  used  in  the  later  part  of  this  thesis  are  briefly  summarized. 
The  first  part  of  the  chapter  is  mainly  concerned  with  the  transport 
phenomena,  and  the  second  part  with  ion  sheath  current  theory  and  with 
the  modified  floating  double  probe  theory  which  applies  to  our 
experimental  situations. 

Chapter  III  The  apparatus  used  in  this  experiment  is  explained 
and  the  diagnostic  techniques  used  are  described.  In  section  3.2.5 
a  newly  developed  technique  of  measuring  the  plasma  conductivity  is 
discussed  in  detail. 

Chapter  IV  The  plasma  used  in  this  experiment  is  examined.  In 
the  first  half  a  new  model  of  the  electromagnetically  driven  shock 
tube  is  used  to  derive  a  theory  for  the  plasma  velocity  in  the  shock 
tube.  In  the  second  half  the  experimentally  determined  plasma 
parameters  using  many  diagnostic  techniques  are  tabulated. 

Chapter  V  Many  models  of  the  sheath  and  electrode  configuration 
are  examined  one  by  one  to  provide  a  solid  ground  for  determining  the 
most  appropriate  theory  to  account  for  the  breakdown  mechanism  between 
electrodes  in  the  plasma. 

Chapter  VI  The  various  breakdown  measurements  are  described. 

Chapter  VII  The  model  of  the  breakdown  mechanism  based  on  the 
theories  in  Chapter  V  is  developed  and  compared  with  the  results  in 
Chapter  VI.  This  is  also  compared  with  existing  theories  and 
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experimental  results. 

The  numerical  solutions  which  appear  frequently  in  this  thesis 
were  carried  out  using  the  IBM  360  computer  at  the  University  of 
Alberta  with  either  Fortran  IV  or  APL  languages. 


(1.1)  A.E.  Robson  and  P.C.  Thonemann:  Proc .  Phys .  Soc .  13_  508  (1959) 
D.W.  George  and  P.H.  Richards:  Brit.  J.  Appl.  Phys.  1_  1171  (1968) 

(1.2)  A.H.  Sharbaugh ,  P.K.  Watson,  D.R.  White,  T.H.  Lee  &  A.  Greenwood: 
A.I.E.E.  Trans.  Pt.  Ill,  P.A.S.  80  333  (1961) 

(1.3)  T.H.  Lee,  A.N.  Greenwood  &  D.R.  White:  I.E.E.E.  Trans.  P.A.S.  84_ 
1116  (1965) 

(1.4)  R.W.  Crowe,  J.K.  Bragg  &  V.G.  Thomas:  Phys.  Rev.  96_  10  (1954) 
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CHAPTER  II  PLASMA  PROPERTIES  AND  ELECTROSTATIC  PROBES 


2.1  ELECTRICAL  CONDUCTIVITY  OF  PLASMA 


2.1.1  Slightly  Ionized  Gas 


Consider  a  Lorentzian  gas  where  ions  are  neglected.  Minority 
particles  (electrons)  lose  momentum  by  colliding  with  the  majority 
particles  (neutrals)  having  a  collision  frequency  of  y  .  If  can  be 
considered  to  be  independent  of  electron  energy,  the  mobility  (ju)  and 
diffusion  constant  (D)  can  be  derived  from  the  Boltzmann  equation  as, 


I1- 


and 


D= 


v‘ 


3  V, 


(2.1) 


\rr\ 


If  the  electrons  have  a  Maxwellian  velocity  distribution,  we  find  that 


the  mean  free  path  of  electrons  in  neutrals  can  be  expressed  as. 


V  _] _  v  1  /  8-feTk  (2.2) 

n0  Qe  y rr\  ym  J  TC  me  * 

where  Qe  and  nQ  are  the  collision  cross  section  of  the  neutrals  for 
scattering  electrons  and  the  neutral  molecule  number  density,  respectively. 
Therefore  Eq .  2.1  becomes, 


_  jw  e _ 

“  ~  noQe 

D_  ftT  = 

/8fig  n0Qe 


(2.3) 

(2.4) 


Now  the  electrical  conductivity  takes  the  form,  <r=nee)U.  This 
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was  calculated  for  the  case  of  a  slightly  ionized  gas  of  rigid  elastic 
spherical  molecules  by  Chapman  and  Cowling  (2.1)  more  rigorously  with 
their  results  that, 


„-_il32*3/rc  ne  e‘  _  ,-a  ne  /(  ,x 

jf  "° QeJ% 


(2.5) 


The  values  of  Qg  for  argon  were  obtained  by  Townsend  and  Bailey  (2.2) 
and  shown  in  Fig.  2.1. 


2.1.2  Fully  Ionized  Gas 


If  an  electron  with  velocity  v  encounters  an  ion,  the  electron 
is  scattered  by  the  Coulomb  force.  The  collision  frequency  for  momentum 
transfer  is  obtained  by  using  the  laws  of  conservation  of  momentum  and 
angular  momentum,  as, 

j  —  In  /I 

m  4*£oWv3  ’  (2-6) 

where  A  is  the  ratio  of  Debye's  shielding  distance,  X^,  and  the  impact 
parameter,  b^.  For  a  Maxwellian  distribution, 

-A*  _  12  Tt  (g.feTP1  ,,  7s 

"to 

When  the  electron  temperature  exceeds  about  4x10^ °K,  A  must  be  reduced 
because  of  quantum  mechanical  effects  by  a  factor  of  (4.2x10  /T)2. 

The  mobility  can  be  found  as  in  the  case  of  a  slightly  ionized  gas, 
except  that  neutral  particles  are  replaced  by  ions,  with  a  different 


collision  frequency  which  is  now  largely  dependent  on  v,  as, 


I 


- 


i  ' 
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„  _  64-J2TC  to  ('kTTr 
r  Zle3ni  U!\ 


(2.8) 


To  obtain  an  accurate  expression  for  a  in  an  ionized  gas,  electron- 
electron  encounters  must  be  taken  into  account.  The  result  given  by 
Spitzer  and  Harm  (2.3)  may  be  expressed  in  the  form, 

^  ~  n  ril  _  2-63x|o2  YeT^7  /i  /  v 

<r  —  Tgneefi  — - “ — -  (nnhos/m)  ,  ^2 . 

2  ImA 

with  the  expression  of  Eq .  2.8,  where  values  of  7,-,  which  depends  on  the 

£j 

ionic  charge ,  Z  ,  are  given  as , 


z 

1 

2 

4 

16 

00 

0.583 

0.683 

0.785 

0.923 

1 

Thus  in  the  important  case  of  Z=1 ,  we  have, 


<r  =  |-53  x/0 


V 

T* 


LA 


( lYihos/m) . 


(2.10) 


The  results  obtained  above  are  based  on  the  assumption  that  the  electric 
field  is  sufficiently  small  so  that  the  potential  energy  gained  across 


one  mean  free  path  is  negligible  compared  to  kT . 


(2.1)  S.  Chapman  and  T.G.  Cowling:  "The  Mathematical  Theory  of  Non- 
Uniform  Gases"  Cambridge  Univ .  Press,  London  Secs.  18.12,  14.12 
and  10.53  (1960) 

(2.2)  J.S.  Townsend  and  V.A.  Bailey:  Phil.  Mag.  43_  593  (1922),  44_  1033  (1922) 

(2.3)  L.  Spitzer  &  R.  Harm:  Phys .  Rev.  89_  977  (1953) 

L.  Spitzer:  "Physics  of  Fully  Ionized  Gases" 


.* 

. 


:  '  '  .  . 

. 


9 


2.1.3  Moderately  Ionized  Gas 


For  the  electrical  conductivity  of  a  moderately  ionized  gas,  the 
approximation  of 


1 


cr 


-L+-L 

<%  <?d 


(2.11) 


can  be  used,  where  a  is  the  conductivity  of  a  slightly  ionized  gas 

c 

(Eq.  2.5)  where  the  close  encounters  between  the  electrons  and  the 

neutral  particles  dictate  the  electron  mobility,  and  is  the 

conductivity  of  a  completely  ionized  gas  (Eq.  2.9)  where  the  distant 

encounters  between  the  ions  and  electrons  dominate  the  situation.  By 

comparing  the  relative  magnitude  of  the  close  and  distant  encounter 

terms  for  argon,  the  two  terms  tr,  and  will  be  equal  at  about  7000°K, 

or  a  degree  of  ionization  of  about  10  ^ .  At  temperatures  higher  than 

-3 

about  8000°K  (o:>10  )  ,  the  close  encounter  term  will  become  so  insignificant 


that  the  resultant  conductivity  will  be  practically  the  same  as  . 
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FIG.  2.1  Effective  elastic 
collision  cross  section  of  argon 
atoms  for  scattering  electrons  as 
a  function  of  electron  temperature 
(2.2). 
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2.2  TRANSPORT  PHENOMENA  OF  PLASMA 


2.2.1  Self -Collision  Time 


The  self -collision  time,  ,  is  defined  by  1 /2i^m  when  the  test 


particle  is  the  same  kind  as  the  field  particles,  and  the  test  particle 


has  a  velocity  of  (3kT/m^)2,  as  (2.4), 

6/3  Tr£?/m  (kT)^ 


tc= 


II4x|06/F  T * 


nZV  U  (Slim -QUty  nZ4  In  A 


(sec) 


(2.12) 


Here  we  have  let  m  equal  to  Am^,  where  m^  is  the  mass  of  unit  atomic 

weight.  tc  provides  a  measure  both  of  the  time  required  to  reduce 

substantially  any  lack  of  isotropy  in  the  velocity  distribution  and  also 

of  the  time  required  for  the  distribution  of  kinetic  energies  to  approach 

the  Maxwellian  distribution.  The  value  of  t  determines  also  the  value 

c 

of  y  to  be  used  for  computing  T  in  an  adiabatic  compression.  If  the 
compression  is  slow  compared  to  t  ,  7  equals  5/3  (three  dimensional 
compression).  For  compression  more  rapid  than  t  ,  Twill  equal  2  or  3 
depending  on  whether  the  compression  is  perpendicular  to  the  lines  of 
magnetic  force  (two  dimensional  compression)  or  parallel  to  magnetic 
force  (one  dimensional  compression) . 


2.2.2  Equiparti tion  Time  for  Charged  Particles 


Let  us  suppose  that  the  test  particles  and  the  field  particles 


(2.4)  S.  Chandrasekhar:  "Principles  of  Stellar  Dynamics"  Univ.  of 
Chicago  Press  (1942) 
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both  have  Maxwellian  velocity  distributions,  but  with  different  kinetic 


temperatures  T  and  T^.  The  rate  at  which  equiparti tion  of  energy  is 
established  between  two  groups  of  particles  is  called  the  time  of 
equiparti tion  and  expressed  by  Spitzer  (2.5)  as, 

>»mf V*  jT  |  Tf 


\m  m<\ 


,6  AAf  (T  ,  Tf) 


=  5.87  x(0 


hfZ’ZfM\A  At/ 


(2.13) 


For  a  proton-electron  mixture, 


tc (proton) 

t  (electron) 
c 


=  43 


^£L 


t  (proton) 
c 


=43 


t  (electron) 
c 


=1836. 

(2.14) 


2.2.3  Equiparti tion  Time  for  Ions  and  Electrons  with  Neutral  Particles 

In  the  case  of  a  center  collision  between  an  ion  and  a  neutral 
particle  of  the  same  kind,  they  exchange  their  total  energies.  If  we 
average  the  energy  exchange  rate  with  respect  to  collision  parameter 
between  ions  and  neutrals,  it  becomes  approximately  half  of  the  value 
of  the  rate  in  the  case  of  center  collisions.  During  At,  an  ion 
experiences  v^At/\^  collisions  with  neutrals,  so  that  the  total  ion 
energy  per  unit  volume  at  t+At  is, 

-f-nj  kt]  (t  +4t)  =  l-nikTi  (t)  - 

and  for  neutral  particles , 

4»»kT,  <t+*t)  =  lnik^4t  • 


(2.15) 

(2.16) 


(2.5)  L.  Spitzer:  Monthly  Notices,  Roy.  Astron.  Soc.  (London)  100 
396  (1940) 
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By  solving  these  equations  with  the  boundary  conditions  that  T^(0)=T\q 
and  Tm(0)=Tm0,  we  obtain, 


T(t)=  i^Tio  +  ,  ~Tio  ~  T^o 

I  +«  i  +<*  ^ 

-  A~fio  __  ^(Tfo-T^o)  p'ixi  ((^)t 

I  +  <X  I  -t-oc  ^ 


(2.17) 


(2.18) 


where  a  is  the  degree  of  ionization.  The  time  of  equiparti tion  for  ions 
and  neutrals  of  the  same  kind  can  then  be  expressed  as , 


2\i 

Vi  (t+<*) 


Al 

l-HX 


I.^Sxlo2-, 


li 

l+<X  • 


(2.19) 


The  last  two  terms  of  the  above  equation  are  valid  only  when  ions  have 
a  Maxwellian  velocity  distribution. 

From  the  classical  theory  the  energy  loss  of  an  electron  due  to 
a  collision  with  a  neutral  particle  is  (2.6), 

Here,  it  is  assumed  that  the  kinetic  energy  of  the  particles  is  much 
smaller  than  the  internal  energy.  By  constructing  the  same  kind  of 
equations  as  Eqs .  2.15  and  2.16  for  electron -neutral  collisions  we  find 
that  the  time  of  equiparti tion  for  electrons  and  neutral  particles  is, 


3  nr>n  Ae 

8me"^  (I-Hx) 


=  o.ni 


O+^fTl  ' 


(2.21) 


2.2.4  Time  Dependent  Ambipolar  Diffusion 


For  the  ambipolar  diffusion  of  ions  and  electrons  in  a  plasma 
consisting  predominantly  of  neutral  particles,  the  following  equations 


(2.6)  A.M.  Cravath:  Phys .  Rev.  36_  248  (1930) 
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can  be  obtained  from  the  Boltzmann  equation, 

m  =  -hiVni+niHi(E 
(T?  =  -BeVne-neHe(E-^-|| 
and  ITT  =  [Je  SIT ,  ni=ne=n. 

By  eliminating^  ,  we  obtain. 


r  = -bav  n  - 


i  sF 


Va  St  7 


where 


Pjfe  +De(*i 


/V 


k(Te  +Ti) 


and 


^me  Mn  i  (He'+HO  ^  ^ 


(2  .22) 
(2  .23) 
(2.24) 


(2.25) 

(2  .26) 

(2.27) 


Here  Da  is  an  ambipolar  diffusion  constant  and  is  an  ambipolar  momentum 
transfer  collision  frequency. 

Plasma  diffusion  in  a  cylinder  of  radius  r^  and  length  L,  as  a 
function  of  position  and  time  is  now  considered.  From  the  law  of 
continui ty , 


V-F  + 


dh 

at 


=  o 


(2  .28) 


and  Eq .  2.25  we  obtain,  in  a  cylindrical  coordinate,  with  d/dd=  0 


\72n  =  J_£!L  4.  a*n 
v  f  ar  T 


_  .  j__  an_ 

31T2-  T  az*  Da  at  i4Da  at2-. 


(2  .29) 


By  putting  n=R(r )Z (z )T( t )  and  with  the  boundary  conditions  that  n=0  at 


z=+L/2  and/or  r=rQ ,  the  above  equation  can  be  solved  as, 


n=noj0(^r)  Q5S(f-2){(«+^)  -  (jiWa)Qp} ,  (2-30) 

where  Jq(x)  is  a  Bessel  function  of  Oth  order  of  the  first  kind,  2.405 
is  the  first  zero  of  Jq(x)  and 

<*l_  -v.  ±  U*  - 

/*J  2 - 


1  -  ■  J.  j  i  f-  err,  j  :  •  :  .  > /"  t  ‘  .  '  ■  ■  •  : 
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If 


(Z&kfcff, 

n-r.M^r)  c°s(j-2)e"h  , 


then  the  equation  is  reduced  to, 


(2 .32) 


where 


fa  = 


(2.33) 


is  the  life  time  of  the  plasma. 


2.2.5  Thermal  Diffusion 


If  nQ  particles  at  r=0 ,  t=0  are  in  a  temperature  gradient  field,  the 
particle  density  at  (r,t)  becomes,  in  the  spherical  coordinate  (2.7), 


n  (r,  t)  = - -  g  4Dt 


(2.34) 


Generally  the  life  time,  r  ,  can  be  approximated  by 


Tt  = 


Lz 


(2.35) 


2fD  , 

where  L  is  the  characteristic  length  and  f  is  the  degree  of  freedom.  The 
root  mean  square  displacement  of  a  particle  from  its  original  position  at 


t  is  , 


fr*  , 


(2.36) 


or,  85%  of  n^  is  confined  within 


(2.7)  H.  Raether:  "Electron  Avalanches  and  Breakdown  in  Gases" 
Butterworth  and  Co.,  Ltd.  (1964) 


. 

.  - 
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2.2.6  Electron -Ion  Volume  Recombination 


The  rate  of  change  of  electron  density  due  to  electron-ion  recombination 
can  be  expressed  as, 


Jme 

dt 


-ofen/ 


-i-  =  <xet  +■  4- 

He  rteo 


(2.37) 

(2.38) 


where  a  is  the  electron-ion  recombination 
e 

density  at  time  zero.  By  defining  the  life 
«e(  Tr)=neo/2  >  then, 

T  _  _! _ 

r  «ene0  ' 

In  the  case  of  argon,  by  using  the  value  of 
by  Kenty  (2.8), 


coefficient  and  ngQ  is  the 
time  of  electrons  ,  r  ,  as 


(2.39) 

16  3 

ae=2xl0  m  /sec  measured 


t  =0.5  ms 
r 

and  rr=0.5  ps 


1 9 

-3 

at 

ne0=1° 

m 

22 

-3 

at 

ne0=1° 

m 

(2.8)  C.  Kenty:  Phys .  Rev.  32  624  (1928) 

The  value  of  ae=3xl0-^  m^/sec  obtained  by  M.A.  Biondi  and  S.C. 
Brown  (Phys.  Rev.  76  1697  (1949))  seems  to  be  too  large. 


- 
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2.3  ION  CURRENT  IN  A  COLLISIONLESS  PLASMA 


The  appropriate  theory  to  use  for  calculating  ion  saturation  currents 
for  cylindrical  probes  in  a  collisionless  plasma  is  shown  in  Table  2.1 
for  various  ranges  of  parameters.  /3(=E^/ZkTe)  is  a  normalized  ion  energy 
and  Xp  (=-eVp/kTe)  is  a  normalized  probe  potential. 

The  probe  voltage  enters  here  because  it  affects  the  distribution  of 
the  electrons.  For  Xp^l  the  electron  distribution  may  be  considered 
Maxwellian.  For  Xp  less  than  about  1,  the  quasi -neutral  solution  holds 
everywhere,  and  probe  theory  is  particularly  simple.  For  l<Xp<5 ,  the 
deviation  from  a  Maxwellian  distribution  due  to  the  loss  of  electrons  to 
the  probe  must  be  taken  into  account.  Although  this  poses  no  problem  in 
principle,  it  complicates  the  computations. 

For  /3<$C1 ,  the  dependence  on  ion  energy  is  slight.  For  (8  >$>1 ,  the 
probe  current  depends  primarily  on  kT^  and  is  given  by  Langmuir's  sheath 
theory  for  a  Maxwellian  distribution.  The  transition  case  /3— 1  is  not  well 
covered  by  any  simple  theory.  The  simple  sheath  theory  fails  because  the 
accelerating  electric  field  outside  the  sheath  is  neglected,  and  the  BR 
theory  fails  because  a  monoenergetic  distribution  is  no  longer  a  good 
approximation . 

To  avoid  the  confusion  the  following  notations  will  be  used  throughout 
I:  current  (Amps) 

J:  current  per  unit  length  (A/m) 


current  density  (A/m  ) 


■ 


' 

■ 


I 

■ 


:T 


18 


Ion  Energy 
Distribution 

P 

rp<<XD 

rp- 

rP»xD 

all  X 

P 

Xp»  1 

V1 

Xp»l 

5>Xp>l 

Xp<i 

Zero  Energy 

0 

LAM 

LAM 

<&-  1 

LO 

BR 

LO 

LAM 

LAM 

LAM 

Mono 

^  1 

LO 

BR 

LO 

LAM 

LAM 

LAM 

Energetic 

^  1 

LO 

BR 

LO 

LS 

LO 

1 

LO 

BR 

LO 

BR 

Maxwellian 

^  1 

LO 

BR 

BR 

»  1 

LO 

BR 

LO 

LS 

LO 

Table  2.1  The  appropriate  theory  to  use  for  calculation  of  ion 
saturation  currents. 

(5  and  Xp  are  normalized  ion  energy  and  normalized  probe 
potential,  respectively.  'LS'  stands  for  Langmuir  sheath 
theory  (2.11)  (Sec.  2.3.1),  'LO'  stands  for  Langmuir  orbital 
motion  theory  (2.10)  (Sec.  2.3.2),  'LAM'  stands  for  Lam's  theory 

(2.9)  (Sec.  2.3.3)  and  'BR'  stands  for  Bernstein  and  Rabinowi tzs ' s 
theory  (2.12). 


(2.9)  S.H.  Lam:  Phys .  Fluids  8  73  (1965) 

(2.10)  H.M.  Motto-Smith  and  I.  Langmuir:  Phys.  Rev.  23.727  (1924) 

(2.11)  I.  Langmuir:  Phys.  Rev.  2_8  72  7  (192  6) 

(2.12)  I.B.  Bernstein  and  I.  Rabinowits:  Phys.  Fluids  2_  112  (1959) 


'  \ 

■* 

I « 

•do'  ci  si  Itiiiofl  b  e  yg'xcin-)  nc  i 
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2.3.1  Langmuir  Sheath  Theory 


The  ion  current  flowing  into  a  cylindrical  probe  per  unit  length 


is  expressed  as , 
J=27rrgjr , 


(2.41) 


where  r  is  the  sheath  radius  and  j  is  the  ion  random  current  density 

S  1C 

which  takes  the  form  shown  in  the  section  "List  of  Symbols".  The  radius 
of  the  sheath  can  be  found  from  the  Chi Id -Langmuir  Law, 


T-  aJ2Tc€ofi~  IVfb  /,  .  2.66  \ 

J  <?  1  ">i  fp/J1  v+  iff  I  ’ 


where  v  =-eV  /kT.  (not  T  )  and  R  is  the  solution  of 
p  p  i  e  K 


(2.42) 


(2.43) 


with  T=ln(r^/rg)  (<0) .  The  value  of  as  a  function  of  rg/r^  was  calculated 
by  Blodgett  (2.13)  and  shown  in  the  Table  III  of  her  paper.  The  1+2 .66/y/rj~^ 
term  arises  from  the  fact  that  the  ions  have  an  initial  energy  given  by  the 
Maxwellian  distribution.  When  the  sheath  thickness,  d,  is  thin  compared 
to  the  probe  radius,  Eq .  2.42  can  be  simplified  to, 


T-  812* &/T  r?  IVrP  /  2.66 \  3.42 WO7 rP  M*  (. .  2.66\ 

J-  i  J<r>i  d1  jx a*  ('  W) 


(2.44) 


and  The  potential  distribution  inside  the  sheath  for  the  case  is, 

(2.45) 

and  the  intensity  of  the  electric  field  at  the  probe  surface  can  be 


VOr)  =-Vp  (-£-)*, 


found  to  be  larger  by  a  factor  of  4/3  than  the  averaged  value  of  the 


(2.13)  K.B.  Blodgett:  Phys .  Rev.  22_  115  (1923) 


t" 

' 


!  • 
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intensity  in  the  sheath  (of  V^/d). 

The  ion  sheath  can  be  formed  within  a  time  corresponding  to  the  ion 

plasma  frequency,  a;  ,  given  in  the  section  "List  of  Symbols".  In  the 

10  -3. 

case  of  an  argon  plasma,  this  time  varies  from  0.4  /us  (ne=10  m  )  to 

2  2 

0.4  ns  (n  =10  m  ^) . 
e 


2.3.2  Langmuir  Orbital  Motion  Theory 


The  ion  current  flowing  into  a  cylindrical  probe  per  unit  length  is 
expressed  as , 


J=27rr  j  F 

pJ  r 


(2.46) 


where 


F=  -jr-erflU  +(l-ei-i\|xP+?)e*P, 

5=  (2.47) 

In  the  limiting  case  of  r  -r  r  ,  F=r  /r  and  we  recover  Eq .  2.41.  In 

&  s  P  p  ’  s  p  n 

the  practical  applications  involving  small  probes  and  positive  ion 

saturation,  we  have  r  »r  and  leV  /kT  ^>1.  With  these  conditions  Eq .  2.46 
*  s  '  p  f  P  e  n 


becomes , 


J=2icrP>-^-erfJl. 


(2.48) 


The  sheath  size  is  determined  by  Eq .  2.42.  By  equating  Eqs .  2.42  and  2.48 
we  obtain  the  ion  saturation  current  reaching  the  probe  as  (2.14), 


J = 2TcrpJv  •  o.7i05V\ 


(2.49) 


(2.14)  G.J. Schulz  and  S.C.  Brown:  Phys .  Rev.  98^  1642  (1955) 


*  ‘ 
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where 


(2.50) 


When  the  electron  temperature  is  much  higher  than  that  of  the  ions,  the 
drift  current  density,  j  ,  reaches  the  sheath  edge,  rather  than  the  random 
current  density,  j  .  The  drift  current  at  the  sheath  can  be  taken  into 
account  by  replacing  j  by  and  correspondingly,  Nr  by  .  With  a  fixed 

probe  and  discharge  conditions  we  obtain, 

^O.SS  .  -0.2  .  ,  0.55  O.fc 

J- 2rcrPjd •  o.7 £  Md  «=VP  ne  .  (2.51) 


In  the  thick  sheath  limit  and  )(?>!,  Eq .  2.47  becomes, 


(2.52) 


and  we  have , 


z 


ill 

J  =2j2  rpneea mj"2  Vp1-  ^  Vp  ne-  (2.53) 

2  2  2  3 

The  slope  of  J  vs  V  yields  8r  n  e  /m-  ,  hence  the  value  of  n  . 


2.3.3  Lam's  theory 


To  obtain  the  ion  saturation  current  when  r ^ An , 

P 

consider  the  potential  distribution  arQund  the  probe. 


equations  are: 


we  have  to 
The  given 


Boltzmann  eq . 

Poisson  eq . 

Energy  eq . 

Angular  momentum  eq . 


a>0-  T  hr)  av  u 


Ej.  =^-m(P+ 
L  =  nr>r26  . 


(2.54) 
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By  assuming  that  the  ion  energy  distribution  is  monoenergetic  (=E^), 
Poisson  equation  reads,  in  a  dimensionless  form, 


l  di  /*•  dX  \ 
Fdf(§df)  = 


7t  JtsJ 

^IF+x 


ZfjiTX 


for 

for  (2.54) 


is  the  normalized  absorption  radius  and  can  be  found  from, 


-jfe(P+x(*°  )_Jlf-)=0'  (255) 


12 

E^/eV  |^C (r^/\^)  ,  Poisson  equation 


is  , 


dn 


(Sdf)  -JX'^+^ex=o. 


(2.56) 


Here  dimensionless  parameters  are  chosen  as , 


k  ~  cAp 
X=-eV/fcTe 
/s  =  Ei/2KTi 
=  TCHAZTe 


Ji 

T  -  Ji  /  mi 
UN  &KXitfkTc]  2He2fo 


(normalized  radius) 

(normalized  potential) 

(normalized  ion  energy) 

(  )(  for  Maxwellian  distribution) 

(ion  current  per  unit  length) 

(normalized  ion  current) 


Jb  =  2  tg  n<?e  f >  /22kle/m{ 
ce 


(Bohm  current) 

(normalized  Bohm  current,  a  function  of  0 

only;  i  =1  when  T.=0  and  1.25  when  T. =ZT  ) 
B  1  ±  G 

(ion  current  increment  ratio  due  to  finite 


sheath  thickness). 


These  equations  were  solved  by  Chen  (2.15)  with  (5-0-0. 1. 


The  V  -J 
P 


(2.15)  F.F.  Chen:  J.  Nucl.  Energy  C7  47  (1965) 


■ 
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characteristics  and  the  shape  of  V(r)  for  the  cylindrical  probes  when 
p=0  are  shown  in  Figs.  10  and  6  of  his  paper,  respectively.  To  determine 
the  plasma  density  one  uses  the  following  formula, 

A(X)  =3.89*)0~6|\/pl  (~4  )*.  (2.57) 

Having  A(r)  from  the  experimental  data,  one  finds  x  from  Fig.  2.2  (page  10) 

Knowing  r,  we  find  n  from  J  . 

e  B 

When  r  and/or  n  is  so  large  that  £  ( J^r^»ZeV^/2m.  )  , 

p  e  P  ip  pi* 

T~1  and  the  ion  saturation  current  becomes  equal  to  the  Bohm  current 

(J.=J_)  which  is  independent  of  V  .  If,  in  addition,  T.=T  ,  then 
l  B  p  i  e 

Ji=JB=3.54rpnee|%p-  (=3.ft  X[6\ne]%  fer  At).  <2-58> 

The  ion  saturation  current  calculated  by  Schulz  and  Brown  (2.14)  is, 

Jd=2*rPn rPnee/^jj ,  (2.59) 

which  is  only  a  factor  of  1.08  larger  than  Eq .  2.58.  When  is  not  too 
large,  Eq .  2.59  can  still  be  used  as  an  order  of  magnitude  check  of  the 


plasma  density. 


•  •  •  •  I 


.r(i  h  ^ 


i 
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2.4  ION  CURRENT  IN  A  COLLISION-DOMINATED  PLASMA 


When  the  ionic  mean  free  path  becomes  shorter  than  the  probe  radius  , 
most  of  the  electric  field  around  the  probe  is  confined  within  the  distance 
of  an  ionic  mean  free  path  from  the  probe  surface.  Accordingly  the  current 
is  essentially  limited  by  the  ambipolar  diffusion  of  charges.  The  ion 
current  can  be  expressed  approximately  as, 


j=eDaWe  —  eba—  , 


(2.60) 


where  5  is  a  characteristic  length  such  as  a  probe  radius.  This  formula 
can  be  modified  using  D^kju^ (Tg+T^ ) /e  ,  M.  =e/m^.=eX7m/v  and  T  (average 
speed)  =  (8kT^/7rm^)  2  , 

j  =  -^neeV-Ji.  (2.6i) 

This  expression  shows  that  nJv/4  ions  enter  the  sheath,  suffer  many 
collisions  and  finally  a  proportion  A^/5  reach  the  probe  surface.  This 
A. J8  term  can  also  be  expressed,  according  to  Schulz  and  Brown  (2.14)  as, 


Ai 


3\\ 


S  Z»it;  2d  1  (2.62) 

where  v  is  the  collision  frequency  for  positive  ions,  and  t^  is  the  time 
spent  by  an  unscattered  ion  in  the  sheath  of  thickness  d. 

The  rigorous  continuum  theory  shows  that  the  ion  saturation  current 


density  is  (2.16), 

J  =  U7irhee 


kTi  Xi 


I**,, 

»S 


=2.B*I018^ 

T> 


(A/m2)  for  (2.63) 


(2.16)  V.M.  Zakharova,  Y.M.  Kagan,  K.S.  Mustafin  and  V.I.  Perel:  Sov. 
Phys .  Tech.  Phys .  5^  411  (1960) 

C. H.  Su  and  R.E.  Kiel:  J.  Appl.  Phys.  37_  4907  (1966) 

D. E.  Weissman,  W.E.  Scharfman  and  H.  Guthart:  Phys.  Fluids  10  464 
(1967) 


' 


......  . 
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where  L  is  the  probe  length  and  r  the  sheath  radius  expressed  as  (2.17) 

s 


and  is  equal  to  0.56  for  L/ r  =300  and  is  0.63  for  L/ r  =100.  The  calculated 

P  P 

results  of  r  /r  as  a  function  of  V  for  a  range  of  electron  densities  are 

s  p 

4 

shown  in  Fig.  2.3  when  T  =T  =10  K,  r  =0.05  mm  and  L=4  cm.  The  sheath 

i  e  p 

size  contains  a  logarithmic  term,  so  that  for  cylindrical  probes  the 
effect  of  variation  in  rg  to  the  probe  current  is  relatively  small. 


2.5  ION  CURRENT  IN  A  COLLISIONLESS  PLASMA  FLOW 


Electrostatic  probes,  when  used  in  a  plasma  flow,  show  different 

characteristics  from  those  in  a  stationary  plasma.  The  electron  temperature 

deduced  from  the  electron  retarding  region,  where  the  electron  current  is 

much  larger  than  the  ion  current,  should  be  correct  as  long  as  r  <\  , 

p  e 

because  of  the  relative  insensitivity  of  electron  collection  to  the 
magnitude  of  the  ion-neutral  and  neutral -neutral  mean  free  paths.  Also 
Tg  can  be  determined  uniquely  with  a  probe  of  any  geometry  and  directional 
alignment  with  respect  to  the  flow.  This  was  experimentally  verified  by 
Dunn  (2.18)  and  Sonin  (2.19). 


(2.17)  W.E.  Scharfman:  "Operation  of  Electrostatic  Flush  Probes  in  This 
Plasma  Layers"  Stanford  Research  Institute  Project  7712  (Oct.  1969) 

(2.18)  M.G.  Dunn  &  J.A.  Lordi:  AIAA  Journal  8_  1077  (1970) 

(2.19)  A. A.  Sonin:  AIAA  Journal  4  1588  (1966) 


. 


.  (  .  '  "  v  : )  f 
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FIG.  2.3  The  ratio  of  the  sheath  radius  to  the  probe  radius  as 
a  function  of  probe  voltage  and  electron  density  in  a  continuum  plasma. 
Tj_=Te=10^  °K ,  Tp=0 . 05  mm  and  L=40  mm. 


V  (V) 


100 


In  the  case  of  determining  the  ion  density  the  situation  is  more 


complicated.  In  the  region  of  the  free -molecular  (r  P<Vi  to  the 

transition  (r~A.  .)  region,  or  more  exactly,  as  long  as  r  <2.5A  , 

p  1-1  o’  °  p—  n-n 

the  ion  current  into  the  probe  (In )  aligned  with  the  flow  direction  has 
the  same  value  as  in  the  case  of  a  stationary  plasma.  With  the  larger 
probes  ion  collection  is  inhibited  by  the  relatively  small  values  of  A. 


1  -n 


and  A^  ^  (2.18),  while  with  the  probes  of  rp/A<3,  ion  current  is  larger 
than  expected  (2 .20) . 

For  probes  whose  axis  makes  an  angle  of  6  with  respect  to  the  flow, 
the  ion  current,  Ij_  ,  is  a  function  of  r  /\.  and  uxsin0/v  ,  where  u  and 

n  J-  s 


v  are  the  flow  velocity  and  the  ion  acoustic  speed  (=J( Y  kT  +y. kT.)/m.). 
When  the  ionic  mean  free  path  is  much  larger  than  the  probe  radius  ,  win 
has  been  determined  by  Kanal  (2.21)  to  be  given  by  the  expression,  for 
cylindrical  probes  , 


(2.65) 


—  fH  , 

p»,x)=J  e  t  at , 


Q-(n,x)  =|  eV  'at, 

J  & 


(2.66) 


If  we  assume  negligible  sheath  thickness,  this  equation  reduces  to 

IT  -  -Ik  -  2-  o ’*■  /  C_  \Y>.  3. 

In  an  approximate  theory  Ijl/I||  can  be  found  as  follows: 

I  =%n  ev  .27rreL 
^  e  s 

I  =n_eusin0.2r  L 
e  s 


(2.67) 


(2.20)  S.  Lederman,  M.H.  Bloom  and  H.F.  Widhope:  AIAA  Journal  6_2133  (1968) 
W.E.  Scharfman  &  W.C.  Taylor:  AIAA  Journal  8^  1067  (1970) 

(2.21)  M.  Kanal:  J.  Appl.  Phys .  35_  1697  (1964) 


0  * 

■ 


. 
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.  _L.  _  _2  u  shB  2.  ^ 

“  i «~^v^r>T=^c-  (2-68 

This  equation  approximates  Eq .  2.67  well  if  C  is  larger  than  2.  (Ij./I=l 
at  C=0.) 


2.6  ION  CURRENT  IN  A  COLLISION-DOMINATED  PLASMA  FLOW 


When  an  electrical  Reynolds  number  defined  by 


_  uL  __  uL 


Pa  H‘(kTe/e) 


(2.69) 


(L  is  the  length  typical  of  probe  dimensions  for  cylindrical  probe, 

u  is  the  undisturbed  plasma  velocity  relative  to  the  probe  and  D  is  the 
ambipolar  diffusion  constant)  exceeds  unity,  the  probe  current  characteristics 
begin  to  be  influenced  by  the  flow  effects  of  the  plasma.  In  addition  when 
the  flow  is  supersonic,  a  shock  will  be  formed  around  the  probe  (see  Fig. 5. 4) 
and  the  probe  current  becomes  dependent  upon  the  superposed  ionization 
and  recombination  processes.  In  such  a  case  the  chemical  and  dynamical 
effects  on  the  probe  current  can  be  considered  separately  if  the  shock 
effects  are  not  very  large.  In  this  section  we  are  mainly  concerned  with 
the  ion  collection  effect  under  the  assumptions  of  frozen  chemistry, 
frictional  heating  in  the  boundary  layer  and  recombination  effects. 

Stagnation  shock  effects  including  new  ionization  will  be  treated  in 


Sec .5.2. 


■ 
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2.6.1  Thin  Sheath 


A  mathematical  treatment  for  thin  sheath  case  was  done  by  Lam  (2.22) 


under  the  assumptions  that 


*=-r«',  x="|r>>,>  eXlE«kTe 


R^  x*<$c  i . 


and 


(2.70) 


(The  last  condition  implies  that  the  sheath  remains  thin  at  moderate 
values  of  probe  bias.)  It  was  found  that  the  electrical  perturbations 
to  the  plasma  can  be  divided  into  three  physically  distinct  and  mathematically 
uncoupled  regions,  namely  the  outer  region,  the  ambipolar  diffusion  region 
and  the  sheath  region.  The  probe  current  can  be  approximated  as  a  current 
diffusing  into  the  electric  boundary  layer  (of  thickness  <5 )  formed  around 


the  probe  ,  i . e .  , 


Ji =2icrPeDdW  *  2tcrPeDay-  =J2  icne Je^rpukTe, 


(2.71) 


where  5  is  expressed  as 


r  2rP  _  /irpPa 


_  |2rpHik~k 
eu 


(2.72) 


2.6.2  Moderate  Sheath 


2  2  2 

The  case  of  <*«1 ,  1 ,  R£oC  «  1  and  X  »1  was  treated  by  Clements 

and  Smy  (2.23).  The  diffusion  current  in  this  case  is  negligible  and  also 
the  sheath  thickness,  d,  becomes  dependent  on  the  angle  6  between  the  radii 


(2.22)  S.H.  Lam:  AIAA  Journal  2_  256  (1964) 

(2.23)  R.M.  Clements  &  P.R.  Smy:  J.  Appl.  Phys .  4^L  3745  (1970) 


' 


j  :  .  . 
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to  the  point  under  consideration  and  the  upstream  edge  of  the  cylinder 
(see  Fig.  5-4).  d  can  be  expressed  as, 


C*19)  =*  Scrp  (jcofe-cose^f 

\  4«eeu  /  Svn  6  . 

The  ion  current  to  the  probe  is, 

_L  . 

J;=2neeu(zal  Sin  8)^  =4  (3pi£0*>)+Cneeu)4  Vp 


(2.73) 


3  , 

,Ti/i 


(2.74) 


2.6.3  Thick  Sheath 


2 

When  *«  i ,  x»  1  and  re**»  1 ,  the  sheath  becomes  much  thicker  than 
the  probe  radius.  The  current  to  the  probe  is. 


=2  r  n  eu  , 


s  e 


(2.73) 


and  so  by  eliminating  r  with  the  space  charge  equation  the  ion  current 


becomes  (2.24), 

(2i4ige)*  (neeuVp)* 

(!» 


(2.76) 


2.6.4  Ion  Current  to  the  Parallel  Mounted  Probe 

When  the  probe  surface  is  parallel  to  the  flow  direction,  the 
expression  for  the  ion  current  should  be  somewhat  changed.  Let  us 
consider  the  thin  sheath  case.  The  problem  can  be  reduced  to  that  of 
finding  the  ion  current  into  a  flat  plate  probe  of  width  27rrp  and  length 
L  placed  in  a  plasma  flow  with  a  flow  direction  parallel  to  the  probe 
surface.  The  current  into  the  sheath  can  be  divided  into  two  components; 


(2.24)  R.M.  Clements  &  P.R.  Smy:  J.  Appl .  Phys .  40  4553  (1969) 


/ 
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ambipolar  diffusion  current  (1^)  and  convection  current  (Ic) .  For 
2  2 

r „orx  4*  1 ,  the  diffusion  current  dominates  the  situation  while  for 

hi 

2  2 

X  ^1  the  convection  component  dominates  the  current. 

The  diffusion  current  can  be  expressed  in  the  form  of  Eq .  2.61  as, 


neev  X 


la  =  2.urPL— — ^  =  Lnej2icrPuekTpi , 


(2.77) 


Here  the  quantity  6  is  taken  as  the  thickness  of  the  electrical 
layer,  because  the  supply  of  ions  is  large  enough  there; 

c  _  irp  _  1  _  Ik2v\ j)> 

JRe  J  eu  J 


boundary 


(2.78) 


The  sheath  thicknes 
a  continuum  plasma  , 

j  =  f-*bPi 


s,  d, 


i  .  e . 


d3. 


can  be  determined  from  Poisson’s  equation  for 
,  by  solving  dE/dy=  J/,£oHiE> 


(2.79) 


Therefore  the  convection  current  can  be  found  to  be, 

I  =27rr  dn  eu.  (2.80) 

The  rigorous  solution  (2.25)  indicates  that  the  convection  current  is  /2 

times  larger  than  the  above  approximate  result,  thus, 

±  3  i 

Ic  =  2rclp(tg0|<iL)*  (neeufVr* 

The  actual  current  to  the  probe  becomes , 

I=Id+Ic.  (2.26) 


(2.81) 


(2.82) 


(2.25)  R.M.  Clements  &  P.R.  Smy:  Private  Communication  (1970) 

(2.26)  By  measuring  the  current  into  the  planar  probe  which  was  aligned 
parallel  to  the  plasma  flow  as  a  function  of  the  applied  voltage 
and  also  by  measuring  n£  and  u  by  the  methods  described  in  Sec.  3.2, 
this  equation  was  tested  and  proved  to  be  correct  for  the  plasma 
density  range  of  lO^^-lO^-^  cm“3. 


.  ■ 
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2.7  FLOATING  DOUBLE  PROBE  THEORY  IN  A  GENERALIZED  CASE 


2.7.1  Potential  Difference  Across  the  Sheath 


When  an  electrode  is  immersed  in  a  plasma,  a  sheath  is  formed  around 

it.  When  the  potential  across  the  sheath  and  the  size  of  the  sheath  reach 

certain  critical  values,  an  electrical  breakdown  in  the  sheath  occurs  which 

is  the  subject  of  this  thesis.  Two  electrodes  of  approximately  the  same 

geometry  were  used  for  this  purpose.  These  were  connected  to  each  other  by 

a  variable  voltage  source,  which  was  completely  electrically  isolated  from 

any  other  systems,  e.g.,  shock  driver  capacitor  bank  system,  oscilloscope 

system.  For  this  reason  the  electrodes  could  be  regarded  as  floating  double 

probes  first  proposed  by  Johnson  and  Malter  (2.27)  (thus  we  refer  to  out 

electrodes  as  probes  hereafter) ,  except  that  in  our  case  the  distance 

between  the  two  probes  was  sometimes  very  large. 

Let  us  find  the  probe  current  and  the  potentials  of  the  two  probes 

with  respect  to  the  plasma.  Put  S,  V  ,  j ^  and  j ^  as  the  sheath  area, 

probe  potential,  ion  and  electron  saturation  current  densities,  and  T  , 

T^  and  Vf  as  the  ion  and  electron  temperatures  and  plasma  potential  outside 

the  sheath,  respectively.  Suffixes  1  and  2  refer  to  the  first  (negative) 

and  the  second  (positive)  probes.  V  =V,1-V£r.  is  the  spatial  potential 

c  1 1  t  z 

difference  of  the  plasma  at  both  sheaths  and  ^cj=^p2~^pi  (^O)  is 

applied  voltage  across  the  probes.  The  potential  difference  across  each 

probe  sheath  is  described  by  V=V  -V  .  For  simplicity  we  put  <£>  =e/kT  , 

t  p  0 

y=e/kT.  and  £  I .  =S.  j  .  _+S  J  .  _ .  Usually  V>|v  I  as  V  is  the  order  of 
l  l  l  il  2  i2  d  1  C'  c 

5  volts. 


(2.27)  E.O.  Johnson  and  L.  Malter:  Phys.  Rev.  J30  58  (1950) 


'  J  o'  *T9:  <  a  rf^)  (X  . 
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Four  different  cases  have  to  be  considered  depending  on  the  amount 
of  electron  saturation  current  to  the  positive  probe,  j e2 • 

(1)  S2  j  e2>  21^ -S^ )  ^^expC -(f>  (V^ -Vc)  )  :  The  potential  diagram  near 

the  probes  is  shown  in  Fig.  2.4-A.  We  have  the  following  basic  equations 


v,+v  =v  +V, 

1  c  2  d 

I  =S,  j  . ,  -S.  j  e 
p  I'll  lei 


~4>1V1 


“Ip  S2^i2"S2^e2e 


~^2V2 


(2.83) 

(2.84) 

(2.85) 


From  these  equations  we  obtain, 


4-, 


_  /  S2  Jez  \  4>2  ~Vc) 

SiOji  “Ip  +  Lp/ 


(2.86) 


This  gives  us  the  relation  between  I  and  V  .  When  V,  is  large  enough  to 

p  P  d 

satisfy  exp(01(Vd-Vc))»  ji2jei/jiije2  and  exp^  (Vd -Vc)  )»  si  j  el/ S2 J*  e2  » 

I  saturates  and  becomes  independent  of  Vd  explicitly  (I  is  dependent 


on  which  may  be  dependent  on  Vd) ,  i.e., 


Vhhi 
VW  -k 


In 


Sa  je2 


Ell 


V 


£1, 


(2.87) 

(2.88) 

(2  .89) 


(2)  s1J*il<  S2J*e2<'  2Ii"SljeieXp('"<^l^Vd"Vc^  “2Ii:  The  Potential  diagram 

in  this  case  is  whosn  in  Fig.  2.4-B  and  the  basic  equations  are, 


0 

0 

> 

> 

1 

r-H 

> 

11 

CM 

> 

(2.90) 

I  — S  i  -Si  G-<^1^1 
p  lJil  blJele 

(2.91) 

T  c  ~  ^9  Vo  O  . 

_Ip-S2 J i2e  2  2 "S2 J  e2 » 

(2.92) 

which  lead  to, 

Sijei  _4>,(Vd-V/c) 

Sl3i|“Ip  \  S2jj2  /  ^ 

(2 .93) 

This  gives  the  relation  between  I  and  Vd .  When  is  large  enough  to 


iJiorl qxs 


' 
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satisfy 


^<{>,(Vd-Vc)  ^  4-SiSi  ji2jei 


(S|0;i  -S2jeif'  * 


(2.94) 


I  reaches  the  saturation  region  and  then 
P 


Ip-Sijji 
V,  =Vd-Vc  --5T  id  -s-t 


^2  Ji; 


^2  SzOez  -S,dn 


Va— i-t 


SzJez  -Sijii 


(2.95) 

(2.96) 

(2.97) 


(3)  S  j  =  S  j  :  Equation  2.93  gives  the  relation  of  I  vs  F..  When 

l  el  1  ll  p  a 


V  is  sufficiently  large  as  shown  in  Eq .  2.94,  then, 
d 

W=s,jn 

4>, M-^)+ 


Sidei 


4>i+?*  S2ja  . 


(2.98) 

(2.99) 

(2.100) 


(4)  S  i  <  S.  ].  •  Equation  2.93  gives  the  relation  of  I  vs  V,.  When 
2  e2  lJil  P  d 


V  is  sufficiently  large  as  shown  in  Eq .  2.94,  then, 
d 

Ip  ~  jez 

=  3) -^n~S~  ■ 

•  i  S2je2  * 

In  addition,  if  ^2^* 62^^!^" il  ’  ^1  Can  as  > 


(2.101) 
(2 .102) 
(2.103) 


T*  Jil  * 


(2.104) 


The  variation  of  V  and  -V2  for  all  the  cases  is  shown  in  Fig.  2.5 

as  a  function  of  S  i  .  When  V,  becomes  large,  the  dif f erencebetween 

2  e2  d 


' 


■ 
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FIG.  2.4  Potential  diagram  near  floating  double  probes.  The 
dashed  and  solid  lines  indicate  the  plasma  potential  (V^)  and  probe 
potential  (V  ) . 


FIG.  2.4-A 


PROBE  1  PROBE  2 


Vf2 

Vp2 


FIG.  2.4-B 


FIG.  2.5  Variation  of  potential  difference  across  the  sheaths  of  floating  double 
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21.  and  S  j  -S  j  iexp(-0i(V  -V  )/2)  becomes  extremely  small  and  so  V. 
1  i  1  -L  e  1  a  c  i 

can  be  approximated  by. 


(2.106) 


(2.105) 


The  calculations  carried  out  above  apply  only  to  the  case  where  the 


collisional  effects  of  ions  and  electrons  with  neutral  particles  can  be 


neglected  in  the  sheath.  However,  even  though  the  ions  (only  ions)  are 
collision  dominated  in  the  sheath,  Eq .  2.106  can  be  proved  to  be  still 
valid . 

2.7.2  Analysis  of  Axially  Positioned  Double  Probe  Current 

One  of  the  interesting  applications  of  the  double  probe  is  to  put 

two  probes  at  different  axial  positions  in  the  shock  tube,  and  Z2  with 

distance  between  probes  d.  Two  kinds  of  plasmas  (precursor  and  shock 

heated  plasma  (6.5  ))  are  created  in  the  tube  and  so  ion  and  electron 

saturation  current  densities  at  each  probe  look  like  Fig.  2.6.  According 

to  the  theory  in  the  previous  section,  the  probe  current  is  equal  to  the 

ion  or  electron  current  at  the  negative  probe  whichever  is  smaller  (when 

V  is  sufficiently  large)  and  so  the  probe  current  can  be  observed  as 
d 

shown  by  the  thicker  lines  in  Fig.  2.6.  Line  A  refers  to  P^  negative,  and 
line  B  refers  to  P^  negative. 

This  type  of  double  probe  can  measure  a  very  low  plasma  density 

(such  as  precursor  plasma)  without  using  a  large,  plasma  disturbing  probe 
7  O 

(10  cm  J  can  be  measured  without  difficulty).  The  other  application  of 


this  probe  is  to  measure  plasma  velocity  very  accurately. 


f  .  BO  6  .  ■  c  >  :  }  j  !  ( 


-  )  f  1  .  '  i 
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FIG.  2.6  The  axially  positioned  probe  current.  Ion  and 
electron  saturation  currents  at  probes  1  (suffix  1)  and  2 
(suffix  2)  change  as  shown  by  the  thin  lines,  and  the  measured 
probe  current  changes  as  shown  by  the  thicker  lines.  Curve  A  is 
when  Probe  1  is  negative  and  curve  B  is  when  probe  2  is  negative. 


at  probes  1  and  2 


plasma  arrives 
at  probe  1 


at  probe  2 
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CHAPTER  III  EXPERIMENTAL  APPARATUS  AND  DIAGNOSTIC  TECHNIQUES 


3.1  EXPERIMENTAL  APPARATUS 


3.1.1  Shock  Driver  and  Shock  Tube 


The  shock  was  created  by  an  electro-magnetic  driver  which  propels  the 
ionized  gas  by  the  jxB  force  and  the  thermal  expansion  of  the  arc  heated 
gas.  The  driver  used  in  this  experiment  consisted  of  the  discharge  chamber 
shown  in  Fig.  3.1.  The  stainless  steel  center  electrode  was  a  button  type 
and  had  a  diameter  of  1"  and  a  curvature  of  3/4".  According  to  Fitch  and 
McCormick  (3.1)  the  volumetric  loss  of  electrodes  in  high  current  discharge 
switches  is  smallest  for  stainless  steel,  with  copper  the  second.  They 
also  found  that  copper  showed  least  profile  distortion  and  least  surface 
roughening.  For  these  reasons  together  with  the  fact  that  they  are  put 
in  a  high  vacuum,  the  electrode  materials  were  selected  as  in  the  figure. 

This  type  of  driver  was  studied  by  Pert  (3.2).  He  found  that  the 
pinched  gas  column  in  the  center  operates  as  a  conventional  metallic 
central  conductor  for  the  coaxial  gun  during  the  first  1.5  ;us,  after 
which  the  plasma  slug  was  not  accelerated  by  the  current,  since  the 
discharge  became  unstable  and  the  head  of  the  pinched  column  could  be 
seen  to  break  away  from  the  outgoing  plasma  slug. 

The  shock  tube  was  made  of  8  foot  long  (changeable),  2"  inner 
diameter  Pyrex  glass  tubes.  At  the  far  downstream  end  the  vacuum  pumps 
were  connected  through  an  expansion  chamber. 

(3.1)  R.A.  Fitch  &  N.R.  McCormick:  Proc .  I.E.E.  106A  117  (1959) 

(3.2)  G.J.  Pert:  Canadian  J.  Phys .  46_  2055  (1968) 
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3.1.2  Gas  Input  and  the  Vacuum  System 

The  gas  used  throughout  the  experiments  was  argon  and  was  obtained 
from  Canadian  Liquid  Air,  who  give  the  following  impurity  content  (%) : 


Ar 

N2 

°2 

H2 

h2o 

99.995 

0.0033 

0.0005 

0.0002 

0.0001 

The  vacuum  system  consisted  of  a  backing  pump,  Zeolite  sorption  trap, 
buffer,  mercury  vapor  diffusion  pump  (Edwards  High  Vacuum  Ltd.  Model  EM2), 
water  cooled  baffle  (CB02),  liquid  nitrogen  cooled  vapor  trap  (NTM2)  and 
butterfly  valve  (QSB2),  with  a  gross  pumping  speed  of  25  liters/sec. 

The  following  vacuum  gauges  were  employed;  a  hot  cathode  ionization 

-9  -3 

gauge  (10  -10  mmHg)  ,  Pirani  gauges  (one  0.0001-0.1  mmHg  and  two  0.01- 

10  mmHg  for  shock  tube  and  backing  system),  Vacustat  mercury  gauges  (0.001- 
1  mmHg  and  0.01-10  mmHg)  and  McLeod  gauge  for  the  calibration  of  all  these 
gauges . 

A  base  pressure  of  2x10  ^  mmHg  was  obtained  with  this  system.  Just 
before  taking  the  measurements  the  pump  was  isolated  and  argon  gas  was 
admitted  through  the  needle  valve  to  bring  the  system  to  the  desired 

pressure.  After  each  firing  the  entire  shock  tube  was  cleaned  to 

-6 

2x10  mmHg. 

3.1.3  Power  Supply 

Figure  3.2  shows  the  block  diagram  of  the  apparatus  used  in  this 
thesis.  The  circuit  parameters  of  the  assembled  bank  discharge  circuit 


. 


t-o\oIrr  8i  ’,u.  mOV 
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FIG.  3.2  Block  diagram  of  the  power  supply  and  diagnostic  system 
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were: 


Bank  capacitance 


240  yF 


Discharge  circuit  inductance 


0.13  yH 


Discharge  circuit  resistance 


10.4  mil 


Maximum  bank  operating  voltage  5  KV 


Maximum  bank  energy 


3  KJ 


Ringing  frequency 


25  KHz. 


3.2  DIAGNOSTIC  TECHNIQUES 

3.2.1  Pressure  Probe 

The  Piezo  electric  pressure  probe  was  flush  mounted  to  the  shock 
tube  wall  for  the  measurements  of  shock  heated  gas  pressure.  The  probe 
used  was  type  LD-25  (Atlantic  Research  Co.)  with  a  sensitivity  of 
0.15  V/psi  (=0.717  mmHg/V)  and  the  rise  time  of  1  /jls  .  The  sensitivity 
was  calibrated  by  using  another  shock  tube  (3.3). 


(3.3)  The  shock  tube  used  for  the  calibration  was  an  electromagnetically 


driven  shock  tube  with  diaphragm  separating  the  driver  and  the 
test  sections.  This  type  of  shock  tube  was  introduced  by 
P.R.  Smy:  Nature  193  969  (1962),  Rev.  Sci.  Instr.  3(3  1334  (1965) 
and  studied  extensively  by 

M.G.R.  Phillips:  Ph.D.  Thesis,  Univ.  of  British  Columbia,  Canada 
(1969)  with  the  conclusion  that  this  type  of  shock  tube  behaves 
in  a  manner  similar  to  cold  gas -driver  diaphragm  shock  tubes, 
in  which  the  Rankine-Hugoniot  relations  are  valid. 
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3.2.2  Plasma  Density  Measurements  by  Electrostatic  Probes 

For  the  measurements  of  plasma  density  the  floating  double  probe 
method  (2.27)  provides  the  most  powerful  tool  in  our  experimental  conditions 
and  plasma  density  region.  The  theory  for  this  probe  is  mentioned  in 
Sec.  2.7,  especially  in  Eq.  2.86.  Since  we  are  not  interested  in  measuring 
the  electron  temperature  (see  Sec.  3.2.4),  the  probe  was  used  only  under 
ion  saturation  conditions.  The  probe  current  in  this  region  has  been 
treated  in  Secs.  2. 3-2. 6. 

The  probes  used  were  mostly  hair  pin  type  tungsten  probes  which  can 
be  heated  and  cleaned  before  each  run  (see  Sec. 5. 3. 7).  These  probes 
had  a  diameter  of  0.1  mm  and  a  length  of  4  cm.  Many  other  types  of  probes 
of  the  same  kind  were  made,  e.g.,  thicker  probes,  shorter  probes  and 
probes  which  could  be  placed  parallel  with  the  plasma  flow  direction. 

For  the  measurements  of  electrical  breakdown  phenomena  the  same  type 
of  probes  were  used,  except  that  they  had  much  longer  tungsten  wires  and 
their  holders  had  long  sleeves.  By  doing  this  the  contaminated  part 
(about  1  mm  long  which  cannot  be  heated  enough)  of  the  tungsten  wire 
adjacent  to  the  probe  holder  was  well  removed  from  the  plasma  and  did 
not  affect  the  measurement  of  the  breakdown  voltage.  Figure  3.3  shows 
the  general  view  of  the  electrode  system. 

The  probe/electrode  power  supply  consisted  of  D.C.  batteries 
paralleled  with  fast  capacitors  mostly  of  capacitance  equal  to  4  P  F. 

(A  capacitor  of  up  to  240 p F  was  tested  and  it  was  found  that  as  long 
as  the  capacitance  was  large  enough,  there  was  no  change  in  the  measured 
probe/electrode  currents.  However,  when  the  capacitance  was  too  large, 
the  electrode  was  destroyed  at  each  breakdown  due  to  the  intense  heat  of 
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the  arc  spots.  )  The  output  impedance  of  the  power  supply  was  0.2  . 

The  probe/electrode  currents  were  measured  by  an  A.C.  current  probe  (Sec. 
3.2.8). 


3.2.3  Plasma  Density  Measurements  by  Microwave  Cavity 


The  resonant  cavity  method  of  measuring  plasma  density  is  a  very 

sensitive  technique  because  it  utilizes  the  high  value  of  Q  of  the  cavity, 

11  -3 

This  technique  is  useful  at  densities  up  to  about  10  cm  ,  before  the 
plasma  frequency  becomes  so  high  that  the  necessary  cavity  dimensions 
become  res trictively  small. 

When  a  low  loss  dielectric  of  permittivity  e  is  inserted  in  a 
microwave  cavity,  the  resonant  frequency  shifts.  From  the  perturbation 
theory  the  shift  of  the  resonant  frequency,  Af ^ ,  from  the  resonant 
frequency,  f q  ,  is  expressed  as  (3.4), 

feVc  ’ 


£] 
fo 


_l_ 

2 


(3.1) 


where  «q  is  the  permittivity  of  the  vacuum  and  Jd-c  is  the  volume  of  the 

cavity.  If  e  does  not  depend  on  position,  the  numerator  of  the  right 

/2 

E  dr  1  ,  where  Jd  r 1 

is  the  volume  of  the  dielectric  in  the  cavity.  The  permittivity  of  the 
plasma  is  expressed  as,  with  the  notations  that 

a 


=  and 


_  cpp  _  ne’ 

^  tu*  *“  6 o  m  <0*  i 


(3.2) 


(3.4)  S.J.  Buchsbaum  &  S.C.  Brown:  Phys .  Rev.  106  196  (1957) 

S.C.  Brown:  Proc .  2nd  Int'l  Conf.  Peaceful  Use  of  Atomic  Energy 
32  394  (1958) 
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since,  the  plasma  under  measurement  is  usually  of  low  density,  low 
collision  plasma  frequency  and  therefore  If  a  plasma  of  the  same 

shape  as  the  dielectric  is  replaced  in  the  cavity,  then  the  resonant 
frequency  shift,  4f,  is  expressed  as, 


_  i  leW 

T“  2  •  (3-4) 

2 

Here,  in  order  to  use  the  perturbation  theory,  must  be  much  smaller 

than  co2  . 

Having  considered  the  actual  plasma  geometry  (2"  diameter,  infinite 
length),  interfering  modes,  cavity  size,  etc.,  a  cylindrical  cavity  of 
radius  R=4M  and  length  4.25"  surrounding  the  shock  tube  (of  radius  a) 
was  built  (3.5).  Further,  the  TMq2q  mode  was  found  to  be  adequate  for 
the  measurements.  For  this  mode  the  electric  field  and  resonant  frequency 
are , 


Ez=AJ0(kz),  EX=E  =0  (3.5) 

clc 

and  f0=-^j~  =2.59  GHz  (S-band)  ,  (3.6) 

where  c  is  the  speed  of  light,  k=S02/R  and  Sq2=5.520  is  the  second  zero  of 
Jq(x).  Integrals  of  Eq .  3.4  then  become. 


_  a1  J0a(ka) +Ji\ka.)  „ 

337 ' 


(3.7) 


Since  an  actual  resonant  cavity  has  two  end  holes  for  the  shock  tube 
to  pass  through,  Eq .  3.7  is  not  exact.  To  minimize  the  consequent  error 
two  improvements  were  made.  One  of  them  was  to  attach  conducting  sleeves 
which  tightly  surrounded  the  plasma  container  on  the  outside  of  both  holes 
on  the  cavity.  These  conducting  sleeves  were  designed  to  be  circular 
waveguides  whose  cut  off  frequency  was  beyond  the  resonant  frequency  of 


(3.5)  E.L.  Ginzton:  "Microwave  Measurements"  McGraw-Hill  Co.  pp357  (1957) 
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the  cavity  so  that  no  radiation  could  escape  and  hence  the  cavity  appeared 
to  act  as  if  it  had  solid  end  plates  (3.6).  These  sleeves  should  be  a 
few  attenuation  distances  long.  The  2.375"  diameter  4"  long  sleeves 
were  attached  on  both  sides  of  the  cavity.  The  sleeves  had  a  minimum 
cut  off  frequency  of  2.93  GHz.  The  resonant  frequency  of  this  assembled 
cavity  was  2407.5  MHz. 

The  second  improvement  was  to  calibrate  the  cavity  by  using  a  2" 
diameter  and  sufficiently  long  dielectric  rod  instead  of  plasma  and  so 
eliminate  the  integrals  of  Eq.  3.4.  The  following  equations  were  used 
for  determining  the  plasma  density. 


or 


n-  -o.oaKi4.rt 


*-<if 


(3.8) 


(3.9) 


To  know  the  precise  value  of  c/cq,  the  same  dielectric  was  coaxially 
mounted  into  a  cavity  with  solid  end  plates  and  the  shift  of  resonant 
frequency  was  measured.  From  the  shift,  was  calculated  by  using 

Eqs .  3.1  and  3.7. 


3.2.4  Electron  Temperature  Measurements  by  Single  Probe  Method 


A  floating  double  probe  collects  only  a  very  small  amount  of 
electrons  of  the  highest  energies.  In  our  plasma  where  there  possibly 
are  small  amounts  of  high  energy  electrons  (e.g.,  run  away  electrons, 
photoionized  electrons)  in  addition  to  the  large  amount  of  Maxwellian 


(3.6)  R.M.  Clements:  Private  Communication  (1969) 
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electrons,  the  double  probe  technique  is  liable  to  produce  erroneous 
measurements  of  bulk  electron  temperature. 

The  original  single  probe  method  proposed  by  Langmuir  and  Mott  -Smith 
(3.7)  collects  electrons  of  the  whole  energy  range,  so  that  this  method 
measures  the  of  the  majority  group.  However,  the  single  probe  technique 
has  disadvantages,  some  of  which  are  now  discussed: 

(1)  Sometimes  it  is  difficult  to  have  a  proper  reference  electrode.  In 
our  shock  tube,  the  plasmoid  becomes  electrically  separated  from  the  driver 
electrodes  and  so  the  driver  electrodes  do  not  provide  a  reference  potential. 

(2)  If  one  uses  a  large  enough  probe  as  a  reference,  then  the  geometries 
of  the  reference  and  the  electron  collecting  probes  are  very  different  and 
one  has  to  be  cautious  about  whether  they  are  in  a  collision  free  or 
collision  dominated  category. 

(3)  If  the  larger  probe  is  collision  dominated  (as  is  usual)  ,  the  current 
to  the  probe  is  mostly  determined  by  ambipolar  diffusion;  i.e.,  the  probe 
potential  extends  over  a  distance  which  is  approximately  equal  to  the 
probe  size.  Frequently  the  smaller  probe  cannot  be  put  outside  this 
ambipolar  diffusion  region.  In  this  ambipolar  field  the  ion  and  electron 
densities  are  still  equal  and  electrons  have  a  Maxwellian  distribution, 
although  their  density  and  potential  are  somewhat  decreased. 

(4)  If  the  larger  probe  is  in  a  collision  dominated  region,  the  current 
to  the  probe  cannot  increase  infinitely  as  the  probe  area  is  increased. 

This  can  be  proved  easily.  From  Eq .  2.60, 

I  =  SeDa-|s1  (3.io) 

where  6  is  the  characteristic  size  and  is  approximately  proportional  to 
S,  the  probe  area.  Therefore  I  becomes  independent  of  S  and  has  a  finite 


(3.7)  I.  Langmuir  &  H.M.  Mott  -Smith:  Gen.  Elec.  Rev.  27_  449 ,  538,  616, 
762,  810  (1924) 
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maximum  value.  In  this  case  the  ion  collection  by  the  larger  probe  is 
not  enough  and  its  potential  decreases  because  of  the  surplus  electron 
accumulation . 

(5)  The  probe  draws  a  heavy  electron  current  from  the  plasma  which 
results  in  a  plasma  perturbation. 

(6)  To  avoid  the  large  electron  currents,  the  electron  collecting  probe 
must  be  very  small.  This  causes  the  probe  geometry  to  be  far  different 
from  the  ideal  planar,  cylindrical  or  spherical  probe  geometry. 

One  of  the  ways  to  avoid  some  of  these  disadvantages  (1,  3  and  4)  is 
to  use  a  third  probe  (dummy  probe)  in  addition  to  the  electron  collecting 
probe  and  reference  probe.  The  current  into  the  electron  collecting  probe 


(of  area  S^)  is  expressed  as, 


=s,jee  kT£M-VdMy)  « 


_  eVf 

=  S|jee  kle  e  *Te  , 


(3.11) 


where  is  the  floating  potential  expressed  by  the  equation  in  Sec.  "List 

of  Symbols",  is  the  potential  of  the  electron  collecting  probe  with 

respect  to  the  reference  probe  and  AV  is  the  potential  float  of  the 

reference  probe.  This  can  readily  be  measured  by  applying  a  voltage 

across  the  dummy  and  the  reference  probes.  The  voltage  at  which  the 

current  into  the  dummy  probe  is  zero  is  equal  to  AV .  The  slope  of  log  I 

vs  V,-AV  when  Vr>V.-4V>0  yields  the  electron  temperature, 
d  t  d 

In  our  temperature  measurement  experiment,  the  reference  probe  was 

2 

made  of  copper  and  had  an  area  of  230  cm  .  Two  types  of  electron  collecting 
probes  were  made  of  cylindrical  tungsten  wire  and  had  dimensions  of  either 
0.1  mm  diameter  and  5  mm  long  or  1  mm  diameter  and  5  mm  long.  The  dummy 
probe  was  1  mm  diameter  and  5  cm  long  copper  wire. 
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3.2.5  Two  Coil  r.f.  Determination  of  Plasma  Conductivity 

This  section  describes  a  new  method  of  measuring  the  electrical 
conductivity  of  a  plasma  by  utilizing  the  dependence  of  the  mutual 
inductance  of  two  coils  upon  the  conductivity  of  the  medium  lying  between 
them  (3.8).  In  addition  to  its  simplicity,  this  method  has  several 
advantages  over  previous  single  coil  measurements  (3.9).  These  are; 

(1)  high  speed  response  (  20  nsec) , 

(2)  it  gives  a  measure  of  the  average  conductivity  in  the  plasma  rather 
than  of  that  very  close  to  the  coil, 

(3)  it  can  be  used  over  a  continuous  and  wide  range  of  frequencies  with 
the  result  that  a  very  wide  range  of  conductivity  can  be  measured 
(to  10^  mhos/m) . 

Electrostatic  effects  can  be  eliminated  by  the  use  of  a  third  coil  which 
avoids  any  necessity  for  electrostatic  screening. 

Consider  a  cylindrical  plasma  (of  radius  a)  with  surrounding  field 
and  search  coils  (of  radius  b)  placed  concentric  with  the  plasma  axis  with 
a  separation  of  Zq.  The  conductivity  of  the  plasma  is  considered  to  be 
constant  throughout,  i.e., 

=  (constant)  when  (Xr^a 

and  <i  =0  when  r>a ,  (3.12) 

where  r  is  the  radial  variable  of  a  cylindrical  coordinate  system  whose 
origin  is  at  the  center  of  the  field  coil.  Let  co  be  the  angular  frequency 
of  the  magnetic  field.  By  solving  Maxwell  equations  the  ratio  of  the 
signals  induced  by  the  search  coil  with  and  without  plasma  (V  and  V  )  is 
found  to  be  (for  the  case  of  a=2.54  cm,  b=3.00  cm). 
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V_=  alSS-^-j^jo^-jX^a]  +  59.Sj,[(Q(-jAzJ4a3 

Vo  (°<.-jX2)\Jo[(a2-jX^a]  -  8.96 J,  [(■*"- jAa;>  a.} 


(3.13) 


2 

where  X  =U)jJO'i  and  are  the  Bessel  functions  of  Oth  and  1st  order  of 
the  first  kind,  respectively.  Here  the  magnetic  field  was  approximated 
to  change  exponentially  along  the  axis  of  the  cylinder  with  the  decay 
constant  of  CX.  (This  is  a  good  approximation  if  z^>b/5.)  The  variation  of 
|v/v  |  as  a  function  of  o~f  was  solved  numerically  and  the  results  shown  in 
Fig.  2  of  reference  3.8. 

Measurements  have  been  performed  upon  electrolytes  of  known  conductivity 

and  plasmas  (whose  conductivity  was  measured  by  Lin's  method  (3.10))  with 

4 

conductivities  from  10  to  10  mhos/m  at  frequencies  from  0.1  to  50  MHz ,  and 
in  every  case  agreement  with  theory  was  excellent. 


3.2.6  Conductivity  Measurements  by  Lin's  Method 


For  conductivity  measurements  of  a  moving  plasma,  the  method  proposed 
by  Lin  et.  al.  (3.10)  provides  a  useful  tool.  This  method  uses  two  coils 
whose  axes  are  aligned  with  that  of  the  shock  tube.  The  downstream  coil 
(field  coil)  provides  an  axially  symmetric  D.C.  magnetic  field,  while  the 
search  coil  slightly  ahead  of  the  field  coil  picks  up  the  electromagnetic 
disturbance  produced  by  the  passage  of  the  plasma  slug  through  the  magnetic 
field . 


(3.8)  S.  Mikoshiba  &  P.R.  Smy:  Rev.  Sci.  Instr.  _40  1187  (1969) 

(3.9)  R.A.  Olson  &  E.C.  Lary:  Rev.  Sci.  Instr.  _33  1350  (1962) 

E.J.  Stubbe:  Proc.  IEEE  56  1483  (1968) 

P.  Savic  &  G.T.  Boult:  J.  Sci.  Instr.  _39  258  (1962) 

H.E.  Koritz  &  J.C.  Keck:  Rev.  Sci.  Instr.  3_5  201  (1964) 

K.B.  Persson:  J.  Appl.  Phys .  _32  2631  (1961) 

M.  Ciampi  &  N.  Talini:  J.  Appl.  Phys.  _38  3771  (1967) 

H.  Tanaca  &  M.  Hagi:  Jap.  J.  Appl.  Phys.  6_  335,  338  (1964) 

K.V.  Donskoi  &  I. A.  Prokof'ev:  Sov.  Phys.  Tech.  Phys.  7_  805  (1963) 
D.D.  Hollister:  AIAA  Journal  _2  1568  (1964) 

(3.10)  S.C.  Lin,  E.L.  Resler  &  A.  Kantrowitz:  J.  Appl.  Phys.  26  95  (1955) 
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The  D.C.  magnetic  field  induced  by  the  field  coil  is  proportional 
to  the  current,  I,  through  the  coil.  When  a  plasma  enters  the  field  with 
a  velocity  u,  a  current  is  induced  in  the  plasma  by  its  diamagnetic 
property.  This  current  has  only  an  azimuthal  component,  ,  and  is 
proportional  tocrul,  where  the  electrical  conductivity  of  the  plasma  is 
considered  to  be  constant  throughout  the  radial  direction  (eq.  3.12).  The 
total  magnetic  flux,  $ ,  due  to  is, 

$  ^  °°  <rUl  *  (3.14) 

therefore  the  voltage  induced  in  the  search  coil,  V,  is  proportional  to 
d<£/dt,  or 

(3.15) 

The  proportionality  constant  can  be  eliminated  by  calibrating  the  system 
with  a  metal  rod  of  known  conductivity  and  velocity.  Finally  we  obtain  the 
equation  for  plasma  conductivity. 


—  Uclc  JVdt  ^ 

IA  I  \\Jc4t  C  ’ 


(3.16) 


where  subscript  c  denotes  the  values  obtained  with  the  calibrating  metal 
rod . 


In  our  system  the  calibration  was  carried  out  by  an  aluminum  rod  of 

2"  diameter  and  6"  long,  with  the  results  that  u  =5.76  m/s  ,  <r  =2.26x10^ 

c  c 

mhos/m,  I  =145  A  and  fv  dt=1.6xl0  5  V.sec.  The  search  coil  signal  was 
c  c 

integrated  by  RC-integrator  with  a  time  constant  of  120  |is .  Equation  3.16 
can  now  be  rewritten  as , 

<r*  =0.944xl09  V  /u,  (3.17) 

L»KU 

where  V  is  the  voltage  measured  by  the  oscilloscope.  The  plasma 
CRO 

velocity  was  measured  simultaneously  by  a  photodiode  technique. 
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3.2.7  Plasma  Luminosity  Measurements 

To  measure  the  radiant  energy  from  the  plasma,  a  photodiode  (Hewlett- 
Packard  Associates,  Type  4205)  was  placed  adjacent  to  the  Pyrex  shock  tube 
perpendicular  to  the  tube  axis.  The  sensitivity  of  the  photodiode  was 
0.75  electrons/photon  at  7700  A,  the  directional  sensitivity  was  confined 
to  within  30°  (this  means  that  the  axial  resolution  was  3  cm)  and  the 
wavelength  response  was  0. 4-1.1  microns. 

The  pictures  of  the  shock  heated  gas  were  taken  with  an  image 
converter  framing  camera  (TRW  Instruments,  Model  ID)  with  a  plug-in 
unit  (Model  26B) . 


3.2.8  Oscilloscope  Assemblies 

The  following  oscilloscope  assemblies  were  used.  The  numbers  in 
the  bracket  show  the  upper  3  db  down  frequency  in  Hz.  The  lower  frequency 
limit  was  D.C.,  unless  otherwise  listed. 

Oscilloscopes:  Tektronix  type  555  (33M)  ,  454  (150M)  ,  549  (30M) ,  321A  (23M) 
Plug-in  units:  1A1  (50M) ,  CA  (24M) ,  1A5  (50M),  D  (2M) ,  G  (20M) ,  0  (23M) 
Probes:  P6006  (33M)  ,  P6012  (33M)  ,  P6047  (150M),  P6015  (50M,  40KV) , 

P6019  (120Hz-60M,  A.C.  current  probe  with  passive  termination) 


Camera : 


C-12  with  Polaroid  films. 
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CHAPTER  IV  PLASMA  PRODUCED  BY  ELECTROMAGNETICALLY  DRIVEN  SHOCK  TUBE 


4.1  PROPERTIES  OF  THE  GAS  BEHIND  THE  SHOCK 


4.1.1  Incident  Shock 


The  relations  between  the  properties  of  the  gas  behind  the  shock  for 
a  one-dimensional  flow  of  gas  are  presented,  taking  into  account  the 
effects  of  ionization.  In  a  reference  frame  moving  with  the  shock  front, 
whose  velocity  is  U,  the  equation  of  conservation  of  mass  is, 


(4.1) 


where  the  flow  quantities  upstream  of  the  interaction  region  are  designated 
by  suffix  "1",  downstream  quantities  by  suffix  "2".  The  equation  for  the 
rate  of  momentum  flow  takes  the  form  of , 

R+p.U*«Pa  +  fr(U-«u£  <4-2> 

The  equation  for  constant  energy  flow  takes  the  form  of  adiabatic  flow, 

Ht*  =  H2+*2  (U-a2)  ?  (4.3) 

where  H  is  the  specific  enthalpy  of  unit  mass  of  the  gas,  expressed  as, 


(4.4) 

Here  we  have  neglected  the  excitation  of  electronic  states  of  the  ions 


and  the  neutral  particles.  The  first  term  of  Eq.4.4  represents  the 
translational  energy  of  the  gas  particles;  the  quantity  1+a  indicates  the 
increase  in  the  number  of  particles  per  unit  mass  by  ionization.  The 
second  term  indicates  the  energy  of  ionization  per  unit  mass.  Together 
these  two  terms  stand  for  the  internal  energy  so  that  the  addition  of  the 
third  term  results  in  an  expression  for  the  enthalpy.  With  the  equation 


of  state. 
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the  expression  for  the  specific  enthalpy  may  be  simplified  to, 

u  _  _X_  P  ,  oieVi 

n~  T-|  J*  T>  • 


(4.5) 


(4.6) 


It  is  customary  to  introduce  the  Mach  number  of  the  shock, 

N  =  — 

ni  q  ’ 

where  c  is  the  sound  speed  in  the  undisturbed  gas , 


=/ 


C»=/T 


h 

J°. 


(4.7) 


(4.8) 


From  Eqs.  4.1,  4.2,  4.3,  4.6,  4.7  and  4.8,  we  obtain. 


o=j2_  (T+l)  Ml1 

1  +rM,2  +  — (^~»)6V;  M,1  , 

(4.9) 

(4.10) 

■£_  i+ei-O-rM,1 

(4.11) 

and 

(4.12) 

Here,  for  convenience,  we  have  normalized  u^  by  c^  and  put  equal  to  M. 
M  must  be  distinguished  from  . 

By  assuming  thermal  equilibrium,  the  expression  for  the  degree  of 
ionization  was  derived  by  Saha  (4.1)  using  the  method  of  statistical 
mechanics  , 


o<2= 


_ _  /0  IcT  -l  I 

3e3i  (2wnVyi  (kT)%  + ' 


Pl> 


eVi 


where  g^ ,  g  and  g^  are  statistical  weights 
and  are,  respectively,  1,  2  and  6  for  argon 
p  expressed  in  mmHg , 


.  (4.13) 

7 

of  neutral ,  ion  and  electron 
Especially  for  argon,  with 


(4.1)  Saha:  A  Treatise  on  Modern  Physics  Allahabad  and  Calcutta  (1934) 


1  ■  •  .  1 


('  .  0 


t  ‘  •  +■  t  • 


l  £  -  «  <*  ■  t  <  £ 


I 


56 


<xz= 


340  pT 


1.82*10* 

T 

4  1  • 

(4.14) 

2*  P2  *  T2 

,  u^  and  a  can  be  determined. 

This 

calculation  was  done  with  T-^=300  K,  M^=0-50,  p^=0. 1-100  mmHg ,  and  the 


results  presented  in  Figs.  4.1~4.5.  For  M^>50  Ar  becomes  significant 


and  at  M^=87  approximately  907,  of  argon  is  Ar  .  When  becomes  much 


+ 


larger,  Ar  becomes  dominant,  while  Ar  and  Ar  decreases (4 .2 ) . 

Especially  when  is  so  small  that  the  ionization  is  negligible, 
Eqs .  4.9-4.12  can  be  simplified  as, 


(t+i) 

A- 

P,  r+ 1 

J2  _  (2YM|-(r-i))(cr-i)M2f2) 


T, 


(T+irMi 


Uz  _  2  (M  7-l) 

c,  CY+I)!^, 


(' 


m*+5 

_  smNi 
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(_  (SHr 

V"  i 


-i)(m43) 


16  M,z 


{  _  3(M  2-l) 

\  4M. 

M  O-Mz  [hr-i)lvlit+2  ( _  I  t+ 

2  J  v  iTTxr 


Cl2 

Cz 

0 


2-(M, -i) 


7((r~0/iM)(2rIi  Ncr-o) 

(T+OM2 


30-tg) 


(- 

(- 


for  Ar)  (4.15) 


for  Ar 


) 


■) 


-) 


(4.16) 


for  Ar)  (4.17) 


for  Ar)  (4.18) 


for  Ar)  (4.19) 


4 


') 


for  Ar)  (4.20) 


■). 


for  Ar)  (4.21) 


For  the  actual  gas  7  changes  its  value  when  T is  raised  to  certain 
degree  by  dissociation  of  the  molecules  and  by  quantum  effects  (Sec.  2.2.1). 


(4.2)  W.A .  Menard  &  T.E.  Horton:  NASA  JPL  Tech.  Rept.  32-1408  (1969) 


rl  ' ■  fTv  r  i •  -  :  <.  1  T  :  ■■  i  ■  •  L 


e  r  i  '  o  oT09  1  3j:  X  “iq  B  \8^H  3*  bl.u 


<  -  jr  -  f  -  .  -  ■  ’ .  .  :  • 


(<  .  >) 

‘ 

' 


- 


. 


57 


M  vs  M-^  for  incident  argon  shock. 


FIG.  4.1 
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FIG.  4.2 


for  incident  argon  shock. 
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cnLt 


Q.  Q- 


FIG.  4.3 


P2/P1  vs  M-^  for  incident  shock. 
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T2/T^  vs 


for  incident  argon  shock. 


FIG.  4.4 


61 


FIG.  4.5  0(  vs  for  incident  argon  shock. 
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Roughly  speaking,  dissociation  becomes  important  at  0.3  eV,  dissociation 
is  about  complete  and  ionization  sets  in  at  1  eV  and  at  2.5  eV,  ionization 
is  of  the  order  of  50%.  The  general  expression  for  specific  heat  at 
constant  volume  can  be  written  as  (4.3), 


(4.22) 


4.1.2  Reflected  Shock 


When  a  primary  or  incident  shock  wave  reaches  the  end  of  a  shock  tube 


it  will  be  reflected  back  into  the  gas  which  has  already  been  heated  by 


the  passage  of  the  incident  shock  (4.4).  This  results  in  a  further 

increase  in  the  temperature,  pressure  and  density  of  the  gas.  To  derive 

the  equations  for  the  quantities  in  the  reflected  shock  it  is  more 

convenient  to  do  so  in  terms  of  the  conditions  in  the  incident  shock. 

From  the  equations  derived,  these  properties  may  then  be  related  to  the 

initial  conditions  by  way  of  the  Mach  number. 

With  a  reference  frame  moving  with  the  reflected  shock  front  whose 

speed  is  U  ,  the  Rankine-Hugoniot  relations  can  be  expressed  in  the  same 
R 

formas  Eqs .  4. 1-4. 6,  except  that  the  subscripts  1  and  2  of  p,  p,  T  and 
H  are  changed  into  2  and  3,  respectively,  and  that  U  and  changed  into 


(4.23) 


The  parameters  behind  the  reflected  shock  wave  can  be  seen  in  Eqs.  4.9- 
4.12  and  4.15-4.21,  with  the  same  changes  in  the  subscripts.  Equations 


(4.3)  J.K.  Wright:  "Shock  Tubes"  John  Willy  &  Sons  Inc.  N.Y.  pp24  (1961) 

(4.4)  A.G.  Gaydon  &  I.R.  Hurle:  "The  Shock  Tube  in  High  Temperature 
Chemical  Physics"  Reinhold  Publishing  Co.  N.Y.  (1963) 
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4.15-4.21  are  applicable  when  there  is  no  further  ionization  in  the 
reflected  shock.  From  Eqs .  4.18  and  its  equivalent  form  in  the  reflected 
shock  region,  we  obtain, 


uz=  2.(n?-i)c  „  2(nM. 

(T+  0  Ml  (T+ 1)  Mr  1 


(4.24) 


or 


,/2rMi  -(r-t) 

0M,z+2  • 


(4.25) 


From  this  equation  it  can  be  seen  that  even  with  M^=oo,  M^=(2Y/t-1)2 
for  argon).  Therefore  the  further  increase  in  the  density  and 
temperature  in  the  reflected  shock  is  not  very  large  and  Eqs.  4.15-4.21 
can  be  used  as  a  good  approximation  for  the  properties  of  the  reflected 
shock.  Also  if  M^>7  in  argon  shock,  then  we  can  approximate  that, 
MR=y5,  P^l  f>2=2  .5 ,  p3/p2=6,  T3/T2=2.4,  u2/c2=3 

and  UR/U=0 .5 .  (4.26) 


4.1.3  Generation  of  Shock  Heated  Gas  by  a  Shock  Tube 


For  some  shock  tube  conditions  a  contact  surface  is  never  established 
and  the  ambient  and  driver  gases  are  thoroughly  mixed.  Several  explanations 
of  this  behavior  have  been  put  forward;  Rayleigh-Taylor  instability  of  the 
contact  surface  (4.5),  inertial  instabilities  (4.6),  and  turbulent  mixing 
across  the  contact  discontinuity  (4.7).  According  to  Pert  (4.8)  the  flow 
is  dominated  by  turbulence  and  the  shock  heated  gas  is  thoroughly  mixed 


(4.5)  J.P.  Barach  &  R.E.  Vermillion:  Phys .  Fluids  8^  1976  (1965) 

G.I.  Taylor:  Proc .  Roy.  Soc .  A201  192  (1950) 

(4.6)  G.D.  Cormack:  Z.  Naturforsch  19a  934  (1964) 

(4.7)  J.P.  Barach  &  T.W.  Mayers;  Phys.  Fluids  11^  89  (1968) 

(4.8)  G.J.  Pert:  J.  Appl .  Phys.  39_  5932  (1968),  J.  Phys.  D,  Appl.  Phys. 
3  203  (1970) 
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behind  the  shock  front  at  low  initial  gas  pressure  (<0.5  mm)  and  high 
shock  speed  (>4  cm/ps),  so  that  no  separated  shock  layer  is  found.  The 
gas  behind  the  front  is  predominantly  due  to  the  collected  driven  gas  (for 
the  axial  distance  much  longer  than  the  length  of  the  driver)  and  accordingly 
the  state  of  the  gas  is  approximately  given  by  the  shock  relations,  if  the 
turbulent  energy  is  not  too  large.  However,  at  high  initial  gas  pressure 
(>1  mm)  and  low  flow  speed  «2 cm/ps)  a  shock  structure  is  formed  with  a 
well  developed  boundary  layer. 

The  conditionfor  onset  of  turbulence  in  the  flow  can  be  described 
in  terms  of  a  hydrodynamic  Reynold's  number, 


(4.27) 


where  V  is  the  viscosity  (values  for  argon  are  given  in  Fig.  4.6)  and  d  is 
the  diameter  of  the  expansion  tube;  i.e.,  with  Re  smaller  than  9000,  the 
shock  heated  gas  is  separated  from  the  driver  gas  (4.8). 


4.1.4  Flow  Duration  in  the  Shock  Tube 


The  length  L  of  the  slug  of  the  shock  processed  gas  contained  between 
the  shock  front  and  the  contact  surface  satisfies  the  relation  P2L=P-^z, 
where  z  is  the  axial  distance  from  the  driver  to  the  specified  point,  so 
that  the  ideal  flow  duration, X,  is  given  by 

•  (4-28) 

It  is  well  known,  however,  that  the  actual  flow  durations  are 
shorter  than  the  ideal  values;  a  factor  of  \  has  been  used  as  a  rule  of 
thumb  which  is  useful  at  higher  pressures  (>5  mmHg) .  At  low  values  of  p^ , 
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FIG.  4.6  Frictional  viscosity  for  argon  as  a  function 
of  temperature  at  1  atm  (4.9). 
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the  factor  can  be  much  less  than  The  severe  effect  at  low  pressures 
was  first  noted  by  Duff  (4.10),  who  found  that,  at  given  p^  ,  there  is  a 
limiting  value  oft,  i.e.,  t  cannot  be  increased  infinitely  with  increasing 
z.  He  gave  a  qualitative  explanation  of  the  loss,  noting  that  the  low 
speed  gas  in  the  boundary  layer  near  the  wall  "leaks"  past  the  contact 
surface.  The  theoretical  values  of  X  as  well  as  the  experimental  values 
are  shown  in  the  report  of  Roshko  (4.11).  These  measured  values  fall  below 
the  theoretical  real-gas  curve.  The  discrepancy  increases  with  increasing 
Mach  number  and  decreasing  initial  pressure. 


4.2  THEORETICAL  DERIVATION  OF  THE  PLASMA  VELOCITY 


The  model  used  for  the  electromagnetically  driven  shock  in  this 
theory  is  the  "finite  energy  snow  plough  with  leakage"  model.  Let  us 
consider  the  energy  balance  equation  so  that  ionization  effects  and 
frictional  energy  losses  to  the  shock  tube  wall  can  easily  be  taken  into 
account.  The  shock  driver  emits  a  driving  gas  slug  (plasmoid)  of  the  mass 
defined  by  mn^L  per  unit  area  of  the  shock  tube  cross  section.  L  is 
determined  experimentally  from  the  product  of  the  plasma  length  at  the 
driver  mouth  and  the  ratio  of  the  particle  densities  of  the  shock  heated 
gas  and  the  unheated  gas  (=n^/n^)  at  the  same  position.  (It  was  experimentally 
found  that  L  is  very  weakly  dependent  on  the  driver  energy  and  ambient  gas 
pressure,  which  makes  this  theory  simple  and  useful.  See  Sec.  4.3.7)  This 

(4.9)  R.S.  Devoto:  Phys .  Fluids  H)  354  (1967) 

(4.10)  R.E.  Duff:  Phys.  Fluids  2  207  (1959) 

(4.11)  A.  Roshko:  Phys.  Fluids  3  835  (1960) 
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plasmoid  works  as  a  piston  and  compresses  the  undisturbed  gas  ahead. 

By  doing  so  the  driving  gas  feeds  its  internal  and  kinetic  energies  into 


the  compressed  gas.  The  shock  heated  gas  is  not  well  separated  from  the 
driving  gas  by  the  contact  surface,  but  tends  to  diffuse  into  it  as 
mentioned  in  Sec.  4.1.4. 

In  the  case  of  the  pressure  driven  shock  tube,  the  actual  length  of 
the  shock  heated  gas  is  reduced  by  a  factor  k  (=0.1 -0.5,  depending  on  the 
initial  pressure  and  Mach  number)  from  the  theoretical  value.  For  our 
electromagnetically  driven  shock  tube  where  the  mixing  is  more  severe 
this  factor  k  is  smaller  and  can  be  regarded  as  a  ratio  of  the  number  of 
particles  which  moves  with  the  driving  gas  slug  to  the  total  number  of 
particles  which  are  affected  by  the  shock  (snow  plough  rate);  i.e.,  only 
kn^z  particles  out  of  n-^z  particles  start  moving  with  the  driving  gas 
slug,  and  the  others  lose  their  energies  by  frictional  force  between  the 
shock  tube  wall  and  by  diffusion  to  the  wall.  The  mass  of  the  driving  gas 
per  unit  area  at  z  can,  therefore,  be  expressed  as  mn^(L+kz). 

By  equating  the  energy  loss  of  driving  plus  shock  heated  gas  and  the 
energy  gain  of  undisturbed  gas  during  the  period  the  shock  proceeds  a 
distance  <dz,  we  obtain, 


(4.29) 


where  H  is  the  internal  energy  expressed  by  Eq .  4.6  with  the  assumption 


of  This  equation  can  be  rearranged  to. 


(4.30) 


By  eliminating  u  and  by  using  Eqs .  4.6,  4.9,  4.10  and  4.12,  we  obtain, 


L±lf2  JM)  -dTh  + 1  _0 


l-fk  dz 


•> 


(4.31) 
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where  JH)  =  fr*  O('l'l)  M !  (M,^  +  f-tf  +  2'rM.'l«-l)  +(W  M,*~ 

(7-i)i3eVi  da 


+ 


C*  m i 


dMl 


and 


fr-i)l*teVi 


m: 


(4.32) 


(4.33) 


By  normalizing  z  as  Z=(l-t-k)z/L,  equation  4.31  is  reduced  to, 

J(M.) 


dM,  n  +  l+k 


-z) 


S(Mi) 


=o. 


(4.34) 


The  first  term  of  J(M^)  arises  from  the  change  in  kinetic  energy,  the 
second  and  third  terms  from  the  change  in  internal  energy,  and  the  fourth 
term  from  the  energy  increase  by  ionization.  In  G(M^)  the  first  term  is 
kinetic  energy  and  the  other  terms  are  internal  energies. 

Equation  4.34  was  solved  numerically  with  the  boundary  conditions 
M=Mio  and  U2=u20  at  t=® »  with  k=0.1-l.  Figure  4.7  shows  how  each  term  in 
Eqs .  4.32  and  4.33  affects  the  plasma  velocity.  L=0.4  m  and  k=0.1  were 
used  here. 


Curve 

A: 

dtf/dM^O 

follows  the  Saha  equation 

Curve 

B: 

d(X/dM  =0 

de-ionization  energy  goes 
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C: 

dH  / dM, =0 
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internal  energy  change  is 
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Curve 

D: 

V° 

internal  energy  is  neglected. 

The  initial  pressure  dependence  was  found  to  be  very  small  for  p^=0.1- 
100  mmHg.  Especially  for  the  case  of  Curve  D,  Eq.  4.34  can  be  solved  as, 

Uo  Uo 


U-j  = 


/2+f  yi^h.2 + 1  * 


(4.35) 


In  deriving  Eq .  4.31  it  was  assumed  that  the  pressure,  density  and 
temperature  of  the  driving  gas  slug  are  the  same  as  those  of  the  shock 
heated  gas  whose  parameters  are  described  by  the  shock  relations.  The 
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validity  of  this  assumption  was  verified  experimentally  by  Pert  (4.8). 
However,  when  the  plasma  density  is  low,  the  loss  of  the  charged  particles 
by  ambipolar  diffusion  to  the  wall  becomes  larger  than  the  loss  by  volume 
recombination,  in  which  case  most  of  the  energy  released  by  de-ionization 
is  fed  into  the  wall,  not  into  the  oncoming  neutral  particles.  Therefore 
in  the  expression  of  we  may  set,  for  the  low  density  plasma, 

do( 

dMi  ■  (4-36) 

The  crossover  takes  place  when  the  volume  recombination  time  constant 

(T\.)  of  Eq.  2.39  is  equal  to  the  ambipolar  diffusion  time  constant  (Xa)  of 

Eq .  2.33.  For  1"  radius  shock  tube  with  the  pressure  of  1  mmHg,  the 

crossover  density  is  n  =2x10  cm  ^ .  If  the  neutral  density  becomes  higher, 

e 

then  Xa  becomes  longer,  recombination  dominates,  the  U2^z)  curve 
approaches  the  curve  A  of  Fig.  4.7  (broken  lines  of  Fig.  4.9),  and  so  the 
plasma  velocity  tends  to  have  a  lower  decay  rate.  On  the  other  hand  if 
the  plasma  density  becomes  lower,  then  *£  becomes  longer,  ambipolar 
diffusion  becomes  dominating  and  ^(z)  curve  approaches  the  curve  B  of 
Fig.  4.7  (solid  lines  of  Fig.  4.9).  It  is  worth  noting  that  the  curve  A 
can  be  approximated  well  by  the  curve  D  of  Eq .  4.35. 

In  Fig.  4.8  is  shown  the  variation  of  u^(z)  with  respect  to  k  when 
do</dM^=0.  Curves  for  u^(z)  with  various  initial  conditions,  M-^q,  are 
shown  in  Fig.  4.9  when  dCX/dM^  is  zero  and  non-zero. 

In  this  theory  the  expansion  wave  emitted  by  the  shock  as  it  slows 
down  was  neglected.  However,  this  effect  becomes  important  when  the  shock 
Mach  number  is  low  and  the  difference  between  U  and  u^  are  appreciable  that 
the  shock  deceleration  is  strong.  Consequently,  the  theory  introduced  here 
might  have  an  appreciable  error  at  low  shock  speeds. 
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FIG.  4.8  The  effect  of  k  on  the  plasma  velocity  decay. 
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FIG.  4.9  The  variation  of  the  plasma  velocity  as  a  function 
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4.3  MEASUREMENTS  IN  THE  SHOCK  TUBE 

4.3.1  Plasma  Velocity  Measurements 

The  velocity  of  the  plasma  was  measured  by  both  photodiode  and  r.f. 
conductivity  probes  with  the  time  of  flight  method.  The  photodiode  measures 
the  time  of  arrival  of  the  peak  luminosity,  i.e.,  the  peak  plasma  density 
(4.12),  and  the  conductivity  probe  measures  the  time  of  arrival  of  the 
peak  plasma  temperature  (Eq.  2.9).  The  experimentally  obtained  plasma 
velocities  by  these  two  methods  agree  well  as  seen  in  Fig.  4.10. 

Some  of  the  experimental  results  are  shown  in  Fig.  4.11  where 
theoretical  curves  are  also  shown.  The  agreement  between  them  is  very 
good  except  for  low  plasma  velocities  (M<5). 


4.3.2  Flow  Duration  Measurements 


Two  methods  of  measuring  the  flow  duration  were  employed;  the  two 
coil  r.f.  conductivity  probes  (Cases  1,3,4  and  5  in  Table  4.1)  and  the 
combination  of  piezo  electric  pressure  probe  and  electrostatic  probe 
(Cases  2  and  6  in  Table  4.1).  The  typical  oscilloscope  traces  for  the 
first  method  are  shown  in  Fig.  4.12.  In  the  second  method  the  pressure 


(4.12)  The  power  density  radiated  by  electrons  which  are  accelerated  by 
the  interaction  with  ions  (Bremss trahlung)  is, 

W=aiSTcVc"^l,  j^l  =  U4  x|6’40Z2ninejTe  O/M,  (4.37) 

which  is  approximately  proportional  to  the  square  of  the  electron 
density . 
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FIG.  4.10 


FIG.  4.12 
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FIG.  4.10  The  plasma  velocity  measurements  by  two  methods.  The 
upper  beam  is  the  signal  of  r.f.  conductivity  probe  (5  mV/vertical  div.) 
and  the  lower  beam  is  the  photodiode  signal  (1  V/vertical  div.)  Time 
base  is  50  ps/div.  Both  signals  were  taken  at  z=100  cm  with  p-^=0.15 
mmHg  in  argon  and  Vp=2.5  KV. 


FIG.  4.12  Typical  oscilloscope  traces  of  r.f.  conductivity  probe 
(case  3  of  Table  4.1).  Vertical  and  horizontal  deflections  are  0.05  V/div. 
and  200  ps/div.  Time  runs  from  left  to  right.  The  upper  trace  was  taken 
at  z=154  cm  from  the  driver  and  the  lower  trace  at  z=146  cm.  The  pinch 
of  the  upper  trace  during  0-100  ps  is  due  to  the  precursor  ionization, 
the  pinch  during  600-800  ps  to  the  shock  heated  gas  and  the  pinch  during 
800-1600  ps  to  the  driving  plasmoid. 
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probe  detects  the  arrival  of  the  shock  front,  while  the  electrostatic 
probe  detects  the  arrival  of  the  driver  gas. 

Some  of  the  results  are  summarized  in  Table  4.1  where  the  measured 
values  are  compared  with  theory.  Considering  the  results  of  Duff  (4.10) 
and  Roshko  (4.11)  the  agreement  can  be  said  to  be  reasonable. 


4.3.3  Pressure  Measurements 


Measurements  of  the  pressure  of  the  shock  heated  gas  give  us 
information  on  the  validity  of  the  shock  relations.  The  experimental 
values  of  p^  obtained  by  the  piezo  electric  pressure  probe  agree  well 
with  that  of  theory  (P2*)  when  Re<9000,  as  seen  in  Table  4.1.  When  Re 
becomes  larger,  the  indicated  values  are  greater  than  the  experimental 
values.  This  is  probably  due  to  the  leakage  of  the  shock  heated  gas 
through  contact  surface  along  the  boundary  layer  on  the  shock  tube  wall. 


4.3.4  Plasma  Density  Measurements 

In  the  first  measurements  with  the  floating  double  probes,  one  of 
the  probes  was  aligned  perpendicular  to  the  plasma  flow,  while  the  other 
was  parallel  to  the  flow.  Some  of  the  results  are  listed  in  Table  4.2  and 
the  typical  V-I  characteristics  for  these  measurements  are  shown  in  Fig. 
4.13  (which  corresponds  to  Cases  1,  2  and  3  of  Table  4.2).  With  a  probe 
voltage,  ,  larger  than  10  volts,  the  probe  currents  can  be  seen  to  be 
proportional  to  V  2,  as  expected  by  the  collision  dominated  theories  (Eq. 


2.74). 
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parallel/perpendicular  floating 
parallel  to  the  flow  while  the 

other  was  perpendicular  to  the  flow. 


FIG.  4.13  Typical 
double  probe.  One  of 


V-I  characteristics  of 
the  probes  was  aligned 
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The  ratio  of  lx / 1 1|  was  calculated  in  three  ways.  Since  the  parallel 
probe  also  was  a  hair  pin  type,  there  still  remained  a  small  portion 
(length  was  perpendicular  to  the  flow,  besides  the  large 

portion  (length  of  the  parallel  part.  The  nominal  length  of  the 

probe  was  L  +1^  ^ .  The  perpendicular  probe  (length  Lu>  did  not  have 
any  parallel  component  to  the  flow.  The  radii  of  both  probes  were  the 
same.  The  three  methods  were: 

(1)  Collisionless  case:  From  Eq .  2.68, 


IjA  La 
In  / 1  Lii  +  L: 


,  ,  ’  (4.38) 

2cL2H  TL*2-l 

where  L  =12  mm,  L  =56  mm,  L0  =4  mm  and  r  =0.05  mm. 

12  2  p 

2  2 

(2)  Collision  dominant,  this  sheath  case  (R^ot  X  <£.  1)  :  From  Eqs .  2.71 
and  2.77, 


ii 


L 


IJ. 


=0.B< f  . 


(4.39) 


.2.,  2. 


I||/2  yjf  L211  +  L21 

(3)  Collision  dominant,  moderate  sheath  case  (RgOOC^^l) :  From  Eqs.  2.74 
and  2.81, 


/  Ii  \  _ 5VS  L  a 

\  In  k  ^6  rp^+L^A4- 


^  +S3L2J.  (4.40) 

The  aim  of  the  second  measurements  was  to  check  the  electron  density 

obtained  by  the  double  probe  method.  Cases  1-4  of  Table  4.2  are  for  long 

ion  mean  free  paths  and  cases  5-11  are  for  short  mean  free  paths.  Since, 

in  all  the  cases  the  ratio  r  /A.  were  not  sufficiently  different  from 

p  i  J 

unity,  electron  densities  were  determined  by  using  three  different  methods 
(1)  Lam's  theory  for  collisionless  plasma  in  which  case  n  are  denoted 


by 


n  (2)  Collision  dominated  theory  with  Eqs.  2.63  and  2.64,  with 


n 


(2) 


(3)  Collision  dominated  sheath  with  electrical  boundary  layer 

(3) 


theory  of  Eq .  2.74,  with  n 


For  large  \i ,  n^1)  agrees  with  n  ^  but  not  with  n 


(3) 


,  since 


k 


. 
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does  not  take  into  account  the  current  due  to  ambipolar  diffusion. 


e 

This  makes  the  calculated  density  larger  than  the  true  value.  However, 

for  short  A  ,  it  is  n  that  does  not  coincide  with  other  n  's. 

i  e  e 

Electron  densities  determined  by  the  conductivity  measurements  ,  shock 
relations  and  Saha  equations  agree  with  n  for  large  A_^  or  n^^  for 

small  ,  as  long  as  M  is  greater  than  8.  With  M<C8  the  plasma  density 
is  too  low  for  the  Saha  equation  to  hold,  i.e.,  density  loss  by  diffusion 
to  the  wall  becomes  greater  than  recombination  loss.  In  addition,  the 
temperature  of  the  shock  heated  gas  is  too  low  for  further  ionization  to 
occur.  Therefore,  for  the  determination  of  electron  density,  the  shock 
relations  are  no  longer  valid  and  the  driving  gas  alone  determines  the 
electron  density. 

Figure  4.14  shows  the  plot  of  the  electron  density  against  time 
measured  from  the  initiation  of  the  driver  discharge  (t  )  for  the  cases 
1-4  of  Table  4.2. 

The  plasma  density  was  measured  by  the  microwave  cavity  method 

in  order  to  check  the  electrostatic  probe  method.  The  cavity  was  set 

at  230  cm  from  the  driver.  With  V  =2.5  KV  and  p  =0.1  mmHg,  we  obtained 

f=0.8  MHz  and  £/£n=2.03,  which  yielded  that  n  =6.3x10^  cm  ^ .  For 

u  e 

8-1  9-1  10  _i 

this  plasma  V  =8x10  sec  ,  U  =4.5x10  sec  and  o)  =  l. 5x10  sec  ,  which 

m  p 

satisfie  the  conditions  that  1«  1  and  r)f4Cl  of  Eq .  3.2.  The  measured 
point  thus  obtained  is  also  plotted  in  Fig.  4.14. 
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FIG.  4.14 


Plasma  density  decay  with  respect  to  time. 
p^=0.1  mmHg  and  Vg=2.5  KV. 
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4.3.5  Electron  Temperature  Measurements 

The  electron  temperature  of  the  argon  plasma  with  Vc=2.5  KV,  p-^= 

0.1  mmHg  and  z=80  cm  from  the  driver  was  measured.  The  results  for 
two  kinds  of  electron  collecting  probes  (Sec.  3.2.4)  are  shown  in  Fig.  4.15 
The  dotted  lines  were  taken  without  a  dummy  probe  and  the  solid  lines 
were  taken  with  a  dummy  probe.  (4V  was  subtracted  from  the  values  of 
dotted  lines  to  obtain  the  solid  lines.)  The  slopes  of  the  curves  yield 
2.0  eV  for  the  0.1  mm  diameter  probe  and  2.1  eV  for  the  1  mm  diameter  probe 
The  electron  density  cannot  be  found  because  of  the  ambiguity  of  the  probe 
area,  since  the  probe  surfaces  were  badly  contaminated  after  a  few  shots. 
Note  that  the  floating  potential  has  a  larger  value  (-20  volts)  than  the 
expected  value  by  calculation  with  Tg=2  eV  (9.5  volts),  because  at  the 
floating  potential  the  probes  collect  only  the  fast  electrons  as  mentioned 
in  Sec .  3.2.4. 

The  speed  of  the  plasma  at  the  probe  position  was  0.75  cm/ps  (M=23.4) 

from  the  photodiode  measurements.  From  the  shock  relations  it  is  predicted 

that  T  =1.4  eV  and  n  /n,=10.  This  result  is  compared  with  the  conductivity 
e  2  1 

measurements  in  Table  4.3.  The  mean  free  path  between  electrons  and 

-3 

neutrals  becomes  4.4x10  m,  which  is  much  larger  than  the  probe  radius. 

The  effect  of  the  plasma  flow  can  be  neglected  as  mentioned  in  Sec.  2.5. 

The  electron  temperature  deduced  from  the  shock  relations  is  30%  lower 
than  that  measured  by  electrostatic  probes.  This  discrepancy  is  considered 
to  arise  from  the  following  factors;  (1)  the  shock  relation  is  not 
exactly  correct  since  the  driving  plasma  slug  still  has  a  large  effect 
on  the  temperature,  (2)  the  probe  current  due  to  faster  electrons  is 


not  negligible. 
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4.3.6  Conductivity  Measurements 

The  plasma  conductivity  was  measured  by  two  methods  mentioned 
previously  (two  coil  method  and  Lin's  method)  and  the  results  are  listed 
in  Table  4.3.  The  calculated  values  of  conductivity  as  well  as  the 
plasma  temperature  from  Mach  number  and  Shock  relations  are  also  listed. 
From  the  table  we  can  find  that  the  conductivity  and  the  temperature  can 
be  determined  from  the  shock  relations  when  M  is  greater  than  8.  With 
M  less  than  8,  driver  gas  rather  than  the  shock  heated  gas  dominates  the 
plasma  conditions  and  soCT  and  hence  T  is  higher  than  expected  from  the 
shock  relations . 


4.3.7  Plasma  Luminosity  Measurements 


The  pictures  of  the  plasma  were  taken  with  the  image  converter 
framing  camera  to  study  the  overall  luminosity.  Typical  pictures  are 
shown  in  Fig.  4.16  when  the  initial  pressures  were  0.1,  1  and  10  mmHg 
and  the  driver  capacitor  bank  voltage  was  3  KV.  From  the  pictures  the 
following  can  be  deduced: 


Pi 

(mmHg) 

z 

(cm) 

M 

(cm) 

M1 

h 

(cm) 

vsx 

(cm/us) 

Re 
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where  L  is 
P 


the  length  of  the  driving  plasma  slug,  L  is  the  length  of 
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TABLE  4.3  Some  results  on  the  plasma  conductivity  and  temperature  measurements. 
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C 


FIG.  4.16  Framing  pictures  of  the  plasma  produced  by  the  electro- 
magnetically  driven  shock  tube.  Time  increases  from  top  to  bottom. 

Each  picture  has  an  exposure  time  of  20  ns  and  the  times  separating  each 
consecutive  picture  are,  10,  12  and  14  ps  for  A,  B  and  C,  respectively. 
The  plasma  proceeds  from  right  to  left.  The  camera  sees  a  length  of 
18.5  cm  of  the  shock  tube.  Initial  pressures  for  A,  B  and  C  are  0.1, 

1  and  10  mmHg.  Other  experimental  conditions  are  listed  in  the  main 


paper . 
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the  shock  produced  plasma  and  v  is  the  expansion  speed  of  the  shock 

sx 

produced  plasma.  It  was  also  evident  that  the  expansion  speed  of  the 

driving  plasma  slug  was  very  small  (i.e.,  L  was  constant),  which  implies 

P 

that  the  snow  plough  rate,  k,  of  Eq .  4.34  was  very  small  (<10.1).  The 
shock  heated  plasma  could  only  be  observed  at  p^=0.1  mmHg.  The  shape  of 
the  contact  surface/shock  front  was  not  always  planar. 

The  length  of  the  driving  plasma,  ,  can  also  be  determined  from 
the  electrostatic  probe  signal  (or  photodiode  signal).  This  was  measured 


with  V  =2.5-4  KV,  p  =0.05-1  mmHg  and  z=50-250  cm  from  the  driver.  The 

! 

results  also  indicated  that  L  was  independent  of  V  ,  pn  and  z. 
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CHAPTER  V  THEORETICAL  CONSIDERATIONS  ON  THE  BREAKDOWN  PHENOMENA 


5.1  EFFECT  OF  THE  SHEATH  SURROUNDING  THE  ELECTRODE 


5.1.1  Floating  Potential  in  the  Collision  Dominated  Plasma 


The  potential  difference  across  the  cathode  sheath,  V  ,  is, 


V  =V,+V., 
p  d  f  * 


(5.1) 


where  is  the  voltage  applied  across  two  floating  electrodes  and  is 
the  floating  potential  of  the  surrounding  plasma  (Sec.  2.7).  Here  Vc  was 
assumed  to  be  zero.  In  an  ion  sheath  around  the  cathode  for  our  experiments 
electrons  experience  no  collisions  with  neutrals  at  all  while  ions  collide 
with  neutrals.  Therefore  the  floating  potential  can  still  be  expressed 


by  Vf  of  Sec.  "List  of  Symbols",  i.e.. 


(5.2) 


where  j  is  the  electron  random  current  density  expressed  by 


je  =  rieef 


ZKme  . 


(5.3) 


j  is  the  ion  saturation  current  density  which  can  be  expressed  by  Eqs 
i 

2.58,  2.59,  2.63,  2.71,  2.74  or  2.76  depending  on  the  plasma/electrode 
conditions . 


for  four  important  cases  will  be  calculated. 

(1)  Ions  are  collision  free  in  the  sheath:  From  Eqs.  5.3,  5.2  and  2.58, 


0; 


(5.4) 
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(2)  Ions  are  collision  dominated  in  the  sheath  and  the  sheath  is 
stationary:  From  Eqs .  5.2  and  2.63, 


Vf  «»-*?. 

O  O 

(3)  Collision  dominated,  flowing  plasma,  Rj(  X  4^1: 

Hi 


and  2.71, 


Vf  =  8.8^. 


(5.5) 

From  Eqs .  5.2 

(5.6) 


(4)  Collision  dominated, 
and  2.74, 


Vf=6.7^-. 


flowing  plasma. 


RjO(2X2»1: 


From  Eqs.  5.2 


(5.7) 


5.1.2  Ionization  of  Neutral  Particles  by  Electron  Impact 

The  values  of  the  Townsend  first  ionization  coef ficient , 0(  ,  defined 

by  the  relative  increase  in  the  electron  density  per  unit  length, 

dn  / n  dr,  can  be  approximated  as  (5.1), 

0  ^ 

_<  -2e 

y=Ae  e  .  (5.8) 

The  constants  A  and  B  can  be  determined  experimentally.  For  argon.  Meek 
and  Craggs  (5.2)  have  given, 

A=13.6  (cm  ^  mmHg  ^) ,  B=235  (V/cm-mmHg), 

while  von  Engel  (5.3)  has  given, 

A=12  (cm--*-  mmHg-^),  B=180  (V/cm-mmHg), 

both  for  E/p=100-600  (V/cm-mmHg).  ot  is  given  in  cm--*-,  p  in  mmHg  and  E 
in  V/cm.  At  higher  E/ p ,  oC  is  given  by  (5.4), 
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The  first  ionization  coefficient  is  only  meaningful  when  there  are 
enough  elastic  collisions  between  electrons  and  neutrals  so  that  the 
electron  mobility  can  be  defined;  with  the  mean  free  path  of  electrons 
comparable  or  smaller  than  the  distance  between  the  cathode  and  sheath 
edge  (anode)  ,  OC.  cannot  be  defined.  In  such  a  case  the  collisional 
phenomenon  has  to  be  treated  microscopically. 

Let  us  denote  the  electron  energy  by  k .  If,  at  x=x  in  a  rectangular 
coordinate,  electrons  have  an  energy  distribution  of  f(x,k),  then  the 


number  density  of  electrons  created  at  x=x~x+dx  is, 


dn?  =n0(%)dxf 


(5.9) 


where  is  an  ionizing  collision  cross  section  (which  is  a  function  of 
electron  energy)  and  nQ  is  the  neutral  particle  density  at  x.  Normalization 

of  f(x,k)  can  be  done  by 

co 

f(x,k)dk  =  (5.10) 


i 

The  function  f(x,k)  follows  the  change  of 

Tjx®  W  +  Ik  lx  =  -^fw^CQe+Ox+Qi), 


(5.11) 


where  Q  and  Q  are  elastic  and  exciting  collision  cross  sections, 
respectively. 

For  electrons  which  are  emitted  from  the  cathode  and  reach  the 
sheath  edge  without  any  appreciable  collisions,  the  following  assumptions 
are  valid: 

(1)  Emitted  electrons  have  small  enough  energies  in  comparison  with  the 


(5.1)  T.  Kihara:  Revs.  Mod.  Phys .  24_  45  (1952) 

(5.2)  J.M.  Meek  &  J.D.  Craggs :  "Electrical  Breakdown  of  Gases"  Oxford, 

Clarendon  Press,  ppl08  (1953) 

(5.3)  A.  von  Engel:  "Ionized  Gases"  Oxford,  Clarendon  Press,  pp  181  (1965) 
A.  von  Engel  &  M.  Steenbeck:  "Electrische  Gasent ladungen"  Springer, 
Berlin  (1934) 

(5.4)  T. L.R„  Ayres:  Phil.  Mag.  4_5  353  (1923) 

D.Q.  Posin:  Phys.  Rev.  5J3  650  (1936) 

M.A.  Harrison:  Phys.  Rev.  105  366  (1957) 

A.  von  Engel:  "Handbuch  der  Physik"  Springer,  Berlin  21  504  (1956) 
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characteristic  voltage  of  the  system  so  that  f(0,k)  can  be  considered  to 
be  a  5-function;  i.e.,  all  the  electrons  are  assumed  to  be  emitted  from 
the  cathode  at  zero  initial  velocity.  (Actual  energy  distribution  of 
secondary  electrons  is  shown  in  (5.2 l)  where  it  can  be  seen  that  the 
secondary  electron  energy  at  the  cathode  surface  is  much  smaller  than 
the  ionization  potential  of  argon  atoms.) 

(2)  Electrons  experience  only  few  collisions  so  that 


=Q 

dx 

Electrons  gain  energy  from  the  electric  field  only, 
The  velocity  of  the  electrons  is  then  expressed  as 
It  follows  from  these  assumptions  that 


(5.12) 

so  that  k(x)=eV(x). 
ve(x)=(2eV(x)/me)^ . 


dne=n0neQi(V)dx. 


(5.13) 


In  a  cylindrical  coordinate  system,  ng(x)  must  be  considered  as  a 
number  of  electrons  per  unit  length;  i.e.,  by  denoting  the  number  of 
electrons  per  unit  length  at  r  as  Ng(r)  , 


or 


dNe=n0NeQ.(V)dr 

CnoOiCv)  dr 

NeW  =  Ne<rp)eirP 


(5.14) 

(5.15) 


Q 

n^  is  still  the  number  density  (in  m  )  of  the  neutral  particles. 


With  the  values  of  E/p  and  d  (interelectrode  distance)  for  the 


normal  breakdown  experiments  where  oi  can  well  be  defined, 

f[0Qidf  >  J^cXdr , 


£ 


(5.16) 

since,  in  such  a  case,  there  are  many  elastic  collisions  which  reduce 

dr  indicates  a 

larger  value  than  the  true  value.  However  for  our  experimental  conditions 


the  speed  of  electrons,  i.e.,  df  (r  ,k)/dr=^0 ,  and  jnQQ^ 

J0 


when  E/p  is  very  large  and  d  (=rg-rp)  is  very  small  with  the  constant 


ratio  of  Ed/p,  electrons  have  to  travel  a  large  portion  of  d  to  obtain 
energy  sufficient  to  ionize  the  neutrals.  Therefore, 
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4  d 

[n0Q?dlr  <  Mr, 

Jo  o  ^ 

d 

and  foCdr  is  larger  than  the  true  value.  The  crossover  point  when 

Jo 

d=l  mm  is  approximately  E/p=1000  V/cm-mmHg.  The  value  of  for  argon 
as  a  function  of  electron  energy  is  shown  in  Fig.  5.1  (5.5). 


5.1.3  BREAKDOWN  between  Cylindrical  Electrodes  in  Neutral  Argon 


Before  proceeding  to  the  study  of  electrical  breakdown  phenomena 
in  a  plasma,  the  breakdown  voltage  in  a  neutral  gas  between  infinitely 
long  cylindrical  electrodes  is  investigated.  By  denoting  the  distance 
between  the  centers  of  electrodes  as  d  and  the  radius  of  both  electrodes 
as  r  ,  and  by  taking  the  x  axis  as  the  line  connecting  the  centers  of 
the  electrodes  and  y  perpendicular  to  x,  with  the  origin  in  the  center  of 
one  of  the  electrodes,  the  potential  distribution  at  P(x,y)  can  be 
expressed  as  , 


Vtw~£ 


-t- 


U 


(*-4  +i<4)vf +/ 


(5.18) 


where  is  the  potential  difference  between  electrodes.  The  radial 

f)  Q  1 

component  of  the  electric  field  at  P  has  the  form,  (r=(xz+yz )^) , 


'W- 


T 


Ei”)" 


(5.19) 


If  d  is  sufficiently  larger  than  r  ,  x  and  y,  then  V  and  Er  at  P  can  be 
approximated  as , 


V(r)«— %-L-L 

il4-  h  ‘'p 


(5.20) 


(5.5)  S.C.  Brown:  "Basic  Data  of  Plasma  Physics"  M.I.T.  Press  (1966) 
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as  a  function  of  electron  energy  (5.5). 
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*4 

2r  H«4  • 

rp 


(5.21) 


as  , 


The  condition  of  breakdown  between  two  electrodes  can  be  expressed 


(5.22) 


^o(dr  c'ln(t+7jr)  =  §.  , 

rf 


(5.23) 


where  0(  and  T  are  the  first  Townsend  ionization  coefficient  and  the  rate 
of  secondary  electron  emission  at  the  cathode  per  positive  ion, 
respectively.  Usually  Y  is  a  weak  function  of  the  incident  ion  energy, 
and  so  4>  is  regarded  as  a  constant.  Since  the  electric  field  is  strong 
only  at  small  r,  the  breakdown  phenomenon  is  dominated  by  the  conditions 
near  the  cathode  and  E(r)  can  be  approximated  by  Eq .  5.21.  By  denoting 
the  breakdown  voltage  by  ,  we  obtain  for  the  condition  of  breakdown, 
from  Eqs .  5.8,  5.21  and  5.23, 


Vb 


2B  Pol  1* 


) 


(5.24) 


From  Eq .  5.24  the  pressure  pm  which  gives  the  minimum  breakdown 
voltage,  Vm  at  particular  values  of  d  and  ^  can  be  determined.  By 
differentiating  with  respect  to  p  and  putting  dV^/dp=0,  we  have, 


w  _  2B3> _  ^  _ 

A  L  Cd/rP )  «i/rp)  'UtCd/rp) 

£  (d/fp^-l  ld/rp)  -| 


(5.25) 


and 


Vm 

2pm(d-rp) 


(=constant) . 


(5.26) 


By  differentiating  Vm  again  with  respect  to  d/rp,  we  obtain  the  absolute 
minimum  breakdown  voltage,  Vmm  for  any  combination  of  p,  d  and  r  .  This 


c  : 


. i  , ?  . r ' f  . s  -i  -j s 


*i  veiI  ,0  qb\  Vb  gflidUnq  bos  q  03  3Deqes  Shiv  ,V'j;  r3i;  l: 


a  - 


.  i  fcj.B  b  ,q  3-  id*ni  na  v*6  -joi  ,eg*3lov  r  obafs®  id  mug  tnlm 


is  given  approximately  as , 
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(e=2,7i?---) 


(5.27) 


when  d-^Tp  and  p-^-ao. 

Equation  5.24  can  be  simplified  when  p  and  d  are  not  too  small  to. 


(5.28) 


This  equation  shows  that  all  the  data  can  be  reduced  to  one  line  if 
V^/ln(d/rp)  is  plotted  against  p.  This  line  corresponds  to  the  well 
known  Paschen  curve  of  vs  pd  for  parallel  plate  breakdown  voltage. 


5.1.4  Electric  Field  at  the  Electrode  Surface  in  the  Sheath 


The  intensity  of  the  electric  field  at  the  electrode  surface  (Ep) 
when  the  electrode  is  working  as  an  ion  collector  will  be  examined.  When 
the  electric  Reynold's  number,  Rg  (Eq.  2.69),  is  smaller  than  one,  the 
sheath  around  the  cylindrical  electrode  is  axially  symmetric  and  Ep  can 
be  found  by  solving, 

J=2Tcrni(r)evi(r)  =  -2TcrniepiE(r)  (5 .29) 

and 


=  -  f  IFte) 


n,-e 

£o  • 


(5.30) 


J  is  the  ion  current  per  unit  length  entering  the  sheath.  From  these 
equations  we  have, 

(5.3i) 

''.-se-MV./fSR  -fW J 


and 


(5.32) 


■v 


r  :  •'  •  (  .  i  .  J  ) 


■ 
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where 


F1-  J 

Kr\  ~ 


2tc£o^ 


(5.33) 


E  can  be  found  from  these  two  equations  by  eliminating  r  .  In  the  thin 
P  s 

sheath  limit  can  be  simplified  to, 


P  l2lc«oPirp;  ? 


(5.34) 


and  in  the  thick  sheath  limit, 

Vp 


EP=  /  Ps  and  Vp=f>Eh<S!«-r 


G 

p 


(5.35) 


When  the  plasma  flows  perpendicular  to  the  electrode  surface  with 

Rg>l ,  the  electric  field  in  the  sheath  is  deformed  from  a  cylindrical  to 

oval  shape.  However,  as  long  as X=eVp/kTe’^>  1 ,  the  direction  of  the  current 

in  the  sheath  is  radial  and  therefore  E  in  Eqs .  5.31,  5.33  and  5.34  is 

P 

still  valid  as  far  as  it  is  considered  to  be  an  average  value. 

The  current  per  unit  area  at  the  stagnation  point  is  larger  than  the 
average  current  density  by  a  factor  of  approximately  2ur„/2r  =1t,  and  so 


the  Ep  at  the  stagnation  point,  Epmax ,  is, 

<C  -  .  =|.46 


(5.36) 


for  the  thin  sheath  limit  and 


£  -  ~ 


I  - 


i.4 


(5.37) 


for  the  thick  sheath  limit. 


,2v2. 


The  exact  value  of  E  for  the  thin  sheath  and  R^of  X  »1  can  be 

pm  ax  & 

obtained  from  Eq .  2.73  as, 

l 


<<»-•> -<±fra4\ 


(5.38) 


■ 
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thus ,  £. 

c  _  iVp  _  /  3*  V,  J  V* 
2d(«.o)“(  4p,'£.rp  ) 


(5.39) 


or 


(5.40) 


5.1.5  Influence  of  Ionizing  Collisions  on  the  Electric  Field  in  the  Sheath 


When  the  neutral  density  in  the  sheath  becomes  high,  or  the 
potential  of  the  cathode  extends  far  from  it,  the  collisional  ionization 
of  neutral  particles  by  secondary  electrons  emitted  from  the  cathode 
becomes  large  and  the  potential  distribution  in  the  sheath  has  to  be 
modified . 

Around  the  electrode  immersed  in  a  plasma  flow  there  is  an  electric 
boundary  layer  in  addition  to  the  sheath  and  hydrodynamic  boundary  layer. 
Usually  the  thicknesses  of  these  boundary  layers  have  the  same  order  of 
magnitude  and  are. 


2rp//R^  for  the  electric  boundary  layer,  and 

2rp//Re  for  the  hydrodynamic  boundary  layer. 

Therefore  the  radial  position  of  the  electric  boundary  layer,  r  , 

a 

defined  by. 


la  =  rp  +■ 


22L 

JWe 


is 


(5.41) 


When  the  sheath  is  thicker  than  the  electric  boundary  layer, 
there  is  no  apparent  change  in  the  electric  field  due  to  the  boundary 
layer.  However,  when  the  sheath  becomes  thinner  than  the  boundary 
layer,  the  situation  is  different.  Most  of  the  electric  field  is  confined 


1  i  1  SB  B  i:  k  '  :  ;-  £  7  )  1 9  ari  f 
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to  the  sheath  region  as  in  the  case  of  the  thick  sheath  although  there 
is  a  weak  field  in  the  region  between  the  sheath  and  the  boundary  layer. 
In  this  region  the  densities  of  ions  and  electrons  are  the  same  so  that 
the  particle  movement  is  dominated  by  ambipolar  diffusion;  therefore 
let  us  call  this  region  an  ambipolar  diffusion  region.  All  the  ionized 
ions  in  this  region  come  back  to  the  sheath  and  affect  its  potential 
distribution . 

Here  we  consider  the  case  when  r  >r  .  For  r  <r  ,  all  r„  should 

cl  S  S  cl  ** 

simply  be  replaced  by  r  .  If  an  ion  current  per  unit  length  equal  to 

o 

J  enters  the  electric  boundary  layer,  then  the  number  of  ions  per  unit 
volume  at  radius  r,  n^(0),  is, 


(5.42) 


These  ions  eventually  hit  the  electrode  and  emit  TJ/e  electrons  per 
second.  These  secondary  electrons,  while  travelling  radially  towards 
the  sheath,  hit  neutral  atoms  and  ionize  some  of  them.  For  our 
experimental  conditions  the  thickness  of  the  boundary  layer  is  much 
smaller  than  the  electron  mean  free  path  and  the  assumptions  made  in 
Sec.  5.1.2  are  valid.  The  total  number  of  ions  per  unit  length  per 
second  created  between  r=r  and  r-fdr  can  be  written  from  Eq.  5.15  as. 


(5.43) 


These  newly  created  ions  begin  to  travel  radially  towards  the  cathode. 
The  total  number  of  these  secondary  ions  per  unit  length  per  second 


passing  through  r=r  is  ilSL  ,  therefore  the  secondary  ion  density  at  r  is, 


r 


(5.44) 
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9  •  .  ■  :  I 


(:  .  ) 


e  '■  i  T  ♦  ;  '  1  f  r  o  a  '  fl 


102 


or,  at  the  electrode  surface, 


n;(,)(rP)  =  — ^ - 


(5-45) 


where 


n0Q,*eJrp  dr. 

JrP 


(5.46) 


(1) 


These  ions  again  hit  the  electrode  surface  and  emi t  T2itrpn^ v  '(rp)v^(rp) 
=TJG/ e  electrons  per  second,  which  happens  to  be  G  times  the  initially 
emitted  electrons.  These  electrons  ionize  neutrals,  and  repeat  the  same 
procedure  again.  Finally  the  total  amount  of  ions  hitting  the  electrode 


is 


and  the  total  ion  density  at  r  is, 

j 


(5.47) 


(5.48) 


Before  solving  Poisson's  equation  we  express  the  ion  velocity  in  the 
sheath  by  v^ (r )= | p^E(r )  =-p^E(r)  (E<0) .  Then, 

W"T  -  TRSaO  -  *\s. 


where 


p  ^  _  J 


2.TC  €0^1 


49) 


(5.50) 


The  contribution  of  the  electrons  to  the  field  is  neglected  because  of 
their  large  mobility.  By  integrating  from  r=r  (E=E,  V=V)  to  r=rg 
(E=0,  V=0)  (5.6),  we  obtain. 


Br)±= (i^lrs>  - 


(5.51) 


(5.6)  This  assumption  of  E=0  and  V=0  at  the  sheath  edge  was  shown  to 

be  reasonable  when  eV  /kT  >?>  1  by  Schulz  and  Brown  (2.14). 

P  6 


c  :J0£«L9  fcriJl  'S 


>■(  p“  V  f:  1  t  v  *o«  n  y;  si 
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(5.52) 


(5.53) 


Since  the  exact  solution  of  these  equations  is  difficult  obtain,  we 
study  two  limiting  cases  depending  on  the  values  of  r  and  r  . 

cl  S 

(1)  r :  In  this  case  the  neutral  density  is  approximately 

equal  to  ^  (suffix  2  indicates  the  quantity  in  the  incident  shock).  As 

V  is  large  (much  larger  than  V. ,  the  ionization  potential  of  the  gas) 

P  1 

for  the  thick  sheath  case,  can  be  considered  to  be  constant  throughout 
the  sheath  (Fig.  5.1)  and 


fn0QiJr  =  n2Q.,(r-rp)  i 

rP 

2  0  2 

where  Q. =2. 8x10  m  .  E  also  takes  the  constant  value  of  E  o 

m  mz 


(5.54) 


Therefore 


G=1  is  the  spark  discharge  breakdown  condition.  Ep ,  Vp  and  rg  are 
computed  with  various  neutral  densities  and  the  results  presented  by 
the  thinner  lines  in  Fig.  5.2  in  terms  of  normalized  quantities,  E  /E^ > 
VrpEm2  and  d/rp=(rs/rp) -1 .  (The  notation  n£  is  changed  to  nw  in  the 
graph,  as  this  graph  can  also  be  used  for  the  case  of  r0^?>rc.)  r  has 

a  o  S 


the  limiting  value  of. 


(es.o 
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rs  =  rP  +  ^Ui^). 


(5.58) 


When  r  is  sufficiently  close  to  the  value,  the  ratio  of  E^/E^  and 

S  r 

Vp/rpEm2  becomes  constant. 

(2)  ra^>rs:  The  neutral  density  is  now  a  function  of  r,  as  seen 
from  Eqs .  5.106  and  5.107, 

13  ns(lc>)  rp) 


n(r)  = 


(rr{)r  -  * 

Since  the  electric  field  outside  the  sheath  is  very  weak,  all  the 

electrons  in  the  ambipolar  diffusion  region  can  be  considered  to  have 

an  energy  of  eV^.  Thus  the  ionizing  cross  section  takes  a  constant 

value  Q-^(Vp)  and  G  of  Eq .  5.46  can  be  calculated  as, 

.Is  rra  *1  M 

^=nif  +7ls  -7tP  Wr,  , 


(5.59) 


(5.60) 


where 


and 


$s  =  r(eniMQi<n~rf)-i) 

-r*  xaQiCndr 


(5.61) 


4a  =  rQi^n^e  V  dr  =  -i] 

=•  r-  i.3M5(3i(r9-rp)jLl((.3-S)< 


(5.62) 


In  the  sheath  the  neutral  particle  density  is  considered  to  be  a  constant, 

n  (suffix  w  indicates  the  quantity  at  very  near  the  electrode  surface), 
w 

Therefore  Eqs.  5.52  and  5.53  can  be  reduced  to  Eqs.  5.56  and  5.57,  with 
the  exceptions  that 

suffix  2  is  changed  to  suffix  w,  and 
G  is  changed  to  Gg+Ga. 

When  the  sheath  is  thin  or  the  neutral  density  is  low  enough,  the 
ionizing  effect  in  the  sheath  can  be  neglected  and  we  have,  from  Eqs.  5.52 


and  5.53, 


. 


<f  !) 


, 


. 
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and 


(5.63) 


(5.64) 


The  graph  of  E  /EL„  and  d/r  against  V  /r  E  is  shown  by  the  thicker 
&  r  p  mw  p  &  p'  p  mw  J 

l. 

lines  in  Fig.  5.2.  (The  (1-G  )2  term  is  included  in  E  in  the  graph, 

since  the  ionization  effect  in  the  ambipolar  diffusion  region  can  be 

interpreted  as  if  J  is  increased  to  J/l-G  .) 

a 

Especially  when  r  is  very  thin,  the  above  equations  can  further  be 

u 

simplified  to, 


p  _  Enn/  j  2<sl  __  Vp  _  / _ 3  J\/p 

tp  =  -/Ral  rp  ~  2  T  -(2^(Jiwrp(i-<5-a) 

and 

W  _  Eyniv  /i?  _  / _ 4  J  ^ 

p  VRU  S<f>  ^tt^piWrp(l-q-ay 


(5.65) 


(5.66) 


5.1.6  Secondary  Ions  Returning  to  Cathode  from  Outside  the  Sheath/Ambipolar 


Diffusion  Region 


Secondary  electrons  emitted  from  the  cathode  gain  an  energy  of 

about  eVp  at  the  sheath  edge.  As  they  leave  the  sheath  they  ionize 

neutrals  for  (V  -V.)/V.  times.  Since  there  is  no  electric  field  outside 

pi  i 

the  sheath,  these  newly  ionized  ions  diffuse 

away  isotropically,  and  0/1C  fraction  of  these  ions  happen  to  fall 

into  the  cathode  sheath,  where  0  =  sin"'*‘(rs/(rp+\e2))  .  The  total  number 
of  secondary  ions  returning  to  the  cathode  per  electron,  f,  is  thus, 

Effective  ion  current  to  the  cathode  is  therefore  increased  to. 


(5.67) 


r  t  •  s  ■;  ■<  if  £  :  ?  '  ’  •  >  o:  ' 

c  H{  ts. 
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(5.68) 


where  T  is  the  secondary  electron  emission  coefficient  at  the  cathode 


surface  per  incident  ion. 

5.1.7  Double  Sheath  Formation  on  the  Cathode 

If  the  cathode  emits  electrons ,  the  potential  distribution  is 
modified  by  the  space  charge  effect  created  by  them,  especially  near 
the  cathode  where  electrons  have  rather  slow  velocities.  This  modification 
will  be  discussed  quantitatively  for  a  collisionless  sheath  below  (5.7). 

Consider  an  infinite  plane  cathode  at  -Vp  «0)  and  a  sheath  edge  at 
zero  potential  at  a  distance  d  from  the  cathode.  At  the  sheath  edge 
there  will  be  a  zero  potential  gradient  and  a  finite  ion  current  j ^  per 
unit  area.  When  there  is  no  electron  emission  from  the  cathode,  j.  can  be 
expressed  by  (Eq.  2.44), 


(5.69) 


where  dQ  is  the  sheath  thickness  when  the  electron  emission  from  the 


cathode  is  zero.  Let  vg  be  the  velocity  of  the  electrons  at  any  point 
P  which  is  at  a  distance  x  from  the  sheath  edge,  ng  the  electron  density 
at  P.  The  corresponding  quantities  for  the  positive  ions  are  denoted  by 
the  subscript  i.  Then, 


(5.70) 


and  assuming  the  ions  have  unit  charge. 


%meve2=e(Vp+V) , 


(5.71) 


where  V  (<0)  is  the  potential  at  the  point  P.  Poisson's  equation  gives. 


.  i'  €>\  v  .  r  a  •  3r  vr  ...  olio©]  >  ©tr* - '  ©b<  rfrlBO  ©r!j 


.  '.t©c-  -r  ©/fa 


d  l/  _  _  SCn-,  -  ne) 

=  to 

We  may  eliminate  n  and  v  and  then  after  substituting 

i 


j  t  V  >»i  /  i  r  vP 


and 


v  -2L 
A  ct  1 
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(5.72) 

(5.73) 


we  obtain. 


j 2±_a_jL  ( 1  *  ^ 

d\*'  9  41  VJ?  JF?/* 


(5.74) 


Since  the  ion  current  is  limited  by  space  charge,  we  impose  the  condition 


dV/dx=0  at  x=0.  Integration  then  gives, 

The  ratio  d/dg  is  found  by  observing  that  V=-V  at  x=d ,  thus. 


(5.75) 


(5.76) 


The  effect  of  electrons  from  the  cathode  in  neutralizing  the  ion  space 
charge  cannot  cause  an  increase  in  for  this  is  fixed  by  the  plasma, 
but  manifests  itself  by  changing  the  sheath  thickness  as  expressed  in  the 
equation  above. 

According  to  Eq .  5.75  the  potential  gradient  at  the  surface  of  the 

cathode  is  proportional  to  (l-0()^  and  so  O^CX^l.  When  0(=1  the  potential 

gradient  at  the  cathode  is  zero  and  the  electron  current  as  well  as  the 

positive  ion  current  is  thus  limited  by  space  charge.  Maximum  electron 

current  is,  therefore,  expressed  by 

JL 

>w -ji(-si)  ,  <s-77> 

and  in  such  a  case  the  curve  representing  the  potential  distribution  in 
the  sheath  is  symmetrical  about  its  central  point. 


(5.7)  I.  Langmuir:  Phys .  Rev.  33.  954  (1929) 

F.W.  Crawford  &  A.B.  Cannara:  J.  Appl.  Phys.  36  3135  (1965) 
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The  intensity  reduction  ratio,  <£*,  of  the  electric  field  at  the 


cathode  surface  by  electron  emission  from  the  cathode  is, 


=  -0)'adf 


(5.78) 


The  numerically  obtained  plot  of  <£  and  d/ d^  against  o(  is  shown  in  Fig.  5.3. 

There  will  not  be  any  noticeable  change  in  electric  field  at  the 
cathode  surface  unless  ok>0 . 1 ,  or  in  the  case  of  argon  plasma,  jeX30  j^. 

If  emission  of  secondary  electrons  is  due  to  ion  bombardment,  then  j  =0.1j.. 
and  the  double  sheath  effect  can  be  completely  neglected.  However,  for 
breakdown  due  to  the  electron  field  emission,  this  condition  appears  to 
be  quite  severe  for  low  density  plasmas  with  a  small  cathode  (I*<10  mA) . 

For  an  argon  plasma,  the  electron  mean  free  path  is  much  longer 
than  the  ion  mean  free  path.  In  our  normal  experimental  conditions  ions 
are  collision  dominated  in  the  sheath,  but  not  electrons.  In  such  a  case 
ionic  charge  contributes  to  the  electric  field  more  than  in  the  collisionless 
case,  and  tends  to  approach  unity  and  again  the  double  sheath  effect  on 
the  cathode  electric  field  is  negligible. 


5.1.8  Ion-Electron  Collision  in  the  Sheath 


Since  electrons  are  emitted  from  the  electrode  which  is  in  a 
positive  sheath,  electrons  may  experience  collisions  not  only  with 
neutral  particles  but  with  ions.  In  this  section  the  collisional  effect 
between  electrons  and  ions  is  considered.  Average  momentum  transfer 
collision  frequency  can  be  calculated  by  using  Eq .  2.6  for  and  Eq .  5.115 
for  electron  velocity  distribution  function,  as. 


r 

•  '  )  • 

' 


Ill 


FIG.  5.3 


The  variation  of  the  sheath  thickness  and  electric 


field  at  the  cathode  surface  due  to  double  sheath  formation  as 

a  function  of  ion  and  electron  current  densities  (i.  and  i  )  to 

Ji  Je 

and  from  the  cathode. 
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6.03x10  n,Z  1/1 .  (5.79) 


This  equation  is  valid  only  when  the  drift  velocity  of  electrons  can  be 
neglected.  For  a  non -Maxwellian  gas,  A  of  Eq .  2.7  takes  a  different  form 
(the  concept  of  temperature  cannot  be  used) .  However  as  it  is  a  logarithmic 
term,  the  result  has  only  a  small  error.  Electron  mobility  in  the  ion 


sheath  is,  from  Boltzmann  equation, 

_ _ !_(jeO.§tjw  _  3ox*  Sc  me 

U  —  /-ao  i  ■—  1 1  -vi /  — 

'  \  -fvu)  Jv 


6SI  xio 
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(fcv)dv)meli»av  n-.Z^/l  - 


(5.80) 


(If  we  assume  a  Maxwellian  distribution  with  Tg=2  eV,  then  we  have, 

(4=5 . 85xl023/niZ2  In  A. ) 

In  the  case  where  the  sheath  has  an  extremely  high  electric  field, 
the  concepts  of  and  cannot  be  applied.  This  is  proved  as  follows. 

Let  us  consider  an  electron  leaving  the  electrode  at  initial  velocity  of 
Vq.  Before  it  collides  with  an  ion  for  the  first  time  at  t-^ ,  it  keeps 
increasing  in  velocity  as, 

V-Vo+^t.  (5.81) 

At  t=  t ,  the  probability  of  collision  with  an  ion  becomes  unity,  so  that 

rti 

we  have  \  \)  dt=l.  By  solving  the  equation  using  Eq .  2.6,  it  is  found  that, 

J0  m 


 WiVo 


J  HiZ* 


me  £E  v. 


r  -I 


(5.82) 


iZ*a3  in  A 

For  t-^  to  exist,  the  term  in  the  square  root  must  be  positive,  i.e., 

<  |.35  xio'S 


IT): 


In  our  typical  conditions  where  E=2xl0^  V/m  and  Vq=1  eV  (both  underestimated) , 


ns  >  5x(o2i  m'3, 


which  is  out  of  our  experimental  range.  Thus,  all  the  electrons  ejected 


' 
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from  the  electrode  "run  away"  and  suffer  no  collisions  with  ions. 

In  an  actual  case  the  electrons  cannot  acquire  infinite  velocity, 
since  they  are  decelerated  by  collisions  with  neutrals.  Again,  by  choosing 
underestimated  values  such  that  (electron-neutral)=  4  m^/V*sec  (at  100 

mmHg)  and  E=2xl0^  V/m,  the  electron-ion  collision  frequency  becomes 
0.94xl07  sec-'*'.  An  electron  proceeds,  between  each  collision  with  an  ion, 

f»eEA4r0.8 

meters,  which  is  many  orders  of  magnitude  larger  than  the 
sheath  thickness;  i.e.,  the  electron  runs  away  again. 


5.1.9  Current  Pinching  Effect 


When  there  is  a  current  through  a  plasma,  the  plasma  may  be  pinched 
by  the  magnetic  pressure  produced  by  the  current  (z -pinch).  The  strength 
of  this  force  can  be  expressed  by  the  ratio  <P>  of  the  plasma  pressure 
which  acts  against  the  pinching  force,  and  the  magnetic  pressure, 


c(flj+ne)kT  


2(*o 


1 


(5.83) 


where  c  is  the  compression  ratio  of  the  plasma  density.  This  ratio  can 
be  defined  as  R  / rz ,  where  R  and  r  are  the  initial  and  final  radii  of  the 
plasma  column.  The  values  of  the  plasma  properties  which  yield  the  smallest 
j9are,  in  our  experimental  conditions,  R=10“^m(=diameter  of  the  electrode), 
T=lcA,  ne=l()16  m"3  and  1=0.01  A  with  ^9=10,  i.e.,  the  magnetic  pressure  is 
not  strong  enough  to  pinch  the  plasma  column.  Therefore  the  pinching 
force  acting  on  the  current  is  always  negligible.  Current  pinching  can 
only  be  formed  in  the  presence  of  some  other  mechanism,  e.g. ,  local 
heating,  non-uniform  insulating  surface  layer,  local  breakdown  of  surface 


layer . 
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3e)  •  ■  -  ■  B  •  V  A  (l£  f.(  r ■!  >‘  If >)  .  3;-  O  18  89 i  ■'  I  ■  iBl  "Jc919! : tf 


. 


y<  booLfco^q  9  ::S  .vc  Brrr  9xf3  yd 


xf-  I  _ry  -O  J  <  1  j3b  J' 


114 


5.2  IONIC  MOBILITY  NEAR  THE  ELECTRODE 


5.2.1  Mobility  of  Ar+  in  the  Parent  Gas 


The  mobility  of  the  argon  ion  at  room  temperature  has  been  measured 

o 

by  many  people  (5.8)  and  the  value  is  consistently  1.6  cm^/V*sec  at 
760  mmHg  ,  0 °C  or, 

pi =0.12  m^/V-sec  at  1  mmHg  and °0  C.  (5.84) 

Theoretical  calculations  of  argon  mobility  (5.9)  support  this  result  well 
Unfortunately  ionic  mobility  at  high  temperature  has  not  been  measured 
except  for  N0+  in  air  (5.10),  where  it  has  been  found  that  the  mobility 
decreases  with  increasing  temperature  slowly  (slower  than  T-^).  Mobility 
is  strictly  inversely  proportional  to  neutral  density. 

In  our  calculations,  the  ionic  version  of  Eq.  2.3, 


u-3Pj;  ljoo_ 

H  0A2  yT.CK) 


(rnYv-sec) 


(5.85) 


will  be  used  for  the  mobility  at  high  temperature.  When  calculating  the 
electric  field  at  the  electrode  surface  from  Eqs .  5.34,  5.66,  etc.,  Ep  is 
approximately  proportional  to  nm  »  i.e.,  Ep  is  almost 

independent  of  T^ .  Therefore  the  approximation  of  Eq .  5.85  is  quite 
reasonable . 


(5.8)  R.M.J.  Munson  &  A.M.  Tyndall:  Proc .  Roy.  Soc .  A117  187  (1940) 
M.A.J.A.  Hornbeck:  Phys .  Rev.  84  615  (1951) 

M.A.  Biondi  &  L.M.  Chanin:  Phys.  Rev.  94.  910  (1954) 

(5.9)  H.R.  Hasse  &  W.R.  Cook:  Phil.  Mag.  12_  554  (1931) 

T.  Holstein:  J.  Phys.  Chem.  56_  832  (1952) 

E.A.  Mason  &  H.W.  Schamp ,  Jr.:  Annals  of  Phys.  4  233  (1958) 

(5.10)  J.K.  Dukowitz:  AIAA  Journal  8  827  (1970) 
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5.2.2  Stagnated  Shock  Wave 


In  the  previous  sections  it  was  found  that  the  knowledge  of  neutral 


When  an  electrode  is  surrounded  by  a  hypersonic  plasma  flow,  a  stationary 
shock  wave  is  formed  around  the  electrode  as  seen  in  Fig.  5.4.  If  we 
denote  two  regions  which  are  separated  by  the  shock  wave  as  2  and  5,  as 
in  the  figure,  then  the  gas  in  the  region  5  is  stagnated  by  the  electrode, 
which  results  in  a  further  increase  in  the  temperature,  pressure  and 
density  of  the  gas.  It  is  more  convenient  to  derive  these  quantities  in 
the  stagnated  shock  region  in  terms  of  the  conditions  in  the  incident  gas, 
rather  than  in  terms  of  the  initial  conditions  where  the  gas  is  stationary. 
The  properties  thus  derived  may  then  be  related  to  the  initial  conditions 
by  way  of  Mach  number. 

For  an  ideal  gas  with  constant  specific  heats,  the  continuity, 
momentum  and  energy  equations  relating  the  flow  variables  on  both  sides 
of  a  shock  wave  can  be  written  as  (5.11), 


and 


i°2^2n  ~ fgUsn 

7  Pa  (  Ui  _  7  fir  ( 


7-\  (V  2  “  7-1  ps  2  ’ 


(5.86) 


by  using  the  notation  introduced  in  Fig.  5.4.  We  can  write  also. 


U*n  -  U2.  s\n  £ 

Uzt  =  U2.  COS  5 

Usn  =  Us  svn(§-J) 
u$t  =u5  cos  (€-<£)  . 


(5.87) 


(5.11)  J.A.  Owczarek:  "Funcamentals  of  Gas  Dynamics"  Int'l  Textbook  Co., 
pp549  (1964) 
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FIG.  5.4 


The  formation  of  a  stagnated  shock  wave  around 


the  stationary  cylindrical  electrode. 


cylindrical  electrode 


.  . 
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The  angle  £  can  take  values  between  90°  and  sin"^!/!^). 


By  solving 


these  equations  we  obtain, 

fg  .  4*1 


or- 1  + 


a 


Ps 


2  nr 


(5.88) 


(5.89) 


where 


M,= 


(r 

9<y  z 

^j(rlx5ini)  - 


fqq  (Mis-mjfcotg  f 

(M2  «!«§)+• 


Uj 


(5.90) 


(5.91) 


(5.92) 


j)^y  and  take  their  maximum  values  at  £=90°  ,  when  the  above  equations 
are  reduced  to  Eqs .  4.15,  4.16  and  4.17  if  the  subscripts  1  and  2  are 
changed  into  2  and  5,  respectively  and  M-^into  M£  . 

There  is  a  limiting  value  for  M2.  From  Eq .  4.15-4.18, 


M*= 


2  (Mi2—  1) 


< 


Jii'f-O M|1+2)(2TM,2'- C'T-0)  J  2 7"fT-0 


(5.93) 


For  argon  plasma  with  M-^>5  and  £=90°  ,  M2  reaches  the  limiting  value  of 
3A/5*,  with 


/V 

A  Pz  ■  ’  rz  S  >  •  •»  *  Cz 

In  Eq.  5.86  it  was  assumed  that  further  ionization  in  the  stagnated  shock 


*  * 1  c'  1  *  -.9  — Tk  _  N/i  _  ^ 

P-  S  »  P2”Z  ’  Tz  3  ?  ? 


Usr  2JF  , 

S’  • 


was  negligible,  since  the  relaxation  time  for  argon  gas  in  the  stagnated 


shock  is  much  longer  than  r  /v,  the  time  spent  by  particles  to  pass  by  the 


.94) 


electrode . 
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5.2.3  Neutral  Particle  Density  in  a  Boundary  Layer 


A  cylindrical  obstacle  in  a  gas  flow  has  a  surrounding  boundary 
layer.  If  we  choose  the  coordinate  system  as  in  Fig.  5.4  (x  direction 
coincides  with  the  surface  of  the  electrode  and  y  direction  with  the 
normal  direction  of  the  surface)  ,  then  the  gas  has  no  motion  at  the  point 
x=y=0,  which  is  called  stagnation  point. 

By  using  this  coordinate  system  the  temperature,  pressure  and 
the  density  distributions  in  the  boundary  layer  can  be  found  by  solving 
the  Navier-Stokes  equations.  In  this  section  we  shall  examine  these 
distributions  on  the  y  axis  only  (x=0 ,  d/dx=0)  ,  i.e.,  examine  the  flow  of 
normal  incidence  to  the  object.  The  following  notations  will  be  used: 
u  =  particle  velocity  with  respect  to  the  obstacle 
^  =  viscosity  (shear  viscosity ,oc  T^) 

V  =  %!  f>  =  kinematic  viscosity 

S  =  L/JTe  =  thickness  of  velocity  boundary  layer 
=  JkL/ c^u  =  thickness  of  thermal  boundary  layer 
L  =  characteristic  length  such  as  electrode  diameter 
k  =  thermal  conductivity 

c  =Tk/(rT-l)m  =  constant  pressure  specific  heat 
P 

Re  =  uL/^>  =  Reynold's  number 

Pr  =  Cp^/k  =  Prandtl  number  (independent  of  temperature,  Pr=0.67  for  Ar) 


Firstly,  the  flow  velocity  change  along  the  y  axis  is  expressed  as, 


U  =  -Jap  W) , 

where  <fX<)  is  the  solution  of  the  differential  equation 

a"1  aj.»  j/2  . 

<f  +<f><p  ~<p  +l=o, 


(5.95) 


(5.96) 
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and  £  = Ja./yy  where  a  is  a  constant.  This  equation  was  solved  numerically 
by  Howarth  (5.12)  and  the  results  shown  in  Fig.  5.5.  If  y  is  large  or 
viscosity  is  small,  u  can  be  approximated  simply  as, 


u=-ay . 

(5.97) 

The  displacement  thickness,  ,  is, 

<J,=  0.648j|. 

(5.98) 

The  temperature  distribution  in  the  boundary  layer  around  the 


adiabatic  wall  is  found  to  be  (5.13), 


(5.99) 


The  subscript  CD  expresses  the  quantity  at  the  place  where  the  stagnation 
effect  is  negligible.  If  the  wall  is  heat  conducting,  then  by  introducing 
Tw,  the  wall  temperature, 


T_  Tv 
Too 


=  r!+('  T»)  Uoo^V 


(5.100) 


This  T(u)  curve  has  the  shape  of  Fig.  5.6.  With  these  values  of  u  and 
T  the  pressure  can  be  determined  from  the  Navier-Stokes  equation.  Especially 


for  the  pressure  on  y  axis , 

dP_  =  _pUM  +2-/  an 
dy  f  9 y  vyKw)' 

The  approximate  solution  to  this  equation  is, 

-P^-zK-Ki+s^). 


(5.101) 


(5.102) 


(This  equation  is  valid  for  a  cylindrical  coordinate  system, 
stagnation  point,  we  have. 


too.)  For  the 


(5.103) 


(5.12)  L.  Howarth:  ARC  R  M  1632  (1935) 

H.  Schlichting:  "Boundary-Layer  Theory"  McGraw-Hill,  6th  ed .  pp88  (1968) 

(5.13)  H.  Schlichting:  ibid.  pp316 
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FIG.  5.6  Gas  temperature  (T)  vs  flow  FIG.  5.5  Flow  velocity  (u)  vs  distance  from 
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CN  *-0  *-0  o 

^  d 


With  the  knowledge  of  the  temperature  and  pressure  we  can  immediately 
deduce  the  gas  density  from  the  relation  of  p=nkT .  The  Prandtl  number 
which  appears  in  Eqs .  5.99  and  5.100  is  expressed  by  the  square  of  the 
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ratio  of  the  thicknesses  of  the  velocity  and  thermal  boundary  layers , 


(5.104) 


Our  experimental  conditions  will  be  considered  for  the  rest  of  the 
section.  The  typical  properties  of  the  gas  in  our  experiments  (at  M=18  , 
T=1 .2  eV)  of 


^=2.2xl0“4  (Kg/m-s) 

(Fig.  4.6) 

L=lxlO'4  (m) 

(electrode  diameter) 

m=40xl. 67x10  (Kg) 

(argon  atom) 

n=3xl023  (m"3) 

(neutral  density  in  the 

stagnated  shock 

with  p^=l  mmHg) 

and 

u=5xl03  (m/s) 

yield 

Re=46  , 

&  =1.5x10  3  (m) 

and 

£t=1 .8x10  3  (m) . 

(5.105) 

With  this  neutral  density  and  wall  temperature  (300°K)  the  mean  free  path 

“6 

of  argon  is  8.4x10  (m)  and  more  than  88%  of  argon  atoms  suffer  at  least 

one  collision  in  the  thermal  boundary  layer  ( l-exp(  )=0 . 88 )  .  Thus, 

a  boundary  layer  can  be  defined  and  the  Navier-Stokes  equations  can  be 
considered  as  good  approximating  equations . 

M  is  Eq.  5.103  is  equal  to  in  Eq .  5.91.  Since  (Mach  number 
of  an  incident  shock)  is  usually  larger  than  5,  M^=3/5  and  pg  in  Eq .  5.103 
becomes,  pg/p00=1.3.  As  p(u)  is  a  slowly  varying  function,  we  can  say 


' 
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approximately  that  p=1.3p00  near  the  wall  of  the  electrode 
temperature  in  the  boundary  layer,  the  condition  of 

«  i— ~| 
c  loo 


For  the 


(5.106) 


shows  that  the  linear  distribution  is  a  good  approximation.  Also  by 
approximating  the  linear  change  of  u  with  respect  to  y,  T(y)  can  be 


expressed  as  , 


TCy)=^y 

JOQ  ~  7w 


y*T»  XttTvv 

Yoo  ~  /vv 


(5.107) 


The  average  neutral  density  can  now  be  obtained  as. 


\kTJ  k(Xo  -K 


i 


iJjw™ 


1*5  Poo  loO 


fcCli-W*"1  Tw  • 


Since  n^p^/kToo  ,  T(X)»TW  and  ln(T<p/Tw)~ln(12000/300)=3 . 7  , 

n<p 

The  density  at  the  wall  is, 

fljV  _  Pw/  Tw  _  .  ^  Too 

Pod  Tw  X 


(5.108) 


(5.109) 


w 


(5.110) 


' 
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5.3  EFFECT  OF  THE  ELECTRODE  SURFACE 


5.3.1  Field  Emission  of  Electrons 


When  an  intense  electric  field  acts  on  the  cathode  surface,  then 
an  electron  current  flows  out.  This  "field  emission"  electron  current 


can  be  expressed  by  the  Fowler-Nordheim  equation,  (5.14) 


—  $3 

r-rr  3  r*2  SlC Jftn  Y 

•uMf  _£_§_e  3hee  _ 

J  p+f  4TC 


6  M/cm2),  <5'm> 


where  E  is  the  intensity  of  an  electric  field  at  the  surface,  <p  is  the 
thermionic  work  function  and  p  is  the  energy  of  the  Fermi  level  inside  the 
metal  relative  to  the  bottom  of  the  conduction  band.  In  the  last  term 
of  the  above  equation  <p  and  p  are  involts  and  E  is  in  V/cm. 


For  tungsten  <p  =4 .5  V  and  the  exact  values  of  j  are, 


E  (V/cm) 

2xl07 

3xl07 

4xl07 

5xl07 

6xl07 

j (A/ cm2) 

4.7x10  ^ 

3.7x10 

l.lxlO4 

3.9xl05 

4.2xl06 

E  (V/cm) 

7xl07 

8xl07 

9xl07 

10xl07 

12xl07 

j (A/ cm2 ) 

2 . 3xl07 

8 . 7xl07 

2 . 5xl08 

5 . 7xl08 

2 . lxlO9 

If  the  temperature  of  the  body  is  very  high,  then  the  current 
emission  is  enhanced  (5.15).  This  is  called  thermo-field  emission  and 
the  approximate  current  takes  the  form  of, 


(5.112) 


(5.14)  R.H.  Fowler  &  L.Nordheim:  Proc .  Roy.  Soc.  A119  173  (1928) 

T.E.  Stern,  B.S.  Gossling  &  R.H.  Fowler:  Proc.  Roy.  Soc.  A124  699(1929) 
"Handbuch  der  Physik"  Springer,  Berlin  21_  188  (1956) 

(5.15)  T.H.  Lee:  J.  Appl .  Phys .  30  166  (1959) 

"Handbuch  der  Physik"  Springer,  Berlin  21_  191  (1956) 
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For  tungsten  at  E=4xl0^  V/cm,  j(T)  is  given  by. 
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T  (  °K) 

300 

1000 

2000 

3000 

j(T)/j(0) 

1.03 

1.5 

10 

144 

5.3.2  Electric  Field  Enhancement  due  to  Surface  Irregularities 


Microgeometry  of  the  irregularities  of  the  cathode  surface  increases 
the  electrical  field  by  a  factor  we  shall  call  i.e., 
actual  surface  field 


h 


(5.113) 


applied  field  # 

Shottky  (5.16)  has  shown  how  large  values  of  may  be  accounted  for  by 
a  superposition  principle.  For  example,  if  the  emitting  area  under 
consideration  is  really  a  hemispherical  boss  on  a  cathode  which  on  the 
scale  of  the  boss  is  effectively  an  infinite  plane,  it  is  easy  to  show 
that  ^  =3.  If  the  area  lies  on  a  second  hemispherical  boss  lying  on  the 
first  and  is  small  compared  with  it,  then  j8=3^=9 ,  and  so  on. 

Since  field  emission  varies  exponentially  with  E,  it  is  evident  that 
an  intensification  as  small  as  10 70  because  of  surface  projections  would 
cause  most  of  the  field  emission  current  to  flow  from  such  points.  This 
results  in  local  temperature  increase  and  enhanced  secondary  effects  which 
lead  to  instability.  Usually,  for  field  emission,  it  has  been  experimentally 
found  that  (5.17), 


(5.16)  Shottky:  Z.  Physi  14_  63  (1923) 

(5.17)  W.S.  Boyle,  P.  Kisliuk  &  L.H.  Germer:  J.  Appl.  Phys .  2_6  720  (1955) 
W.J.R.  Calvert:  Proc .  Phys.  Soc .  B69  651  (1956) 

C.  Germain  &  F.  Rohrback:  Vacuum  1^8  371  (1968) 

R.N.  Bloomer  &  B.M.  Cox:  Vacuum  18  379  (1968) 
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^9=3-  100  (5.114) 

and  the  electron  emitting  area  is  10~^^-10"^  cm^  for  reasonably  clean 
electrodes . 


5.3.3  Argon  Ion  Interaction  with  a  Tungsten  Surface 

An  argon  ion  with  a  kinetic  energy  less  than  about  1  KeV  incident  on 
a  metal  surface  undergoes  two  processes;  (1)  it  becomes  neutralized  by 
an  Auger  transition  a  short  distance  in  front  of  the  surface,  (2)  the 
resulting  gas  atom  transfers  momentum  to  one  or  more  metal  atoms  in 
predominantly  elastic  collisions.  During  both  processes  electron  ejection 
can  occur,  while  during  the  latter  neutral  metal  atoms  which  have  absorbed 
sufficient  momentum  may  leave  their  lattice  sites  producing  disorder, 
damage,  and/or  sputtering.  The  gas  atom  finally  either  comes  to  rest 
within  the  metal  (i.e.  becomes  trapped)  or  escapes  through  the  surface  with 
reduced  kinetic  energy. 

Argon  ion  trapping  efficiency  into  tungsten  at  20°C  is,  10"^  at  ion 
energy  of  72  eV  (experimental  threshold  of  trapping),  0.0007  at  100  eV, 

0.01  at  250  eV  and  0.04  at  500  eV  (5.18).  Sputtering  yield  for  tungsten 
bombarded  by  argon  ion  is,  0.001,  0.1,  0.2,  0.5,  0.7  and  0.8  atoms/ion  at 
ion  energy  of  50,  100,  150,  200,  250  and  300  eV,  respectively  (5.19). 


(5.18)  G.  Carter  &  J.S.  Colligon:  "Ion  Bombardment  of  Solids"  Heinemann 
Educational  Books  Ltd. ,  London  pp363  (1968) 

(5.19)  G.  Carter  &  J.S.  Colligon:  ibid.  pp321 


V:  .•  3  r  '  ■  c  .1  r,  ;  :  :  t 


• '  '  :  ■'  .1  .  .1  :  .  ,  (i  .  -  ) 


126 


5.3.4  Secondary  Electron  Emission  from  the  Metal  Surface 


Secondary  electrons  can  only  be  ejected  from  the  surface  of  a  metal 
by  a  positive  ion  if  the  sum  of  its  kinetic  (K^)  and  potential  energy  (eV^) 
exceeds  2e  times  the  work  function  <p  of  the  metal;  this  follows  from  the 
conservation  of  energy  and  the  fact  that  for  each  ejected  electron  another 
electron  has  to  escape  to  neutralize  the  positive  ion.  In  the  case  of  Ar+ 
falling  on  the  negatively  biased  probe  made  of  tungsten,  the  above  equation 
is  generally  satisfied  and  the  secondary  electron  emission  coefficient  'Y 
for  various  impinging  ion  energies  is,  (5.20), 


V.  (eV) 

10 

100 

1000 

'Y 

0.096 

0.095 

0.099 

It  is  important  to  know  the  energy  distribution  of  the  emitted 
electrons.  These  electrons  fall  into  two  distinct  groups; 

(1)  small  energy  of  1-3  eV  irrespective  of  primary  ion  energy,  and 

(2)  a  relatively  narrow  group  of  greater  energy  dependent  on  the 
nature  of  the  target  material  and  the  energy  of  the  incident  ions 
( energy  ~ KH-eV^ -2 e^)  . 

Mainly  electron  emission  is  caused  by  Auger  de-excitation  of  an  excited 
atom  (first  group)  and  the  theoretical  electron  energy  distribution  for 
this  group  was  obtained  by  Hagstrum  (5.21).  In  the  case  of  Ar"*~  interaction 
with  tungsten,  the  theoretical  electron  yield  can  be  approximated  by  a 
velocity  distribution  function  of 


A  -7  ,  -2.6lx/o“l2Vz 

f(v)=r|.<93x/C?  V/3e 


(5.115) 


(5.20)  H.D.  Hagstrum:  Phys .  Rev.  104  317  (1956) 

(5.21)  H.D.  Hagstrum:  Phys.  Rev.  96_  336  (1954) 
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In  the  presense  of  intense  electric  field  on  the  metal  surface, 
has  been  observed  to  change  its  value.  For  instance,  if  hits  steel, 
increases  to  about  twice  as  large  with  increasing  E  field  while  for  H+ 
and  steel  there  is  not  an  appreciable  change  in  If  (5 .22)  .  The  variation 
of  Of  for  Ar+-W  is  not  known. 

It  has  been  found  that  may  become  larger  with  ion  energy  of  1000  eV 
or  more  when  the  metal  surface  has  adsorbed  gas  layers.  This  can  be 
explained  as  follows.  Since  the  ion  has  to  expend  energy  in  penetrating 
the  adsorbed  layer  it  has  more  chance  of  losing  energy  to  atoms  near  the 
surface  of  the  metal.  The  electrons  are  thus  formed  nearer  to  the  surface 
and  their  diffusion  through  the  adsorbed  layer  is  easier  than  through  the 
corresponding  thickness  of  metal  so  that  their  escape  is  easier.  In  our 
case,  however,  the  voltage  applied  across  the  sheath  is  usually  the  order 
of  200  volts  and  so  the  gas  adsorption  does  not  have  much  effect  to  the 
breakdown. 


5.3.5  Ionic  Charge  Accumulation  in  the  Sheath 

If  the  charge  exchange  efficiency,  ^ ,  between  impinging  ions  and 
the  electrode  surface  is  not  unity,  the  ionic  charge  is  accumulated  and 
the  ion  density  increases  in  the  sheath.  If  we  denote  the  ion  density  and 
velocity  at  the  electrode  surface  and  sheath  edge  as  n£(rp),  v^Crp),  n^(rs) 
and  v^(rs),  respectively,  then  (for  a  planar  sheath), 

ni  (rs>vi  (rs>=*Lni  (rp)vi(rp>  (5.116) 

and  so  it  can  be  understood  that  the  ion  density  in  the  sheath  at 


(5.22)  E.W.  Webster,  R.J.  Van  de  Graaff  &  J.G.  Trump:  J.  Appl.  Phys . 
23  264  (1952) 
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equilibrium  is  l/l^  times  the  ion  density  when’ll.  The  equilibrium 
is  attained  exponentially  with  a  time  constant  of  2d//^”^  ,  where  d  is 
the  sheath  thickness  and  v\  is  the  mean  ion  speed  in  the  sheath.  (The 
time  constant  is  of  the  order  of  one  to  ten  nanoseconds  for  our  experimental 
conditions . ) 

The  charge  exchange  efficiency  of  Ar+  on  W  is  approximately  807o  (5.23) 
and  so  the  ion  density  should  be  considered  to  be  1.25  times  larger  than 
the  value  obtained  withf[=l.  If  W  is  covered  by  an  oxide  surface  layer, 

^  becomes  smaller,  and  the  density  becomes  higher.  In  the  course  of 
calculation  of  the  electric  field  in  the  sheath,  the  charge  accumulation 
effect  can  in  effect  be  considered  as  if  the  ion  current  flowing  into 
the  sheath  is  increased  by  a  factor  of  i/i- 


5.3.6  Adsorption  and  Chemisorption  of  Gases  on  the  Metal  Surface 

The  number  of  monolayers,  N,  formed  on  tungsten  is  (5.24), 

N=7*lo7pj==C RT  (5.U7) 

for  N<1.  Here  p  is  the  pressure  in  mmHg ,  A  is  the  molecular  weight,  T 

is  the  surface  temperature  and  is  the  heat  of  adsorption/chemisorption . 

-12 

For  10  mmHg  argon,  E^=1.9  Kcal/mol ,  I"Qexp(E(j/RT)  =  10  sec  (soj ourntime) 

and  so  N=6.4xl0-^  monolayers.  In  other  words,  6.4x10  ^  %  of  the  electrode 

9 

surface  is  covered  by  3.2x10  argon  atoms.  Therefore  the  physical 
adsorption  of  argon  is  negligible  except  when  T  is  extremely  low. 

(5.23)  R.  Lawson:  Private  Communication  (1970) 

(5.24)  G.  Lewis:  "Fundamentals  of  Vacuum  Science  and  Technology" 
McGraw-Hill,  pp45  (1965) 
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In  the  case  of  an  oxygen  molecule,  £^=194  Kcal/mol  and  so  even 

at  p=10“l®  mmHg,  N2s>l,  which  means  that  every  oxygen  molecules  hitting 

the  tungsten  become  trapped  by  chemisorption  until  a  monolayer  is  formed. 

If  the  system  is  evacuated  to  10"  mmHg,  the  density  of  oxygen  in  the 

system  is  6x10^  cm"^.  At  300°K,  2.7x10^  oxygen  molecules  hit  an  area  of 
o 

1  cm  per  second.  If  it  takes  10  seconds  to  fire  the  shock  driver  after 

ceasing  to  heat  the  electrode  (actually  it  takes  more  than  that),  2.7x10^ 

o 

oxygen  molecules  per  1  cm  hit  the  tungsten.  For  tungsten  trioxide  (WO3) 
monolayer  to  form,  about  10  oxygen  molecules/cmz  are  necessary  and  so 
after  10  seconds  tungsten  electrode  has  a  complete  oxide  monolayer.  After 
the  monolayer  is  formed,  decreases  and  the  rate  of  chemisorption 
decreases  exponentially.  It  is  hard  to  find  the  number  of  layers  formed, 
but  it  should  be  around  3-30  layers  (5.25). 

of  H2 ,  CO,  CO2  and  N2  on  tungsten  are  46,  100,  122  and  85  Kcal/mol, 
respectively  and  the  chemisorption  of  these  gases  can  be  neglected  in 
comparison  with  that  of  03.  Adsorptions  are,  of  course,  much  smaller 
than  chemisorptions . 

The  surface  structure  of  WO3  formed  by  the  chemisorption  of  oxygen  on 
tungsten  is  (5.26), 


WO2  is  unstable  at  the  room  temperature  and  is  not  formed. 


(5.25)  R.  Lawson:  Private  Communication  (1970) 

(5.26)  N.K.  Adam:  "The  Physics  and  Chemistry  of  Surfaces"  Oxford  (1941) 
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5.3.7  Cleaning  of  Electrode  Surface  by  Heating 

When  the  shock  tube  is  driven  by  an  intense  arc  discharge,  small 
parts  of  the  electrode  and  insulator  are  melted,  accelerated  downstream 
and  contaminate  the  surface  of  the  electrodes  (or  probes).  Besides  the 
exposure  of  the  electrode  to  any  gas  leads  to  an  adsorption  or  chemisorption 
of  gas  on  the  surface.  The  surface  layers  thus  made  form  an  electrical 
insulator  along  the  electrode  surface  and  affect  the  electrode  functioning 
in  many  ways ; 

(1)  Contamination  affects  the  area  of  the  electrode  surface.  The 
reduction  of  the  effective  area  by  1/10  or  more  was  experienced. 

(2)  As  will  be  discussed  in  Sec.  6.3,  the  electrode  currents  tend  to  have 
spikes  even  at  low  electrode  biases  (as  low  as  15  volts,  depending  on 

the  surrounding  plasma  density  and  temperature)  ,  since  on  the  surface 
of  the  contamination  are  accumulated  the  electrical  charges  which,  with 
sufficient  quantity,  result  in  electron  avalanches  inside  the  contamination. 

Most  of  this  surface  contamination  can  readily  be  removed  by  pre¬ 
heating  in  vacuum  the  electrode  to  a  sufficiently  high  temperature.  (It 
is  for  this  purpose  that  tungsten  was  chosen  as  an  electrode  material.) 
Tungsten  can  be  repeatedly  heated,  cooled  and  air  admitted  without  damage 
to  its  properties,  although  heating  renders  the  tungsten  brittle.  Espe 
and  Knoll  (5.27)  have  given  for  tungsten  wire  of  circular  cross  section 
and  of  various  diameters  the  current  required  to  bring  the  tungsten  to 
various  temperatures.  As  the  melting  points  of  tungsten  and  tungsten 
trioxide  are  3370°C  and  2130°C,  respectively,  the  heating  to  about  2230°C 
(1.3  A  for  0.1  mm  diameter  wire)  cleans  the  surface  perfectly.  The  life 


(5.27)  Espe  &  Knoll:  "Werks tof fkunde  der  Hochvakuumtechnik"  Springer, 
Berlin  (1936) 
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time  of  tungsten  at  2230 °C  is  200  hours  for  0.1  mm  diameter  wire  and 
150  hours  for  0.05  mm  diameter  wire  (5.28).  The  electrode  was  heated 
under  high  vacuum  before  each  run  for  about  one  minute.  The  electrode  of 
0.1  mm  diameter,  therefore,  can  be  used  for  at  least  1000  times  without  any 
reduction  in  its  diameter. 

There  is  another  advantage  in  heating  the  electrode,  i.e.,  smoothing 
its  surface.  At  about  a  quarter  of  the  absolute  melting  temperature, 
surface  migration  allows  the  atoms  on  the  rough  tip  to  rearrange  themselves, 
reducing  the  free  surface  energy.  Perfect  smoothness  is  achieved  after 
a  few  minutes  at  2/3  of  the  absolute  melting  temperature,  provided  that 
no  oxides  or  other  contaminations  still  stable  at  this  temperature  are 
present . 

The  hair  pin  type  electrodes  as  seen  in  Fig.  3.3  were  made  so  that 
they  can  be  heated  white  hot.  McLaren  et.  al.  (5.29)  reported  the 
successful  use  of  this  type  of  probe  but  with  shorter  tungsten  wires.  If 
the  probe  is  long,  it  disturbes  the  plasma  and  also  degrades  the  spatial 
resolution.  On  the  other  hand,  the  short  probe  has  more  ambiguity  in  the 
probe  area  because  of  the  small  unheated  part  of  the  tungsten  feed-through 
j  oint . 

5.3.8  Cleaning  of  Electrode  Surface  by  Ion  Sputtering 

When  the  electrode  is  surrounded  by  ionized  gas,  high  energy  ions 
hit  the  negative  electrode,  which  results  in  the  sputtering  of  oxide 

(5.28)  M.  Pirani  &  J.  Yarwood:  "Principles  of  Vacuum  Engineering" 

Reinhold  Publishing  Co.  (1961) 

(5.29)  T.I.  McLaren,  J.N.  Fox  &  R.M.  Hobson:  6th  Int'l  Conf.  Ionization 
Phen.  in  Gases  4  311  (1963) 
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surface  layer  on  the  electrode.  An  empirical  formulation  for  the  time 
during  which  N  layers  of  tungsten  are  destroyed  by  argon  ions  is, 

t  =  I0-8  4^-  <sec)  (5.118) 

vJ  r 

where  j  is  the  current  density  into  the  electrode  in  A/ cm  and  p  is  the 
ambient  pressure  in  mmHg.  By  substituting  our  typical  experimental 
conditions  of  N=10,  j=l  (A/cm  )  and  p=10  (mmHg),  we  obtain  t=0.1  psec. 

On  the  other  hand,  oxidation  of  the  surface  continues  during  this 
period.  This  rate,  however,  is  only  1.5x10”^  monolayers  per  sec  with 
shock  processed  oxygen  molecules  of  10^ °K  and  10-^  mmHg,  and  so  further 
oxidation  can  be  completely  neglected.  Thus  with  the  arrival  of  plasma  at 
the  electrode,  the  WO^  surface  layer  on  the  surface  is  destroyed  within 
about  1  jjsec  and  after  that  the  electrode  can  be  considered  to  have  a 
tungsten  metal  surface. 

After  a  long  period  of  sputtering  (a  time  long  enough  for  over  100 
layers  to  be  removed)  ,  there  appear  on  the  electrode  surface  many 
irregularities  which  have  a  strong  influence  on  the  electric  field  at  the 
surface  (Sec.  5.3.2). 


5.3.9  Breakdown  of  WO3  Surface  Layer 

Although  most  of  the  WO3  will  be  destroyed  by  ion  sputtering,  there 
still  may  remain  some  on  the  electrode  surface.  The  conductivity ,  cr  ,  of 

O  fa 

WO3  was  measured  to  be  5x10"°  mhos/m  at  an  electric  field  of  10  V/cm  in 
WO3  (with  higher  field,  the  conductivity  becomes  lower)  (5.30).  By  denoting 
the  current  density  to  the  WO3  layer  as  j ,  the  electric  field  in  the  WO3 

r\ 

can  be  expressed  as  j /<T.  If  j  exceeds  0.01  A/cm  ,  the  field  becomes  higher 


(5.30)  A.  Guntherschulze  &  H.  Betz:  Z.  Physik  1\_  106  (1931) 
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than  2x10^  V/cm  which  is  sufficient  for  field  electron  emission.  In 
such  a  case  the  spikes  may  be  observed  in  the  electrode  current  and  this 
may  sometimes  lead  to  the  full  path  breakdown. 


5.4  EFFECT  OF  THE  TEMPERATURE  OF  THE  ELECTRODE  SURFACE 


5.4.1  Thermionic  Electron  Emission 


Thermionic  electron  current  is  expressed  by  the  Richardson-Dushman 


equation , 


j  =4T2  e~  kT 


(A/cnf )  1 


(5.119) 


where  A  is  the  thermionic  emission  constant  and  <p  is  the  thermionic  work 
function.  A=70  A/cm^°K^  and<P=4.5  V  for  tungsten  (5.31)  and  A=5xl0^ ^A/cm^ °K^ 
and  (p=9 .22  V  for  WO^  or  W-0  film  (5.32). 


5.4.2  Basic  Theory  of  Hot  Cathode  Arc  Discharge 


In  the  case  of  an  arc  discharge  the  secondary  electron  emission 
constant  due  to  positive  ion  impact,'/,  must  be  replaced  by  which 

gives  the  number  of  electrons,  thermionic  and  secondary  by  bombardment, 
liberated  per  positive  ion  at  the  cathode.  Thus  we  can  write. 


(5.31)  D.A.  Wright:  Proc .  Inst.  Elect.  Engrs .  100  pt.  3,  125  (1953) 

(5.32)  V.S.  Fourenko:  "Handbook  of  Thermionic  Properties"  Plenum 
Press,  Data  Division  ppl29  (1966) 
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ry'_  Ie 
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(5.120) 


where  Ie  and  1^  are  the  electron  and  ion  currents  at  the  cathode.  If 
we  take  It  as  the  current  due  to  thermionic  emission  specified  by 
Richardson  equation  (Eq.  5.119),  the  total  electron  current  at  the  cathode 
is  then,  I^riH-It ,  and  the  total  current,  I,  is,  I=Ie+I^=(l+'7)I^+It .  From 
this  we  can  evaluate  rf/  as  , 


V=y+  {+qr 


it  1 


(5.121) 


If  a  current  Irv  is  emitted  from  the  anode  (or  in  our  case  if  I 


■0 

enters  from  the  cathode  sheath  edge),  then, 


0 


«s  rrs 

/  (<xc/r  v  -j  °<dr 

Li=Io+le(e'l>  -•  «i0e* 


/s 

-  ^r\ 
+  i(i-eJ*  ). 


(5.122) 


By  eliminating  I.  and  I  we  obtain, 

1  r*  e 

erp 

at  + 


i+t 


i  -‘rCefif— •*-»)  l-r(e(^r-i) 


I 


O  . 


(5.123) 


The  condition  of  breakdown  is  determined  by  I/Iq-M)^  or, 

£«dr 


L  = 


£  rP 


l  -T(e^0(c,r-i) 


1 1. 


(5.124) 


By  using  ^  of  E^.  5.121  this  condition  can  be  rewritten  as, 

-i)  =  i. 


(5.125) 


When  ionization  effects  in  the  sheath  are  negligible,  Eq .  5.123  can 
be  reduced  to  I=It+lQ  and  the  condition  of  breakdown  is  simply, 


i=it. 


(5.126) 


1  rfv 
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5.4.3  Self  Sustaining  Arc  Discharge 


To  see  how  a  self  sustained  arc  can  be  established,  we  will  examine 
the  simplified  case  where  the  energy  input  to  the  cathode  is  due  to  ion 
bombardment  (=I^V,  V=potential  difference  across  cathode  sheath)  and  the 
output  is  due  to  the  radiant  energy  emission.  More  precise  discussion  of 
the  energy  balance  will  be  done  later.  The  basic  mechanism  of  the 
discharge,  however,  does  not  change  from  this  simplified  case.  The 
cathode  temperature  can  be  determined  by  equating  the  input  and  output  as  , 
IiV=So'T4  , 

where  black  body  radiation  was  assumed,  o'  is  the  Stefan-Boltzmann 
constant  (=5.67x10“-^  W/cm^°K4)  and  S  is  the  area  of  the  cathode  surface 

involved.  The  Richardson  equation  takes  the  form, 

oUJ  .  ,,  i 

It=F(I;V)=SATV*T  -SA(§9  e 

ao.g 

-  3.03  X |07 JsTy  6  (5.127) 


If  ionization  effects  in  the  sheath  are  very  small, 

I«lo  +  F[loV*IV$dr]# 


(5.129) 


Equation  5.128  will  be  examined  further.  Let  us  put  the  left  hand  side 
of  the  equation  as  y^  and  the  right  hand  side  as  y2  .  y-^  and  y2  as 
functions  of  I  have  the  shapes  of  Fig.  5.7.  Fig.  5. 7 -A  is  valid  when  Iq 
is  small.  Initially  I  stays  at  1^  (which  is  a  stable  solution).  If  Iq 


r  f#BO  otra  I  '  (  V  SS8  oT 


136 


■FIG.  5. 7 -A  The  right  and 

the  left  hand  side  of  Eq . 
5.128.  The  solutions  are 
I-j^  and  I^. 


FIG.  5.7-B 


The  condition 


of  breakdown  for  Iq  (=Iob^ ' 
I  jumps  from  1^  to  1^. 


FIG.  5.9  V-I  character¬ 
istics  of  an  arc  discharge 
of  Eq.  5.128. 


■  . 


. 
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is  increased,  y£(I)  curve  floats  up  and  with  1q=Iqb  (Fig.  5.7-B)  I  jumps 
from  Ig(=I-^=l2)  to  Iq  (Iq  is  also  stable  but  I2  is  unstable).  Now  if  Iq 
is  decreased,  y^Cl)  curve  sinks  down  to  Fig.  5. 7 -A  with  I  staying  at  Iq, 
and  even  with  Iq=0  ,  I  does  not  become  zero.  This  indicates  a  self 
sustaining  discharge. 

The  relation  of  I  vs  Iq  is  shown  in  Fig.  5. 8 -A,  B  and  C  for  three 
cases  of  V. 

(1)  small  V:  With  increasing  Iq,  I(Iq)  curve  follows  0->E,  and  E->0  with 
decreasing  Iq.  I  is  increased  slightly  more  than  Iq  by  the  thermionic 
electron  current. 

(2)  intermediate  V:  If  Iq  is  increased  from  zero,  I  vs  Iq  curve  follows 

0->A->B  and  jumps  to  D,  then  to  E .  If  Iq  is  decreased  from  that  point, 

I  changes  as  E->D-»C  -»( jump)->A-»0  and  I  cannot  be  self  sustained.  However, 

10  2 

I  at  point  D  usually  exceeds  10  A/cm  ,  in  which  case  I  is  restricted  by 
an  external  circuit  and  is  practically  considered  to  be  a  breakdown. 

(3)  large  V:  If  Iq  is  increased  I  follows  0->B->( jump)->D-^E  and  if  it 
is  decreased,  E-»D->C.  At  point  C,  a  self  sustained  arc  is  being  formed 
and  the  condition  of  Eq.  5.125  is  satisfied. 

2  2 

Let  us  denote  the  values  of  I  and  y2  which  give  d  y£/dl  -0  when  Iq=0 
as  I^j  and  y2m*  Then, 


(5.130) 


and 
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FIG.  5.8  I  vs  IQ  of  Eq.  5.128 
at  various  values  of  V. 


ii  \ 

A.  Small  V,  'T  (e r?  -|)<l. 

No  breakdown  occurs  even  with  Iq=00 


Is 

B.  Intermediate  V,  /Y^rP  —  1^0 
Current  jump  occurs  at  Iq^.  This 
is  not  a  self  sustaining  discharge. 


since  I  becomes  zero  with  lg=0. 
However,  if  1^  is  larger  than  the 
current  available  from  the  external 


source,  then  this  can  be  considered 


as  an  actual  breakdown. 


C.  Large  V,  'T  -l)  =  l. 

Self  sustaining  breakdown  occurs  if 
Iq  is  increased  up  to  Iq^  and  then 
decreased  to  zero. 


.  .OJ'i  • 
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The  last  terms  of  these  equations  were  calculated  for  the  case  of  a  tungsten 
electrode.  The  condition  for  a  self  sustaining  discharge  is  that  Eq .  5.128 
has  three  solutions  when  Iq=0.  By  considering  that  F(I^V)  is  an  exponential 
function,  this  condition  can  be  replaced  by  the  condition  that  F(I^V)  becomes 
larger  than  I  at  some  value  of  I.  Therefore  as  a  sufficient  condition. 


>2m  -23Sx|0  >  Im  -  £|_£-J£oidr)y 
(l-e^V  >4;22x|0'4. 


(5.133) 


(5.134) 


Even  with  the  severest  conditions  that  O(=0.09  cm~^  (at  100  mmHg)  and 
rs-rp=5xl0"^  cm, 

V >  10  volts 

is  the  condition  for  a  self  sustaining  discharge.  This  is  always  satisfied t 

since  V  never  becomes  smaller  than  the  ionization  potential,  . 

The  V-I  characteristics  of  an  arc  discharge  of  Eq .  5.128  are  seen 

in  Fig.  5.9  with  Iq  as  step  functions  of  various  amplitudes  (not  impulses). 

5  o 

Since  Iy[  in  the  figure  is  the  order  of  10  A/cm  ,  the  practical  experimental 
range  lies  on  the  left  hand  side  of  1^  and  the  N  shaped  curve  cannot  be 
observed . 


5.4.4  Energy  Balance  at  the  Cathode 

The  cathode  temperature  can  be  determined  by  equating  the  input 
energy  to  the  output  energy  at  the  cathode,  where  (5.33): 


(5.33)  K.T.  Compton:  Phys .  Rev.  37  1077  (1931) 
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Energy  input  =  1(1 -f ) (aVp+Vi -f0)+PJ+PH+PFf PN 
Energy  output  =  If^+Pg+Pg+Pj-yf  Pp 

and  f  =  fraction  of  current  carried  by  electron  at  cathode 

a  =  accomodation  coefficient  (fraction  of  kinetic  energy  of  ion 
delivered  to  cathode,  a=0.9-l) 

V  =  cathode  potential  drop 


V. 

1 


=  ionization  potential  of  gas 


%  =  normal  work  function  (work  function  at  E=0) 

<Pi 


=  effective  work  function  (work  function  at  E  present,  for 


PT  = 


P  T  T  = 


H 


thermionic  electron  emission,  and  for  field  emission  <Pi=°> 

Joule  heating 

F(f(VVoM4)'(1+S)(  l-f)V^-Ve)  =  heating  by  energy  derived  by 
electrons  in  the  cathode  fall  space  and  indirectly  returned  to 
cathode 


ip  =  V.-lfj+L  =  heat  of  neutralization  of  ion 

L  =  heat  of  condensation  of  neutrals  of  cathode 

F  =  fraction  of  the  energy  acquired  by  unelectrified  carriers  from 

the  electrons  which  have  moved  through  the  fall  space,  which 

returns  to  the  cathode  as  by  radiation,  metastable  state,  etc. 

Ve  =  electron  mean  kinetic  energy  in  plasma 

Pp  =  radiation  heating 

Pjvj  =  heating  by  neutral  particles 

Pg  =  cooling  by  evaporation 

Pg  =  cooling  by  conduction  through  cathode 

Pp  =  cooling  by  gas  conduction  and  convection 

P  =  cooling  by  radiation. 

K 


. 
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5.4.5  Heat  Conduction  into  the  Metal 


and 

where 


and 


(5.135) 


(5.136) 


The  basic  equations  of  heat  transfer  are, 

H--A  WT 
-^  =  -avT+-^ 

at  cj°  ’ 

T(r,t)  =  temperature  (°F) 

A  =  area  (ft^) 

q  =  rate  of  heat  flow  (Btu/hr) 

O 

W  =  heat  generation  per  unit  volume  per  unit  time  (Btu/hr • ft0) 
mass  density  (=1208  lbm/ft^  for  tungsten) 
c  =  specific  heat  (=0.032  Btu/lbm°F  for  tungsten) 
k  =  thermal  conductivity  (=87  Btu/hr* ft*  F  for  tungsten) 
a  =  k/cp  =  thermal  diffusivity  (=2.430  ft^/hr  for  tungsten). 


(1  watt/cm^=3171  Btu/hr*ft2,  1  cal/ sec • cm • °C  =  242  Btu/hr • ft • °F) 
Equations  5.135  and  5.136  are  solved  for  the  four  cases  below. 


(1)  Constant  heat  flux  coming  into  a  semi -inf inite  body  (one  dimensional) 


(5.34)  . 

T  (x,t) 


n  Hl  *-4Se 


Ak 


Ire 


(5.137) 


(2)  Q  is  suddenly  generated  on  each  unit  area  of  a  plane  in  an  infinite 
body  at  x=0 ,  t=0  (one  dimensional). 


T«,t)  = 


Q 


2k 


.(X 

Jk t 


x 

e  ^ 


(5.138) 


(5.34)  L.R.  Ingersoll,  O.J.  Zobel  &  A.C.  Ingersoll:  "Heat  Conduction" 
McGraw-Hill  (1954) 
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(3) 


(4) 


Constant  heat  is  generated  inside  the  body  (uniform). 


Constant  heat  flux  is  emitted  at  r=0  (three  dimensional) 


(5.139) 

(5.35) . 

(5.140) 


For  our  electrode,  when  immersed  in  a  plasma,  ions  of  energy  eV^ 
hit  the  cathode  surface  and  transfer  an  energy  of  about  j.Vp.  In  the 

O 

typical  case  of  j^=10  A/cm  and  Vp=100  V,  the  temperature  rise  of  the 
tungsten  electrode  at  the  surface  is  from  Eq.  5.137,  only  19 °C  after 
1  msec,  which  is  far  too  low  for  thermionic  electrons  to  be  emitted. 

If  the  surface  of  the  electrode  is  covered  by  a  layer  of  W0^ ,  then 
P  =448  lbm/ft3,  c=0.073  Btu/lbm°F,  k=0.1-l  (5.36)  and  Eq.  5.137  yields 
a  temperature  rise  of  9-30°C. 


5.4.6  Local  Heating  due  to  Current  Focusing  on  the  Electrode 

When  the  sheath  breakdown  voltage  is  about  100  V,  the  current  into 
the  electrode  is  the  order  of  1  amperes.  If  this  current  is  focused 
into  a  small  area  of  the  electrode,  for  say  10  |Usec,  then  the  temperature 
rise  around  the  point  is,  from  Eq.  5.140, 

Ten  =  (i  —  erf  2o2 r) ,  (5.141) 

o 

where  r  is  in  cm  and  T  is  in  °C.  At  r=1000  A  from  the  current  injection 
point,  the  temperature  rises  to  T=5 . 37x10^ °C . 

(5.35)  M.  Jakob:  "Heat  Transfer  I"  John  Wiley  &  Sons,  N.Y.  pp339  (1949) 

(5.36)  "Gmelins  Handbuch  der  anorganischen  Chemie"  Verlag  Chemie ,  Bergstrasse. 
54  (1962) 
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At  the  cathode  spot  in  the  arc  discharge  the  current  density  was 
shown  to  be  between  a  few  thousand  and  more  than  a  million  A/ cmz  (5.37). 
The  temperature  obtained  above  is  sufficiently  high  to  emit  that  amount 
of  current  density  (Eq.  5.119). 


5.4.7  Local  Heating  of  Electrode  by  Ion  Bombardment 


The  rate  of  the  heat  flow  in  our  experimental  conditions  is  enormous 
(even  though  the  temperature  rise  is  very  small)  and  it  is  quite  possible 
that  thermal  equilibrium  is  not  attained  in  the  metal.  If  it  is  so,  then 
a  microscopic  treatment  should  be  used. 

The  average  depth,  L,  of  atoms  (of  atomic  number  Z^,  mass  m^  and 
initial  energy  E)  penetrating  into  metal  (of  atomic  number  >  mass  and 
density  n)  is  expressed  as  (5.38), 

1  3.06 


L  = 


l  +  i-at-  4icn5i;Z,'Zie*  » 

2/3,  „  2/3. %  __J  _  j.2  t  2  „  „„  ,„-8 


(5.142) 


where  a=0 . 8853ag/ (Z^  +Z2  )2  and  aQ=1i  /m^e  =0.523x10  cm  is  Bohr's 

critical  radius.  The  l-hiUj/Bm^  term  is  for  the  projected  path  length. 

o 

In  the  case  of  100  eV  argon  (atomic  diameter=3 .82  A)  penetrating  into 

0 

tungsten  (atomic  diameter=2 .82  A), 

T  =  0.478  A. 

The  relaxation  time  for  phonon-phonon  collision  is  (5.39), 

*Cpp=10 sec 

and  that  of  electrons  for  electron-phonon  collisions  is, 

-13 

7**  =10  sec 

ep 

ft  =10’U  sec. 
tpe 


(5.37)  J.D.  Cobine  &  E.E.  Burger:  J.  Appl.  Phys .  26.  895  (1955) 
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If  the  heat  flux  input  is,  as  before,  100  W/cm^  which  is  confined  into 

o 

a  0.478  A  thick  layer  along  the  metal  surface,  then  the  temperature  rise 
of  the  layer  in  10”^  sec  is,  from  Eq .  5.139,  0.9°C.  After  this  time  the 
energy  of  the  layer  begins  to  be  diffused  into  the  metal  by  phonons,  and 
T  is  specified  by  Eq .  5.137.  Again  this  temperature  is  very  low. 

If  an  incident  atom  hits  only  one  metal  atom,  then  the  atom  is 

TO 

heated  to  the  temperature  of  eVp/k  for  10”  J  seconds,  after  which  it 
begins  to  cool  off.  If  ion  current  with  ion  energy  of  100  eV  flows  into 
an  electrode,  then  each  ion  heats  electrode  atoms  to  1. 16x10^ °K.  Thermionic 
electron  current  from  one  atom  of  Ttx(l  .41x10  8)^x9. 7x10"^  (Amps)  flows  for 
10”  sec  (1.41x10  is  the  radius  of  tungsten  atom),  so  that  the  ratio  of 
secondary  electron  current  to  primary  ion  current  becomes  3.8x10^. 

Although  this  current  increase  rate  looks  quite  impressive  for 
the  explanation  of  breakdown  phenomena,  we  have  to  give  up  this  idea  because  of 
the  fact  that  the  current  ratio  thus  obtained  is  independent  of  the  plasma 
density  (i.e.,  ion  current  into  the  electrode).  From  the  experimental 
results,  the  breakdown  voltage  is  a  function  of  plasma  density.  The  fault 
of  this  theory  might  be  the  use  of  Richardson's  equation,  which  is  not 
proved  to  be  correct  in  the  thermally  non-equilibrium  state. 


(5.38)  G.  Carter  &  J.S.  Colligon:  "Ion  Bombardment  of  Solids"  Heinemann 
Educational  Books  Ltd.,  London  pp97  (1968) 

(5.39)  C.  Kittel:  "Introduction  to  Solid  State  Physics"  John  Wiley  &  Sons 
Inc.  N.Y.  (1967) 
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CHAPTER  VI  EXPERIMENTS  ON  THE  BREAKDOWN  BETWEEN  CYLINDRICAL  ELECTRODES 

IN  A  PLASMA 


6.1  BREAKDOWN  VOLTAGE  IN  A  NEUTRAL  GAS 

The  breakdown  voltage  (V^)  between  two  cylindrical  tungsten  electrodes 
in  an  unionized  argon  gas  was  measured.  The  radius  (rp)  and  the  length  of 
both  electrodes  were  0.05  mm  and  35  mm,  respectively.  The  interelectrode 
distance  (d)  and  the  ambient  pressure  (p)  were  varied.  In  Fig.  6.1  is 
shown  the  measured  breakdown  voltages.  Theoretical  V^  obtained  from  Eq.5.24 
and  A=13.6  ( cm “ ^mmHg " ^ )  ,  B=235  (V/cm-mmHg)  (5.2)  andY  =  0.1  are  shown  in 
Fig.  6.2. 

Curves  of  Fig.  6.1  are  replotted  in  terms  of  V^/lnCd/rp)  in  Fig.  6.3, 
together  with  the  theoretical  values  obtained  from  Eq .  5.24.  As  expected 
from  Eq .  5.28,  all  the  curves  tend  to  coincide  with  one  line  except  for 
small  p  and  large  d.  When  p  is  small  the  last  term  of  Eq .  5.24  becomes 
significant  and  Eq .  5.28  is  no  longer  valid;  when  d  is  large  the  streamer 
mechanism  takes  place  and  changes  the  picture  of  the  breakdown  phenomena 
greatly.  Also  in  Fig.  6.3  the  results  of  Meek  and  Craggs  (6.1)  are  plotted. 
In  their  case  a  cylindrical  wire  and  concentric  cylindrical  electrodes 
were  used;  the  wire  electrode  was  negative  (electrode  material  unknown). 

The  constancy  of  V  /2p  (d-r  )  of  Eq .  5.26  from  our  data  is  examined 
J  m  rm N  p  n 

in  the  Table  6.1.  The  experimentally  obtained  B  is  approximately  30% 


(6.1)  J.M.  Meek  &  J.D.  Craggs:  "Electrical  Breakdown  of  Gases"  Oxford, 
Clarendon  Press  pp60  (1953) 
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smaller  than  the  value  given  by  Meek  and  Craggs  (=235)  and  77,  smaller 
than  the  value  given  by  von  Engel  (5.3)  (=180)  except  for  large  d.  When 


d  is  large  V  and  hence  B  is  reduced  by  a  streamer  mechanism.  In  addition, 
b 

the  35  mm  long  wire  electrode  has  to  be  approximated  by  a  point  electrode 

rather  than  a  cylindrical  electrode.  The  theoretical  absolute  minimum 

voltage,  V  can  be  calculated  from  Eq .  5.27  to  be  225  volts  for  Meek's 
°  ’  mm  n 

values  of  A  and  B,  and  195  volts  for  von  Engel's  values.  The  minimum 
sparking  voltage  for  parallel  plate  electrodes  was  found  to  be  195-275 


volts  (6.2). 


It  should  be  noted  that  contrary  to  the  remarks  of  von  Engel  (6.3) 
that  in  the  case  of  breakdown  between  two  concentric  cylindrical  electrodes 
the  equivalent  formula  for  Eq .  5.24  does  not  agree  with  observations,  the 
experimental  results  in  our  case  agreed  reasonably  well  with  the  theory 
developed  here. 


6.2  SHORT  GAP  BREAKDOWN  EXPERIMENTS 


The  hair  pin  type  electrodes  used  in  the  plasma  breakdown  experiments 
were  mounted  face  to  face  in  a  specially  made  electrode  holder  with  which 
the  distance  between  the  tips  of  the  electrodes  can  be  adjusted  within 
an  accuracy  of  lxlO-^  inches.  This  device  was  put  into  the  vacuum  jar 
(pressure=0 . 05  mmHg)  and  the  breakdown  voltage  between  these  two  electrodes 
was  measured. 


(6.2)  L.B.  Loeb:  "Fundamental  Processes  of  Electrical  Discharge  in  Gases" 
John  Wiley  &  Sons  Inc.  N.Y.  pp414  (1939) 

(6.3)  A.  von  Engel:  "Ionized  Gases"  Oxford,  Clarendon  Press  ppl91  (1965) 
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d(cm) 

S’ 

/'■"N 

33 

0Q 

V 

v  (v) 

m N 

V  /2p  (d-r  ) 

m  Mn v  p 

B  from  (6.1) 

0.1 

5 

255 

260 

235 

0.25 

2.5 

265 

210 

If 

0.5 

2 

340 

170 

VI 

1.1 

1 

350 

160 

It 

2.1 

0.5 

380 

180 

ii 

5 

0.3 

410 

140 

If 

10 

0.18 

450 

130 

11 

20 

0.14 

600 

110 

11 

50 

0.1 

780 

80 

II 

TABLE  6.1  The  minimum  breakdown  voltage,  Vm,  and  which 

gives  V  at  various  values  of  d.  In  the  table  the  experimentally 
m 

obtained  B  (=Vm/2pm(d-rp) )  and  the  value  of  B  taken  from  Meek  (6.4) 


are  also  compared. 
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The  electrodes  used  were  bombarded  by  argon  ions  in  the  shock  tube 
for  various  times  beforehand.  The  measured  breakdown  voltage  as  a 

_3 

function  of  electrode  separation  (d)  yielded  straight  lines  (d=0-3xl0  cm 
V^=0-800  V)  which  showed  that  the  electric  field  at  the  electrode 
surface,  ,  were  constant.  The  value  of  Ep  changed  with  the  degree  of 
bombardment  and  was  ,  from  Eq .  5.19, 

Ep  =  2  -  4x1  O'*  V/cm. 

The  same  kind  of  short  gap  vacuum  breakdown  has  been  investigated  by 
many  people,  especially  Kisliuk  and  others  (7.4),  and  the  breakdown  was 
found  to  be  triggered  by  electron  field  emission  from  the  cathode  surface. 
This  theory  applies  to  our  case.  For  our  hair  pin  type  electrodes, 
E=2-4xl0^  V/cm  is  necessary  to  extract  electrons  from  the  cathode.  The 
field  required  for  the  field  emission  is,  from  Sec.  5.3.1,  4x10^  V/cm, 
and  the  discrepancy  is  attributed  to  the  field  intensification  by  the 
surface  irregularities.  By  using  the  notation  of  Sec.  5.3.2  p  for  our 
case  is  calculated  as  , 

p  =  100  -  200. 

$  is  large  due  to  ion  bombardment. 


6.3  THE  EFFECT  OF  ELECTRODE  CONTAMINATION  ON  THE  BREAKDOWN  VOLTAGE 


Two  cylindrical  electrodes  were  immersed  in  the  plasma  and  the 
electrical  breakdown  voltage  between  these  two  was  measured.  The  plasma 
was  produced  by  the  electromagnetically  driven  shock  tube  mentioned  in  the 
previous  sections. 

After  some  trials  it  was  found  that  in  addition  to  very  bad  reproduci 
bility  of  the  breakdown  voltage,  there  were  no  correlations  between 
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breakdown  voltage  and  plasma  density,  neutral  particle  density,  electrode 
separation,  etc.  However,  by  integrating  the  electrode  current  with 
respect  to  time  until  breakdown  occurs  (t.)  it  was  found  that 

A 

\  Idt  =2-4  (  Mcoul) 

J0 

for  most  of  the  cases  at  various  experimental  conditions.  Detailed 
examples  are  listed  in  Table  6. 2 -A.  The  reason  of  the  constancy  of  the 
integral  can  be  considered  as  follows.  Due  to  the  intense  arc  discharge 
in  the  shock  tube  driver,  small  portions  of  electrodes  and  insulator 
materials  of  the  driver  are  melted,  pushed  downstream,  and  some  of  them 
stick  to  the  surface  of  the  electrode,  forming  an  insulating  surface  layer. 
In  the  presence  of  oxygen,  the  oxidation  of  the  electrode  surface  also 
makes  an  insulating  layer.  Now,  if  a  potential  difference,  ,  is  applied 
on  two  electrodes  in  a  plasma,  these  electrodes  act  as  a  floating  double 
electrostatic  probes  and  begin  to  collect  ions  and  electrons.  However, 
because  of  the  insulating  surface  layer,  the  electric  charge  starts 

accumulating  on  the  surface  of  the  layer,  making  an  electric  field  of  the 

/-"tb 

order  of  Idt)/£QS  inside  the  layer,  where  S  is  the  electrode  surface 
area.  When  this  electric  field  becomes  large  enough  (  -  3x10^  V/cm)  an 
electron  avalanche  starts  inside  the  layer,  which  produces  enough  electrons 
to  trigger  the  breakdown  between  electrodes.  The  same  effect  was  reported 
by  Hancox  (6.4),  who  made  an  artificial  surface  layer  and  found  that 

A 

\  Idt  was  constant. 

Jo 

To  eliminate  this  effect,  the  electrodes  were  heated  to  approximately 
2300°C  for  one  minute  before  each  measurement  (Sec.  5.3.7).  With  cleaned 

a 

electrodes  j^Idt  was  measured  again  with  the  results  in  Table  6.2-B.  In 
this  case  the  integral  did  not  take  constant  values.  For  contaminated 


(6.4)  R.  Hancox:  Brit.  J.  Appl.  Phys .  11_  468  (1960) 
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Date 

z  (cm) 

VC(KV) 

P^(mmHg) 

> 

'w' 

> 

\  Idt  (A-fisec) 
J0 

Mar.  18 

68 

2.5 

0.1 

110 

2 

Jun.  25 

ii 

ti 

0.1 

130 

2  .5 

II 

II 

ii 

0.05 

150 

3 

II 

II 

it 

0.5 

125 

4 

July  22 

it 

It 

1 

350 

2 

it 

II 

ti 

10 

350 

2.5 

TABLE  6. 2 -A 

\  Idt 
Jo 

with  contaminated  electrodes. 

DATE 

z  (cm) 

VC(KV) 

p^mmHg) 

Vb(V) 

/■’tb 

l  Idt  (A-^sec) 

J0 

July  22 

68 

2.5 

0.1 

600 

0.2 

it 

II 

II 

0.1 

200 

20* 

ti 

ii 

It 

1 

600 

0.7 

II 

ti 

ii 

10 

600 

1.3 

ii 

111 

II 

0.1 

600 

0.1 

it 

ii 

it 

0.1 

200 

20* 

II 

ii 

ii 

1 

600 

0.56 

II 

ii 

it 

5 

600 

0.6 

TABLE  6.2-B 


with  clean  electrodes. 


*  indicates  that  no  breakdown  was  observed. 
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electrodes  to  breakdown,  the  integral,  not  ,  is  important  and  can  be 
much  smaller  than  the  true  breakdown  voltage. 

Besides  leading  to  an  interelectrode  breakdown,  electric  charge 
accumulation  often  causes  small,  successive  electron  avalanches  in  the 
contamination  layer  of  the  electrode  surface.  This  effect  can  be  seen 
in  Figs.  6.4  and  6.5.  In  the  following  sections  electrodes  were  always 
cleaned  before  taking  measurements,  unless  otherwise  mentioned. 


6.4  PRE-BREAKDOWN  CURRENT  CHARACTERISTICS 


If  the  electrode  potential,  V,  is  increased  without  changing  the 
plasma  properties,  the  electrode  current,  I,  experiences  three  stages, 
i  .  e . : 

Stage  (1)  When  V  is  far  below  the  breakdown  voltage,  ,  the  electrode 
works  as  a  Langmuir  probe,  and  I  can  be  expressed  by  the  equations  in 
Sec .  2.7. 

Stage  (2)  When  V  is  increased,  electron  avalanches  begin  to  take  place 
in  the  cathode  sheath,  which  cause  an  increase  in  I. 

Stage  (3)  At  V^V^,  I  is  dominated  mainly  by  the  external  circuit. 

In  actual  cases,  it  has  been  known  from  the  experiments  that  there 
are  two  types  of  I  vs  V  characteristics,  as  seen  in  Figs.  6. 6 -A  and  B. 

Three  stages  are  clearly  shown  in  Fig.  6.6-A,  whereas  in  Fig.  6.6-B,  stage 
(2)  is  unnoticeable  and  stage  (3)  follows  immediately  after  stage  (1). 

In  Fig.  6.7  are  shown  the  C.R.O.  traces  of  Ip  with  various  values  of 
V.  Constant  V  was  applied  from  t=0  to  infinity.  The  plasma  surrounding 
the  electrodes  was  a  weak  precursor  plasma  (6.5).  It  reaches  the  electrodes 
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FIG.  6.4 


FIG.  6.5 
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FIG.  6.4  Oscilloscope  traces  of  the  currents  flowing  into  a 
contaminated  electrode.  p^=0.1  mmHg ,  V^=2.5  KV  and  the  negative 
and  positive  electrode  positions  =  70  cm  and  140  cm  from  the  driver, 
respectively.  Vertical  and  horizontal  deflections  =  0.2  A/div.  and 
20  psec/div.  Time  runs  from  left  to  right.  Upper  beam:  V^=80  V, 
middle  beam:  7^=70  V  and  lower  beam:  V(j=50  V. 


FIG.  6.5  The  change  of  the  electrode  current  with  different 
degree  of  contamination.  V^=400  V.  Vertical  and  horizontal  deflections 
=  10  mA/div.  and  20  psec/div.  Upper  beam:  clean  electrodes,  middle 
beam:  lightly  contaminated  electrodes  (after  2  shots  without  cleaning) 
and  lower  beam:  heavily  contaminated  electrodes  (after  10  shots 
without  cleaning). 
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FIG.  6.6 


Pre-breakdown  V-I  characteristics. 


A.  Stages  1,  2  and  3  can  be 

clearly  distinguished. 

Stage  1:  V-I  curve  follows 
P 

the  probe  theory. 

Stage  2:  Electron  avalanches 
increase  the  electrode  current. 


Stage  3:  Ip  is  mainly  dominated 
by  external  circuit. 


B.  Stage  2  cannot  be  observed, 
or  hidden  by  an  experimental 
error . 
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within  1  |Asec  and  changes  its  density  as  shown  approximately  in  the 
V=200  volts  wave  form,  where  the  electrodes  were  working  as  floating  double 
probes.  The  current  gradually  increases  with  increasing  V  and  when  V=610  V, 
a  large  secondary  process  occurs  with  a  time  lag  of  24  psec  (6.6)  from  the 
time  of  the  first  peak  current  (  at  6  psec)  .  At  this  time  (V=610  V,  t=30^s), 
the  relation  of  holds,  since  I/Iq  approaches  infinity,  where  Iq 

is  the  current  when  electron  avalanche  effect  can  be  negligible.  With 
higher  V,  I  goes  to  infinity  with  a  certain  time  lag,  depending  on  the 
magnitude  of  overvoltage  and  the  plasma  density. 

An  approximate  theory  expressing  the  relation  of  I  and  V  will  be 
presented  when  the  ionizing  effect  in  the  sheath  by  secondary  electrons 


is  considerable.  Since  V  is  applied  across  the  cathode  sheath,  I  can  be 
considered  as  a  current  flowing  between  concentric  cylindrical  electrodes 
with  inner  and  outer  electrode  radii  of  r^  (electrode  radius)  and  rs  (sheath 
radius),  respectively.  I  can  be  expressed  as, 

(  o<dr 

He  eV 


l  -'r(e^0<dr-i) 


where 


rCs  A\J  f  -  gprP-^-rP 


i(e 


Bp  •"s  (in  Tp 


V 


-C 


) 


(6.1) 


(6.2) 


(6.5)  The  precursor  plasma  in  our  shock  tube  is  produced  mainly  by  the 

photoionization  by  the  intense  radiation  of  the  arc  discharge  in  the 
driver.  This  precursor  plasma  has  a  high  velocity  (10-100  times  the 
velocity  of  the  shock  produced  plasma),  a  low  density  (lO^-lO^  cm-3) 
and  a  slow  density  decay  along  the  axis  of  the  shock  tube  (  -  10  times 
the  decay  time  of  the  shock  produced  plasma) . 


(6.6)  When  the  breakdown  voltage  was  measured  with  shock  produced  plasma, 
the  breakdown  occured  at  or  before  the  density  peak  arrives  at  the 
electrode.  The  time  lag  mentioned  here  was  observed  only  in  the  case 
of  precursor  plasma.  This  may  be  due  to  the  incomplete  surface  layer 

cleaning  by  ion  sputtering. 
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and  I  is  the  initially  emitted  electron  current  from  the  cathode.  Here 


j  =Ae  e 


£X 


and 


EV)  = 


(6.3) 


were  employed.  Also  it  was  assumed  that  the  field  perturbation  by  space 

charge  was  negligible.  Equation  6.1  was  computed  with  A=13.6  cm"^mmHg  ^ , 

B=235  V/cm*mmHg  (5.2),  r  =0.05  mm,  r  =0.01-100  cm,  V=0-2000  V  and  p=0 . 01 - 

P  s 

100  mmHg  and  the  results  presented  in  Figs.  6.8  and  6.9.  From  Fig.  6.8 
we  can  see  that  I/Ie  becomes  a  maximum  at  a  certain  combination  of  p  and 
rg  given  by 


— ; - r  =  B  =  constant. 

P(rs-rP)  rs 


(6.4) 

does  not  reach  the  value  of 


At  small  rg  in  Fig.  6.9,  exp(j°<clr) 

1+1  Ay  (i.e.,  breakdown)  even  with  V=Q0  .  However,  at  large  V,  the  electric 

field  at  the  electrode  surface  becomes  very  intense  and  electron  field 

emission  begins  to  take  place  which  leads  to  the  voltage  collapse  in  the 

sheath.  In  this  case  of  field  emission,  stage  (2)  in  Fig.  6.6-A  will 

never  be  observed,  since  the  field  emission  current  is  a  very  strong 

function  of  electric  field  and  so  V-I  curve  takes  the  form  of  Fig.  6.6-B. 

In  other  words,  V-I  curve  is  of  type  Fig.  6.6-A  when  p  is  large  and/or 

r  is  large  (n  is  small)  and  of  type  Fig.  6.6-B  when  p  is  small  and/or 
s  ^ 

r  is  small  (n^  is  large). 

In  our  experiments  V  was  kept  constant,  while  ng  was  increased  with 
respect  to  time,  i.e.,  rg  was  decreased.  In  this  case,  too,  I-rg  curves 
experience  qualitatively  the  same  curves  as  shown  in  Figs.  6.6-A  and  B. 
These  two  types  of  the  current  change,  observed  in  the  experiments  at 
various  values  of  p  and  ng ,  were  in  good  correspondense  with  the  theory. 
Fig.  6.7  belongs  to  the  Fig.  6.6-A  type.  In  most  of  the  cases,  however, 
the  type  of  Fig.  6.6-B  was  observed. 
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6.5  MEASUREMENTS  OF  BREAKDOWN  VOLTAGE  vs  ELECTRODE  SEPARATION 

According  to  the  floating  probe  theory  mentioned  before,  it  is 
shown  that 

vl=vd+vf-vc  (6.5) 

at  sufficiently  large  V,  (Sec.  2.7).  If  V  is  increased,  then  only  V-, 

u  d  1 

is  increased  while  V2  and  stay  constant.  Therefore  the  breakdown 

should  always  occur  in  the  sheath  of  the  negative  electrode  no  matter 

what  the  positive  electrode  condition  is. 

To  examine  this  statement  the  breakdown  voltage  was  measured  with 

cathode  set  stationary  while  the  anode  was  moved  along  the  shock  tube 

axis.  Under  these  conditions  the  plasma  density  and  the  temperature  at 

the  anode  decreased  as  it  was  moved  away  from  the  shock  driver. 

The  measurements  were  carried  out  with  the  initial  pressure  of  the 

shock  tube  of  0.15  mmHg  in  argon  and  with  the  cathode  kept  at  a  distance 

12  -3 

of  100  cm  from  the  driver,  where  the  plasma  density  was  2.5x10  cm  .  By 
increasing  while  observing  the  change  in  electrode  current,  the 
breakdown  voltage  was  measured  as  a  function  of  electrode  separation,  with 
the  results  shown  in  Fig.  6.10.  From  the  figure  it  is  evident  that  the 
breakdown  voltage  was  always  constant  and  independent  of  the  electrode 
separation  and  anode  condition.  In  the  figure  the  breakdown  voltage  in 
the  absense  of  plasma  is  also  plotted.  On  the  other  hand,  when  the  cathode 
was  moved,  the  breakdown  voltage  showed  a  change  depending  on  the  plasma 
density,  etc.  at  the  cathode.  These  facts  clearly  show  that  the  breakdown 
took  place  at  the  cathode. 

When  the  electrical  conductivity  of  the  plasma  becomes  low,  there 
exists  a  voltage  drop  in  the  plasma  and  V-^  is  decreased  by  some  amount. 
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Let  us  examine  this  effect  by  taking  the  example  of  the  severest  case,  as, 

Oi  (degree  of  ionization)  =  10"  , 

T  =  1000  K, 

Qe  =  3xl0-2 1  m2  (Fig.  2.1) 

and  n  =  10^  cm  ^ , 

e  ’ 

then  from  Eq.  2.5  the  plasma  conductivity  becomes, 

<T*  =  40  mhos/m. 

If  the  electrode  distance  is  10  cm  and  the  effective  cross  section  of 

O 

the  plasma  cylinder  is  10  cm^ ,  then  the  resistance  of  the  plasma  between 
two  electrodes  is , 

R  =  2.5 

With  the  electrode  current  of  0.1  amperes  the  voltage  drop  is  only  0.25  volts, 
which  can  be  neglected. 


6.6  MEASUREMENTS  OF  BREAKDOWN  VOLTAGE  vs  PLASMA  DENSITY 


Breakdown  voltages  between  two  electrodes  were  measured  at  various 
plasma  and  neutral  densities.  Electrode  currents  at  the  same  ambient 
conditions  as  well  as  the  plasma  velocities  were  also  measured  to 
facilitate  the  calculation  of  the  plasma  and  neutral  densities  (Sec.  7.2). 
The  electrode  current  variation  with  respect  to  applied  voltage  at  various 
shock  driver  conditions  is  shown  in  Fig.  6.11.  The  breakdown  voltage 
against  plasma  density  at  initial  pressures  of  0.1,  1  and  10  mmHg  in  argon 
is  also  shown  in  Fig.  6.12  (see  Table  7.1). 

It  should  be  noted  that  in  Figs.  6.11  and  6.12,  decreases  to  a 
value  much  lower  than  the  Paschen  minimum  voltage  in  Sec.  6.1  and  Eq .  5.27 
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of  V  =200  volts.  This  shows  that  the  breakdown  phenomena  is  not  as 

mm  1 


mm 


simple  as  the  breakdown  in  a  neutral  gas.  The  surrounding  plasma  is 
causing  a  large  influence  on  the  breakdown  voltage. 

The  shot-to-shot  method  was  used  for  the  voltage  and  current 
measurements.  The  reproducibility  was  sometimes  not  very  good.  After 
some  trials,  however,  it  was  found  that  if  the  data  was  interpreted  as 
a  function  of  time  after  firing  the  driver,  instead  of  as  a  function  of 
axial  position  of  the  shock  tube,  the  reproducibility  was  fair  to  good. 

This  indicated  that  d/dt  of  the  plasma  properties  (e.g.  recombination, 
ambipolar  diffusion)  was  quite  reproducible  while  v(d/dt)  of  the  plasma 
properties  (e.g.  heating  the  oncoming  neutrals  by  the  shock,  frictional 
force  between  the  wall  and  plasmoid)  was  not.  Therefore  the  data  was  taken 
as  a  function  of  time,  and  then  normalized  into  a  function  of  axial 
distance  whenever  necessary. 

6.7  EFFECT  OF  THE  ELECTRODE  EDGE  ON  THE  BREAKDOWN  EXPERIMENTS 

With  our  electrode  configuration  there  was  a  limit  to  measurements 
of  the  breakdown  voltage.  Since  the  electrodes  used  here  traversed  the 
shock  tube  as  seen  in  Fig.  3.3,  the  tip  of  the  electrode  obviously 
measured  the  highest  plasma  density,  while  the  bottom  part  of  the  electrode 
measured  somewhat  lower  plasma  density.  This  density  variation  can  be 
approximated  by  (Eq.  2.30), 


(6.6) 
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order  of  the  first  kind. 

In  Fig.  6.13  are  shown  the  V,  -n  curves  which  would  have  been 

be 

measured  at  various  positions  of  the  cathode  if  we  could  have  separated 

the  cathode  into  small  segments.  For  instance  at  R=0  (cathode  tip), 

follows  the  curve  F-E  — C  as  ng  is  increased  from  zero,  coinciding  with 

the  Paschen  curve  of  Eq.  5.24.  With  higher  ng,  electron  emission  from 

the  cathode  becomes  dominant  over  electron  creation  by  collision  and 

becomes  smaller  again,  as  C-B  — A.  However,  at  the  density  of  n£ ,  for 

instance,  Vu  at  the  position  of  R=R  on  the  cathode  is  only  V  .  and  so 
’  b  r  c  mm 

all  we  can  detect  as  a  breakdown  voltage  is  Vm^n,  not  Vc .  Consequently  if 
the  plasma  density  is  increased  gradually,  measured  changes  along  the 
thick  line  as  F-  E-D-B-A.  The  E-C-B  curve  can  never  be  observed. 
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FIG.  6.13  The  breakdown  voltage  at  various  positions  of  the 

cathode  as  a  function  of  plasma  density.  R=0  at  the  cathode  tip 

and  R=Rq  at  the  cathode  bottom  (Rq  is  the  radius  of  the  shock 

tube).  Measured  V,  changes  as  A-B-D-E-F  as  n  is  decreased. 

b  e 
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CHAPTER  VII  DISCUSSIONS  ON  THE  BREAKDOWN  THEORY  AND  EXPERIMENTS 


7.1  PROPOSED  MODEL  OF  THE  BREAKDOWN  MECHANISM 


By  having  checked  all  the  theoretical  possibilities  which  the  author 
could  have  thought  of  on  the  breakdown  mechanism  and  by  having  done  many 
experimental  investigations,  the  most  likely  mechanism  of  the  breakdown 
between  the  electrodes  immersed  in  a  plasma  is  now  proposed. 

THE  MODEL  OF  THE  BREAKDOWN  MECHANISM:  There  are  two  independent  kinds 

of  mechanisms  for  the  breakdown  between  two  electrodes  in  a  plasma;  one  is 

the  "Cold  cathode  arc  discharge"  and  the  other  is  the  "Spark  discharge". 

(1)  Cold  cathode  arc  discharge  takes  place  at  relatively  high  plasma 

densities.  A  very  thin  sheath  is  formed  around  the  cathode  due  to  the 
high  plasma  density  and  low  ionic  mobility.  Also  the  current  enhancement 
by  the  ionizing  effect  of  electrons  inside  the  sheath  causes  the  reduction 
of  the  sheath  thickness.  Because  of  the  thin  sheath  thus  created  the 
electric  field  at  the  cathode  surface  becomes  strong  enough  to  cause 
electron  ejection  by  field  emission,  which  triggers  the  discharge.  After 
the  triggering  the  current  is  concentrated  into  a  small  area  of  the  cathode 
and,  due  to  its  local  heating,  thermionic  electrons  dominate  the  field 
electrons  as  a  main  current  source  for  the  full  path  breakdown. 

(2)  Spark  discharge  occurs  at  lower  plasma  densities.  The  secondary 
electrons  emitted  from  the  cathode  ionize  neutral  particles  in  the  cathode 
sheath  and  this  ionization  causes  the  instability  of  the  electrode  system. 
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When  the  anode  is  contained  in  the  sheath,  the  breakdown  mechanism  is 
similar  to  that  for  zero  plasma  density  (which  is  the  extreme  case  in 
which  the  sheath  is  infinitely  large) ,  except  that  due  to  the  space 
charge  accumulation  first  ionization  coefficient  becomes  larger  and  the 
breakdown  voltage  decreases. 

In  the  following  part  of  this  section  the  equations  which  should  be 
used  in  executing  the  proposed  models  will  be  listed.  For  the  investigation 
of  the  cold  cathode  arc  mechanism,  the  following  equations  and  figures 
should  be  employed: 

(1)  Applied  voltage  across  the  cathode  sheath 

Vf:  Eqs .  5. 1-5. 7 

Vp:  Eqs.  2.83-2.106,  Fig.  2.5 

(2)  Ionic  mobility  near  the  cathode 


Y±:  Eq.  5.85 

n:  Fig.  4.2,  Eqs.  5.88  and  5.110 
T:  Fig.  4.4,  Eq .  5.90 


(3)  Cathode  sheath  and  electrical  boundary  layer  thickness 

r  :  Eqs.  5.52,  5.53,  5.56,  5.57,  5.63,  5.64,  5.65  and  5.66 
s 

ra:  Eq.  5.41 

(4)  Effective  current  in  the  sheath 

^(charge  accumulation):  Sec.  5.3.5 
G:  Eqs.  5.46,  5.55,  5.60,  5.61  and  5.62 
Qi:  Fig.  5.1 

(5)  Electric  field  at  the  cathode  surface 

Ep:  Eqs.  5.52,  5.56,  5.63  and  5.66 

|3(field  enhancement  factor  by  surface  irregularities):  Secs.  5.3.2 
and  6.2,  Eq .  5.114 

£ (field  enhancement  factor  by  stagnation):  Sec.  5.14,  Eqs.  5.36, 


5.37  and  5.40. 
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(6)  Field  emission  current 

j:  Eq.  5.111 

(7)  Cathode  spot  temperature  by  local  heating 

T:  Eqs.5.140  and  5.141 

(8)  Thermionic  electron  emission  current 

j:  Eq.  5.119 

(9)  Self  sustaining  arc  discharge 

Eqs.  5.120-5.134. 

As  there  is  no  collision  between  ions  and  electrons  in  the  sheath 
(Sec.  5.1.8),  the  ion  density  in  the  sheath  is  not  decreased  due  to  the 
boundary  layer  effect  although  its  temperature  is  decreased.  This  effect, 
therefore,  need  not  be  considered. 

For  the  spark  discharge,  G  and  rg  are  given  by  Eqs.  5.55-5.57  and 
Fig.  5.2.  In  the  figure,  if  Vp/rpEm  becomes  large  and  d/rp  becomes  close 
enough  to  its  final  value,  then  even  a  very  small  perturbation  to  the 
sheath  thickness  can  lead  the  system  to  the  breakdown  (i.e.,  G=l);  therefore 
this  is  the  condition  for  breakdown.  If  rg  becomes  larger  than  the  inter¬ 
electrode  separation,  then  the  modified  form  of  Eqs.  5.22-5.24  should  be 
used.  (There  is  a  change  in  potential  distribution  due  to  the  formation 
of  sheath.)  In  our  experimental  conditions,  however,  rg  is  much  smaller 
than  the  spark  discharge  criterion  and  the  discussion  on  these  mechanisms 
will  not  be  conducted  any  further. 


7.2  COMPARISON  OF  THE  BREAKDOWN  MODEL  WITH  THE  EXPERIMENTAL  RESULTS 


In  this  section  the  experimentally  obtained  values  will  be  examined 
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according  to  the  model  theory  of  the  previous  section.  If  the  sheath 

thickness  is  small  for  the  spark  discharge  criterion,  then  the  electric 

field  at  the  cathode  surface  should  be  strong  enough  to  cause  field  emission, 

or  in  other  words  the  applied  field  at  the  cathode  surface  should  be 

2-4x10^  V/cm  (Secs.  5.3.1,  5.3.2  and  6.2).  Therefore  the  intensity  of  the 

electric  field  at  the  cathode  surface  will  also  be  calculated. 

In  Table  7.1  the  experimentally  obtained  values  of  M,  and  many 

other  calculated  values  are  listed.  (ne  vs  curve  is  shown  in  Fig.  6.12.) 

The  relation  of  •,3Cd«A  T  in  the  table  implies  that,  in  the  sheath,  the 

iw  ew 

ions  are  collision  dominated  while  electrons  are  not.  For  the  ion  saturation 
current  to  the  electrode,  collision  dominated  theory,  especially  Eq .  2.74 
should  be  used  since  RgOC^X^^l  (not  listed).  For  the  floating  potential 
of  electrodes,  the  argument  of  Sec.  5.1.1  is  correct.  is  an 

important  relation,  which  indicates  that  Navier-Stokes  equation  is 
applicable  to  our  case  and  the  boundary  layer  is  formed  around  the  electrode. 
In  all  the  cases  G„«l  and  also  1/1  (not  listed)  show  that  the 

a 

ionization  effect  in  the  ambipolar  diffusion  region  as  well  as  that  outside 
the  ambipolar  diffusion  and  sheath  regions  is  negligible. 

Ep2  and  ^2^rp  are  0frtai-necl  with  the  assumptions  that  the  boundary 
layer  effect  is  totally  negligible,  i.e.,  <^E<^'\W  anc^  ^  d  •  On  the  other 
hand  EpW  and  dw/ are  obtained  with  the  assumptions  that  the  boundary 
layer  effect  is  dominating,  i.e.,  and  d.  From  the  table  it 

can  be  found  that  and  d ,  which  means  that  EpW  and  dw  are 

more  correct  than  Ep2  and  d2 .  Especially  for  the  calculation  of  Ep ,  the 
assumption  of  (S'g'^d  is  not  very  significant  to  the  value  of  Ep ,  since  the 
Poisson's  equation  was  integrated  in  the  cathode  region  and  also  the 
charge  density  is  the  largest  at  the  cathode  surface  (Eq.  5.48).  However, 
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TABLE  7.1  Calculation  of  the  electric  field  at  the  cathode  surface 

Experimental  and  calculated  results  according  to  Sec.  7.1  are  listed  here. 

The  notations  used  are: 

=  initial  pressure  in  the  shock  tube 
z  =  axial  distance  from  the  driver  mouth 
=  capacitor  bank  charging  voltage 
M  =  ^/c-L  =  normalized  plasma  velocity 
tarr=  arrival  time  of  plasmoid  at  the  electrodes 

=  electrode  current  just  before  the  breakdown  (7.1) 

=  breakdown  voltage 

n2/no  =  the  ratio  of  the  neutral  particle  density  in  the  incident  shock 
to  the  neutral  density  at  1  mmHg  and  0  C  (nQ=3 .22x10^  m-^). 
nw/ng  =  the  ratio  of  the  neutral  particle  density  near  the  cathode 
to  the  neutral  density  at  1  mmHg  and  0  C 
=  ionic  mean  free  path  near  the  cathode 
\ew  =  electron  mean  free  path  near  the  cathode 
Rg  =  electrical  Reynold's  number 

/r  =  the  ratio  of  the  thickness  of  electrical  boundary  layer  to 
the  electrode  radius 

G_  =  ionizing  effect  in  the  ambipolar  diffusion  (Eq.  5.62) 
ng  =  electron  density  in  the  incident  shock 
T0  =  electron  temperature  in  the  incident  shock 

=  electrical  conductivity  in  the  incident  shock 


E- 


176 


Ep2  =  electric  field  at  the  cathode  surface  when  rs^>ra  (Eq.  5.56) 
d£/rp  =  ratio  of  the  sheath  thickness  to  the  electrode  radius 
when  rg^>ra  (Eqs.  5.55-5.57) 

EpW  =  electric  field  at  the  cathode  surface  when  ra^>rg  (Eq.  5.56 

with  G  and  suffix  2  changed  to  Gg+Ga  and  suffix  w,  respectively) 
dw/rp  =  ratio  of  the  sheath  thickness  to  the  electrode  radius 

when  r J$>r  (Eqs.  5.55-5.57  with  G  and  suffix  2  changed  to 
G  +G  and  suffix  w,  respectively), 

o  cl 


(7.1)  "Just  before  the  breakdown"  means  that  the  current  at  the  end  of 
stage  (1)  of  Figs.  6.6-A  and  B.  In  most  of  the  experiments,  the 
V-I  characteristics  were  of  type  B  and  could  be  measured 
clearly. 
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case 

PI 

(mmHg) 

z 

(ft) 

vc 

(KV) 

M 

<255 

(A/m) 

vb 

(V) 

1 

0.1 

3 

1.6 

6.2 

300 

0.025 

610 

2 

ft 

5 

2  .9 

9.4 

280 

0.1 

530 

3 

II 

5 

3.2 

15.5 

250 

0.25 

400 

4 

II 

3 

1.7 

7.7 

240 

0.7 

230 

5 

II 

3 

2.0 

9.3 

170 

2.5 

125 

6 

It 

3 

3.1 

25 

150 

4 

70 

7 

II 

1 

1.5 

28 

80 

20 

50 

8 

II 

1 

4.0 

31 

38 

200 

55 

9 

1 

3 

1.5 

2.5 

900 

0.0015 

340 

10 

n 

3 

2.3 

4.7 

800 

0.02 

310 

11 

II 

3 

3.0 

9.4 

480 

0.1 

265 

12 

It 

1 

1.3 

3.1 

340 

0.2 

260 

13 

II 

3 

4.2 

15.5 

260 

0.5 

207 

14 

II 

1 

1.5 

3.1 

250 

1 

135 

15 

m 

1 

2.1 

9.4 

131 

5 

75 

16 

II 

1 

2.5 

12 

120 

20 

65 

17 

ti 

1 

2.8 

14 

96 

40 

60 

18 

ii 

1 

3.5 

22 

85 

100 

95 

19 

II 

1 

4.5 

31 

68 

200 

90 

20 

10 

1 

2  .5 

1 

260 

0.00667 

175 

21 

11 

1 

2  .9 

1 

220 

0.0166 

95 

22 

ii 

1 

3.2 

2.1 

250 

0.05 

65 

23 

11 

1 

3.3 

2.1 

300 

0.1 

65 

24 

II 

1 

3.5 

2.65 

260 

0.5 

60 

25 

It 

0 

2.8 

7.2 

62 

5 

40 

26 

II 

0 

4.0 

9.35 

50 

30 

60 

(m\  A) 

II 


£.€ 

* 
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case 

n2 

nw 

X  iw 
(m) 

Xew 

(m) 

Re 

Ga 

n0 

n0 

rp 

1 

0.39 

22.3 

3 . 77xl0_6 

17.9xl0'5 

5.10 

0.887 

1 . 11x10 “2 

2 

0.56 

43.7 

1.92  " 

9.15  " 

9.47 

0.650 

1.26  " 

3 

1.02 

98.1 

0.856  " 

4.08  " 

25.8 

0.394 

1.47  " 

4 

0.46 

33.5 

2.51  " 

11.9  " 

6.57 

0.780 

1.22  " 

5 

0.56 

45.1 

1.86  " 

8.87  " 

9.18 

0.660 

1.29  " 

6 

1.47 

175.7 

0.478  " 

2.28  " 

54.1 

0.272 

1.54  " 

7 

1.51 

196.4 

0.428  " 

2.04  " 

59.2 

0.260 

1.54  " 

8 

0.40 

230 

0.365  " 

1.74  " 

5.31 

0.868 

3.06  " 

9 

3.2 

131 

0.643xl0"6 

3 . 07xl0~5 

30.2 

0.364 

2 . llxlO-2 

10 

3.8 

198 

0.424  " 

2.02  " 

46.9 

0.292 

3.11  " 

11 

5.3 

454.7 

0.185  " 

0.880  " 

85.7 

0.216 

4.07  " 

12 

3.4 

148 

0.570  " 

2.71  " 

42.7 

0.306 

1.87  " 

13 

9.0 

983 

0.0855  " 

0.407  " 

211 

0.138 

4.68  " 

14 

3.4 

148 

0.570  " 

2.71  " 

42.7 

0.306 

1.87  " 

15 

5.3 

454.7 

0.185  " 

0.880  " 

85.7 

0.216 

4.07  " 

16 

6.9 

663.4 

0.127  " 

0.603  " 

133 

0.173 

4.38  " 

17 

8.2 

852.4 

0.0985  " 

0.469  " 
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0.150 

4.58  " 

18 

12  .2 

1554 

0.0541  " 

0.257  " 

379 

0.103 

4.89  " 

19 

12  .4 

2256 

0.0372  " 

0.177  " 

453 

0.094 

4.92  " 

20 

23 

617 

0.136xl0~6 

0.649xl0"5 

157 

0.160 

2  .23x10 "2 

21 

23 

617 

0.136  " 

0.649  " 

157 

0.160 

2.23  " 

22 

31 

1200 

0.0702  " 

0.334  " 

320 

0.112 

4.66  " 

23 

31 

1200 

0.0702  " 

0.334  " 

320 

0.112 

4.66  " 

24 

33 

1380 

0.0610  " 

0.290  " 

352 

0.107 

6.16  " 

25 

41 

2982 

0.0282  " 

0.134  " 

548 

0.0854 

11.9  " 

26 

50 

4680 

0.0179  " 

0.0855  " 

767 

0.0722 

13.1  " 

9E£0 

. 

it 

: 

it 

. 
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case 

ne, 

(m-3) 

le 

(°K) 

cr 

(moh/m) 

Ep2 

(V/m) 

d2 

rp 

Epw 

(V/m) 

dw 

rp 

1 

l.llxlO16 

0.66xl04 

1.18x10 

5.30xl06 

38 

2 . 12xl07  2.3 

2 

6.15  " 

0.90 

h 

2.46  " 

5.92  " 

15.7 

4.91  " 

1.3 

3 

1.93xl017 

1.11 

ii 

3.01  " 

6.18  " 

7.0 

4.28  " 

0.53 

4 

1.67xl018 

0.84 

II 

5 . 78x  102 

4.59  " 

4.6 

1.50  " 

1.3 

5 

1 . 19xl019 

0.93 

ii 

1.42xl03 

4.48  " 

1.8 

1.24  " 

0.67 

6 

1.83  " 

1.83 

It 

0.942  " 

5.47  " 

0.68 

1.55  " 

0.235 

7 

1.75xl020 

1.50 

ii 

2.72  " 

7.75  " 

0.32 

2.01  " 

0.116 

8 

2.53xl021 

2.0 

ft 

4.32  " 

12.9  " 

0.21 

2.49  " 

0.0145 

9 

1.61xl015 

0.6 

It 

2 . 59xlO-1 

4 . OlxlO6 

20 

4.86x10 

0.44 

10 

3.42xl016 

0.6 

ii 

4.60 

4.45  " 

9.8 

6.06  " 

0.265 

11 

2 . 10xl017 

0.99 

ti 

7.65 

5.60  " 

4.2 

11.3  " 

0.13 

12 

1 . 06xl018 

0.6 

If 

1 . 35xl02 

5.10  " 

5.0 

4.15  " 

0.37 

13 

1.56  " 

1.26 

IV 

2.16x10 

7.53  " 

1.68 

10.7  " 

0.06 

14 

1.40xl019 

0.6 

If 

6.81xl02 

4.83  " 

1.88 

2.10  " 

0.36 

15 

8.93  " 

0.99 

II 

1.29xl03 

8.23  " 

0.46 

3.52  " 

0.11 

16 

5.42xl020 

1.11 

If 

2.25  " 

13.2  " 

0.24 

4.55  " 

0.072 

17 

21 

1.32xlOz 

1.2 

II 

2.85  " 

16.9  " 

0.167 

5.7 

0.052 

18 

2.49  " 

1.47 

ti 

3.57  " 

31.6  " 

0.142 

14.1  " 

0.032 

19 

4.77  " 

2.1 

II 

5.22  " 

42.7  " 

0.10 

20.8  " 

0.022 

20 

7 .30xl016 

0.6  xlO4 

1.63 

5 . 85xl06 

2.4 

10.8x10 

7  0.09 

21 

3 . 67xl017 

0.6 

II 

8.13 

3.22  " 

2.3 

5.93  " 

0.09 

22 

1.21xl018 

0.6 

II 

1.96x10 

3.16  " 

1.3 

8.43  " 

0.047 

23 

3.04  " 

0.6 

II 

4.78  " 

3.68  " 

1.08 

8.23  " 

0.047 

24 

2 . 34xl039 

0.6 

ii 

2 . 71xl02 

5.41  " 

0.61 

8.13  " 

0.040 

25 

20 

3.42xlOZU 

0.84 

II 

1 . 02xl03 

11.8  " 

0.16 

9.80  " 

0.020 

26 

2 .40xl021 

1.08 

II 

2.10  " 

28.1  " 

0.10 

27.9  " 

0.012 
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the  invalid  assumption  of  *£jsn  X.  does  affect  the  values  of  E 

E  iw  p 

significantly  and  E  gives  smaller  values  than  the  true  values.  Therefore 

pz, 

the  values  of  E  can  be  regarded  to  be  close  to  the  true  values.  There 

pw 

is  one  defect  in  calculating  E  ;  that  is,  E  has  to  be  calculated  with 

pw  pw 

a  few  more  steps  than  E^*  Consequently  E^  might  have  large  calculation 

errors.  (An  error  of  30%  in  M  gives  6%  error  to  E  .  and  40%  error  to  E  .) 

p2  pw 

In  Figs.  7.1  and  7.2,  E  _  and  E  are  shown  as  functions  of  electrode 

p2  pw 

current.  It  can  be  understood  that  E  (c^  5x10^  V/cm)  is  not  only 

pw 

sufficiently  large  for  field  emission  of  electrons  (Secs.  6.2  and  7.4) 
but  also  is  independent  of  the  values  of  and  n^  over  more  than  five 
orders  of  magnitude  as  well  as  of  the  values  of  initial  pressure  of  0.1- 
10  mmHg.  These  facts  strongly  support  the  validity  of  the  proposed  model 


of  the  breakdown  mechanism. 


7.3  POST  FIELD  ELECTRON  EMISSION  PHENOMENA 


When  the  field  emission  takes  place  at  a  particular  spot  (whose  area 
-14  -7  2 

is  usually  10  -10  cm  because  of  the  size  of  the  surface  irregularities, 

Sec.  5.3.2)  the  ions  coming  from  the  plasma  tend  to  be  concentrated  into 
that  spot  in  order  to  neutralize  the  electronic  space  charge.  Also  ions 
are  produced  near  the  area  due  to  the  ionization  of  neutral  particles  by 
field  emitted  electrons.  These  ionic  effects  increase  the  electric  field 
at  that  spot.  On  the  other  hand,  the  emitted  electrons  form  a  double 
sheath  around  the  emitted  area  (Sec.  5.1.7)  which  reduces  the  intensity 
of  the  field  at  the  area.  If,  during  the  time  electrons  are  emitted,  the 
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spot  is  heated  enough  by  the  bombardment  of  these  energetic  ions  according 
to  Eqs .  5.140  and  5.141,  then  the  thermionic  electrons  can  take  the  place 
of  the  field  electrons  even  though  the  field  electrons  are  suppressed  by 
the  double  sheath  after  that  period.  Once  a  hot  cathode  spot  is  established, 
the  stable  self  sustaining  arc  discharge  is  formed  as  mentioned  in  Sec.  5.4.3, 
and  this  arc  discharge  leads  the  electrode  system  to  a  full  path  breakdown. 


7.4  COMPARISON  OF  THE  RESULTS  WITH  EXISTING  THEORIES  AND  EXPERIMENTS 


Varney  et.  al.  (7.2)  as  well  as  Crowe  et .  al.(7.3)  carried  out  a 
one-dimensional  analysis  on  the  influence  of  the  space  charge  on  the 
spark  discharge.  In  their  calculations  (Crowe)  the  space  charge  is  so 
small  that  the  cathode  sheath  extends  much  further  than  the  anode.  This 
means  that  the  upper  limit  of  the  integral  of  the  Poisson's  equation  is 
the  anode  position,  not  the  cathode  sheath  edge.  (Sec.  7.1,  spark  discharge) 
They  have  found  that  the  influence  of  the  space  charge  on  the  threshold 
voltage  is  negligible  unless  the  initial  current  is  significantly  larger 
than  10  ^  A/cm^ . 

The  cathode  mechanism  of  extremely  short  arcs  at  atmospheric  pressure 
(Ae»d)  has  been  studied  by  Kisliuk,  Boyle  and  others  (7.4).  They 
observed  a  breakdown  voltage  as  low  as  30  volts  at  small  electrode 

O 

separations  of  the  order  of  2000  A.  The  phenomena  are  shown  to  be 
consistent  with  the  hypothesis  that  (1)  the  arcs  are  initiated  by  field 
emission  electrons,  (2)  after  the  initiation  the  anode  material  is 
evaporated  by  the  secondary  electrons  and  then  ionized,  and  (3)  the  space 
charge  thus  formed  enhances  the  electric  field  at  the  cathode,  and  increases 
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the  field  emission  current  until  the  breakdown  occurs.  Because  of  the 
steep  increase  in  current  with  field,  this  process  is  able  to  build  up 
to  breakdown  even  when  the  probability  of  any  particular  electron 
having  an  ionizing  collision  is  extremely  small. 

Their  case  is  very  similar  to  ours,  i.e.,  breakdown  is  triggered  by 
field  emission  electrons,  which  modifie  the  potential  distribution 
between  the  electrodes  and  cause  instability.  The  main  difference  is, 
of  course,  that  we  have  a  sheath  edge  instead  of  their  anode.  Also  the 
vaporization  of  anode  material  (this  makes  it  difficult  to  calculate  the 
neutral  particle  density  between  the  electrodes)  is  substituted  for  the 
dense  neutral  gas  in  the  boundary  layer  of  the  stagnated  shock. 

The  modification  of  their  breakdown  theory  to  fit  our  case  will  be 
presented  briefly.  Consider  the  field  emission  equation  (Eq .  5.111), 

jG-4E2e'^,  (7-0 

where  A  and  B  are  constants.  Let  E=Ep+E ' ,  where  Ep  is  the  initial  field 
and  E'  is  the  field  due  to  the  additional  positive  ion  space  charge  (ions 
created  by  ionizing  collisions  of  field  emission  electrons)  measured  at 
the  cathode  surface.  Because  of  the  steepness  of  the  dependence  of 
current  on  field,  E'^\Ep  and  by  expansion  we  find, 

_Jt  M 

je  =  aeV  6p  e  v .  (7-2) 


From  the  Poisson's  equation. 


n.e 


Ji  +  J' 


dr 


^(r2ePe')  =  re^, 


(7.3) 
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where  E^^=J  72**0^  ,  is  the  initial  ion  current  per  unit  length  of 
the  cathode  and  J'  is  the  extra  positive  ion  current  arising  from  . 


Since  the  change  in  sheath  thickness  is  negligible, 


c«2 
^ 


2£P 


&->) 


(7.5) 


or  at  the  cathode  surface, 

p*2. 


(7.6) 


where  C=^rg/r  y  -iy41C£0(4  ^Ep.  The  breakdown  condition  is 


Je(e 


n0Q,-(^-rp) 


-i)  -  j'. 


From  Eqs .  7.2,  7.6  and  7.7,  we  obtain, 

®  §£(  n,(?i(rs-rp) 

2TcrpAEP  e  sf  e  ^  (6 


(7.7) 


(7.8) 


Breakdown  occurs  when  the  plots  of  the  right  and  left  hand  side  of  Eq.  7.8 
are  tangent,  i.e.. 


ABrPe(e 


i  t'o^i  £p 


(7.9) 


Equation  7.9  is  a  very  strong  function  of  Ep  and  the  exact  value  of  Ep 
required.  In  our  experiment,  however,  Ep  is  obtained  after  lengthy 
calculations  and  might  have  a  considerable  amount  of  error.  In  addition, 
this  theory  neglects  the  double  sheath  effect.  Therefore  it  does  not  seem 
worth  while  to  test  our  results  with  the  condition  of  breakdown  of  Eq .  7.9. 

The  effect  of  temperature  in  the  range  of  room  temperature  to  the 
neighborhood  of  1000°C  was  studied  by  Bowder  (7.5)  and  Alston  (7.6).  In 
this  temperature  range,  if  the  gas  density  is  kept  constant  and  if  the 
electrical  field  is  uniform,  the  breakdown  voltage  was  rather  insensitive 
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to  the  gas  temperature. 

In  1961  and  1964,  Sharbaugh,  Lee  and  others  (7.7)  measured  the 
breakdown  voltage  of  a  plasma  produced  by  a  pressure  driven  shock  tube, 
and  observed  a  breakdown  voltage  as  low  as  50%  (at  4000°K)  of  that 

predicted  by  Paschen's  Law.  Exact  theoretical  analysis  has  not  been  made. 

10  -3 

The  electron  density  was  the  order  of  10  cm  which  indicated  that  the 
thickness  of  the  cathode  sheath  across  which  most  of  the  potential 
difference  was  applied  became  comparable  to  the  hydrodynamic  boundary 
layer  formed  along  the  shock  tube  wall.  However,  the  neutral  particle 
density  (to  which  first  ionization  coefficient  is  proportional)  inside 
and  outside  the  boundary  layer  was  not  known. 

Lee  et.  al.  also  have  concluded  that  for  a  high  density  plasma 
10  -3 

(.>2x10  cm-J)  ,  thermal  breakdown  takes  place  instead  of  dielectric 
breakdown,  i.e.,  after  the  initiation  of  breakdown  and  while  the  breakdown 
process  is  developing,  a  great  amount  of  energy  (  -  100  KW)  is  supplied  to 
the  gas  and  further  ionization  takes  place  which  leads  to  the  complete 
full  path  breakdown. 

The  electric  field  at  the  cathode  surface  of  -  2x10^  V/cm  may  seem 
a  bit  low  for  the  field  emission.  This  field  strength,  however,  is 
enough  as  supported  by  existing  experimental  results  by  many  people  other 
than  the  results  obtained  in  Sec.  6.2.  The  value  of  electric  field  of 
the  order  of  10^-10^  V/cm  was  reported  by  Llewellyn  Jones  and  Morgan  (7.8) 
and  Morgan  and  Harcombe  (7.9)  in  studies  of  enhanced  electron  emission 
from  cold  electrodes  in  gases,  in  agreement  with  the  estimates  made  on  the 
basis  of  the  Fowler -Nordheim  equation  by  Kerner  and  Raether  (7.10)  from 
their  measurements  of  electron  emission  from  cold  metals  under  electric 
fields  of  the  order  of  10^  V/cm.  Calvert  (7.11)  also  reported  a  cold 
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emission  of  electrons  under  electric  fields  of  3-8x10^  V/cm.  His 
estimated  emitting  area  of  the  electrons  was  10_i^  cm  . Penning  and  Mulder 

(7.12)  and  Lamar  and  Compton  (7.13)  observed  the  field  emission  at 
2.5x10^  and  7.6x10^  V/cm  or  more,  respectively. 


(7.2)  R.N.  Varney,  H.J.  White,  L.B.  Loeb  &  D.Q.  Posomin:  Phys .  Rev.  48^ 

818  (1935) 

(7.3)  R.W.  Crowe,  J.K.  Bragg  &  V.G.  Thomas:  Phys.  Rev.  96.  10  (1954) 

(7.4)  L.H.  Germer  &  F.E.  Haworth:  Phys.  Rev.  7_3  1121  (1948) 

L.H.  Germer  &  F.E.  Haworth:  J.  Appl.  Phys.  20_  1085  (1949) 

L.H.  Germer:  J.  Appl.  Phys.  22.  955  (1951) 

L.H.  Germer  &  J.L.  Smith:  J.  Appl  Phys.  2J3  553  (1952) 

W.P.  Dyke  &  J.K.  Troian:  Phys.  Rev.  89  799  (1953) 

P.  Kisliuk:  Phys.  Rev.  92.  847  (1953) 

P.  Kisliuk:  J.  Appl.  Phys.  2_5  897  (1954) 

W.S.  Boyle  &  P.  Kisliuk:  Phys.  Rev.  97.255  (1955) 

W.S.  Boyle,  P.  Kisliuk  &  L.H.  Germer:  J.  Appl.  Phys.  2_6  720  (1955) 
W.S.  Boyle  &  F.E.  Haworth:  Phys.  Rev.  101  935  (1956) 

R.N.  Bloomer  &  B.M.  Cox:  Vacuum  L8  379  (1968) 

(7.5)  H.C.  Bowder:  Proc.  Phys.  Soc .  43  96  (1931) 

(7.6)  L.L.  Alston:  Proc.  IEE  (london)  105  pt.  A  549  (1958) 

(7.7)  A.H.  Sharbaugh,  P.K.  Watson,  D.R.  White,  T.H.  Lee  &  A.  Greenwood: 

AIEE  Trans,  pt .  3  PAS  80  333  (1961) 

T.H.  Lee,  A. N.  Greenwood  &  D.R.  White:  IEEE  Trans.  PAS  84  1116  (1965) 

(7.8)  F.  Llewellyn  Jones  &  C.G.  Morgan:  Proc.  Roy.  Soc.  A2 18  88  (1953) 

(7.9)  C.G.  Morgan  &  D.  Harcombe:  Proc.  Phys.  Soc.  B66  665  (1953) 

(7.10)  K.  Kerner  &  H.  Raether:  Z.  angew.  Phys.  6.212  (1954) 

(7.11)  W.J.R.  Calvert:  Proc.  Phys.  Soc.  B69  651  (1956) 

(7.12)  F.M.  Penning  &  J.G.W.  Mulder:  Physica  2_  724  (1935) 

(7.13)  E.S.  Lamar  &  K.T.  Compton:  Phys.  Rev.  yi_  1069  (1931) 


1 .  )  *.i>  1 1  f  i  •  ;  - 


X 


•  •  '  .  ,  .  •  •} xiW  . L.H  (£.$) 

.  .  •  :  .  ,>.L  .  (  .  \) 

... 

.  .  •  •  •'  .  . 

.  I  :  ■ 

•  /  .  ■  "/>  •  1  :  '  .■  •  Sr.  .h.J 

:  r  .L  '  , 

:  .  r;  . 

■  .t  :  ... 

,  -  ’* 

:  .  -  v  . 

.  .  .  •  i  i  :i 

'>  ' 

•  '  -  :  .  .  .  :  . 

•  •  ’  1  •  •  •  ...  ,  [,  •  (•  .  v) 

•  .  .  ...  • 

•  ’•  • 

•  ”  •  .  .  .  ;  .  v) 

:  .  ,J  •  ;  '  .  • 

('  r  .  •.  •  -■  .  .  , t  ,  i  (n . 

:  .  .r  ,i  3  •  .  ,  ;  ; 

exriS  uiclqtaoD  .T,a  £>  ijatnaJ  .8.3  (? 


188 


CHAPTER  VIII  CONCLUSION 


The  electrical  breakdown  phenomenon  between  two  electrodes  changes 
drastically  when  there  is  a  plasma  surrounding  these  electrodes.  One  of 
the  results  of  this  change  is  the  lowering  of  the  breakdown  voltage  (it 
decreased  to  40  volts  or  less  which  was  much  less  than  Paschen’s  minimum 
breakdown  voltage  (Secs.  6.1  and  6.6)).  It  was  assumed  that  this  breakdown 
phenomenon  was  due  to  electrode  surface  contamination.  Experimental  work 
proved  that  this  assumption  was  correct  (Sec.  6.3).  If  the  total  charge 
accumulated  on  the  electrically  insulating  surface  layer  reached  a 
certain  amount,  then  the  field  in  the  layer  became  intense  enough  to  cause 
the  electron  avalanche,  which  led  the  system  to  the  full  path  breakdown. 
Therefore  the  breakdown  condition  was  that  the  integral  of  the  pre-breakdown 
electrode  current  from  time  zero  to  the  time  when  breakdown  occured  was  some 
constant  value  and  was  independent  of  the  applied  voltage. 

Consequently  if  the  electrode  surface  had  been  clean  the  breakdown 
voltage  should  have  been  higher  than  Paschen's  minimum  value.  Therefore 

f 

the  electrode  was  cleaned  well  by  heating  to  2300°C  (Sec.  5.3.7)  and  the 
measurements  were  carried  out;  the  breakdown  occured  at  40  volts  again. 

There  must  be  some  other  mechanism  to  cause  the  breakdown.  The  next 
hypothesis  was  that  this  phenomenon  was  due  to  the  hot  cathode  arc 
discharge.  The  plasma  density  (which  is  roughly  proportional  to  electrode 
current)  against  breakdown  voltage  curve  of  Fig.  6.12  nicely  coincided 
with  the  well  known  V-I  characteristic  of  the  breakdown,  especially  with 
the  transition  region  from  abnormal  glow  discharge  to  arc  discharge.  If 
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this  was  the  case  then  the  negative  electrode  (or  cathode)  should  have 
been  heated  by  the  ion  bombardment  to  a  temperature  high  enough  for 
thermionic  electron  emission.  In  some  cases  the  energy  input  rate,  I-^Vp , 
into  the  0.1  mm  diameter  cylindrical  electrode  reached  1000  watts,  which 
appeared  abnormally  large.  After  some  calculations,  however,  it  was  found 
that  the  energy  loss  by  heat  conduction  into  the  metal  was  so  large  that 
the  cathode  temperature  rise  was  only  10°C  at  most  (Sec.  5.4.5).  Thus 
the  hot  cathode  arc  hypothesis  was  proved  invalid. 

Many  effects  happening  in  the  sheath  region  such  as  current  pinching, 
oxide  surface  layer,  local  heating  by  ion  bombardments,  etc.  were  tested 
(Secs.  5. 1.6, 5. 1.8,  5.1.9,  5.3.9,  5.4.6,  5.4.7)  to  explain  this  peculiar 
breakdown  phenomena,  in  vain.  The  last  hypothesis  was  to  assume  cold- 
cathode  arc  where  the  field  emission  of  electrons  from  the  cathode 
dominates  the  situation.  From  rough  calculations,  the  intensity  of  the 
field  at  the  cathode  surface  was  only  approximately  10^  V/cm,  which  was 
too  low  for  field  emission.  However,  after  careful  examination  of  the 
neutral  density  change  in  the  incident  and  stagnated  shock  as  well  as  in 
the  boundary  layer  around  the  electrode,  the  ionic  mobility  was  found  to 
be  considerably  lower  than  that  used  in  the  rough  calculations ,  and  this 
increased  the  field  intensity  by  approximately  ten  times.  Also  this  increase 
in  neutral  particle  density  enhanced  the  ionizing  collisions  of  neutrals 
by  secondary  electrons  and  brought  the  resultant  field  to  about  5x10"*  V/cm. 
This  intensity  could  be  considered  to  be  strong  enough  for  field  emission 
by  looking  into  the  supplementary  experiments  (Sec.  6.2)  and  many  reports 
of  vacuum  breakdown  experiments  (Sec.  7.4).  In  addition,  this  electric 
field  strength  at  the  cathode  was  constant  over  a  wide  range  of  plasma 
density  (10^-10^  cm”^)  and  initial  pressure  (0.1-10  mmHg) .  Thus  the 
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breakdown  mechanism  was  finally  made  clear. 

The  application  of  the  theory  of  the  breakdown  mechanism  developed 
in  this  thesis  is  virtually  limitless.  The  most  important  application 
arises,  as  mentioned  in  Chapter  I,  to  the  reignition  phenomena  after 
current  zero  of  an  alternating-current  circuit  breaker  where  the  electrodes 
are  surrounded  by  the  high  density  post -arc  plasma.  Also  this  theory  is 
applicable  to  the  mechanisms  of  the  unipolar  arc  discharge  as  well  as 
the  vacuum  breakdown  and  the  analysis  of  electrostatic  probe  current  in 
a  dense  plasma. 

During  the  course  of  this  research  work  many  other  valuable  results 
were  also  obtained.  One  of  them  is  the  development  of  a  new  technique  of 
determining  the  electrical  conductivity  of  plasma  by  using  radio  frequency 
magnetic  fields  (Sec.  3.2.5).  In  section  4.2  a  new  theory  which  predicts 
the  variation  of  the  plasma  velocity  is  introduced,  with  good  agreement 
with  the  experimentally  obtained  values  (Sec.  4.3.1).  Also  in  sections 
4. 3. 2 -4. 3. 7  many  properties  of  the  plasma  produced  by  the  electromagnetically 
driven  shock  tube  were  carefully  examined  by  many  diagnostic  techniques 
and  the  validity  of  the  shock  relations  was  tested. 

Section  2.7  describes  a  generalized  floating  double  probe  theory  as 
well  as  an  interesting  application  of  the  double  probe  to  measure  a  low 
density  plasma  (-10'  cm  ) .  The  ion  current  to  the  probe  which  was 
aligned  parallel  to  the  plasma  flow  was  measured  in  the  case  of  continuum 
high  speed  plasma  flow  and  an  explanatory  theory  was  developed  (Sec.  2.6.4). 
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